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g-RANDOMIZED ROBINSON-SCHENSTED-KNUTH CORRESPONDENCES 

AND RANDOM POLYMERS 

KONSTANTIN MATVEEV AND LEONID PETROV 


Abstract. We introduce and study g-randomized Robinson-Schensted-Knuth (RSK) corre¬ 
spondences which interpolate between the classical (g = 0) and geometric (g 1) RSK cor¬ 
respondences (the latter ones are sometimes also called tropical). 

For 0 < g < 1 our correspondences are randomized, i.e., the result of an insertion is a certain 
probability distribution on semistandard Young tableaux. Because of this randomness, we use 
the language of discrete time Markov dynamics on two-dimensional interlacing particle arrays 
(these arrays are in a natural bijection with semistandard tableaux). Our dynamics act nicely on 
a certain class of probability measures on arrays, namely, on g-Whittaker processes (which are 
t = 0 versions of Macdonald processes of Borodin-Corwin El). We present four Markov dynamics 
which for g = 0 reduce to the classical row or column RSK correspondences applied to a random 
input matrix with independent geometric or Bernoulli entries. 

Our new two-dimensional discrete time dynamics generalize and extend several known con¬ 
structions: (1) The discrete time g-TASEPs studied by Borodin-Corwin [S] arise as one-dimen¬ 
sional marginals of our “column” dynamics. In a similar way, our “row” dynamics lead to discrete 
time g-PushTASEPs — new integrable particle systems in the Kardar-Parisi-Zhang universality 
class. We employ these new one-dimensional discrete time systems to establish a Fredholm deter- 
minantal formula for the two-sided continuous time g-PushASEP conjectured by Corwin-Petrov 
[23]. (2) In a certain Poisson-type limit (from discrete to continuous time), our two-dimensional 
dynamics reduce to the g-randomized column and row Robinson-Schensted correspondences in¬ 
troduced by O’Connell-Pei [59] and Borodin-Petrov [16], respectively. (3) In a scaling limit as 
g ^ 1, two of our four dynamics on interlacing arrays turn into the geometric RSK correspon¬ 
dences associated with log-Gamma (introduced by Seppalainen [7D]) or strict-weak (introduced 
independently by O’Connell-Ortmann [58] and Corwin-Seppalainen-Shen [25]) directed random 
lattice polymers. 
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1. Introduction 


1.1. Overview. The classical Robinson-Schensted-Knuth (RSK) correspondence associates to 
an integer matrix a pair of semistandard Young tableaux of the same shape [l5j, [35], nn, m- 
It is informative to view an integer matrix M = (Mij) as a configuration of points (“balls”) in 
cells of the lattice with Mij balls in the (i, j)-th cell (see Fig.lll left). There are also simpler 
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Figure 1. Left: an integer matrix as an input to the RSK. Center: an integer word 
as an input into the RS viewed as a matrix with the continuous j coordinate (at most 
one ball at a given horizontal position is allowed; the word encodes vertical positions 
of consecutive balls). Right: a permutation viewed as a matrix with both continuous 
coordinates (at most one ball at a given horizontal or vertical position is allowed; balls 
represent the graph of the permutation). 

correspondences obtained from the RSK if one makes one or both dimensions of the input con¬ 
tinuous, see Fig. center and right. In particular, the Robinson-Schensted (RS) correspondence 
maps integer words into pairs of Young tableaux of the same shape, but now one of them is 
standard. 

The idea of applying the RSK correspondence to a random input can be traced back to m 
where it was used in connection with the asymptotic theory of characters of the infinite symmetric 
group (see also m)- Together with combinatorial properties of the RSK this idea has been 
extensively employed in studying various stochastic systems, e.g., TASEP (= totally asymmetric 
simple exclusion process), the last-passage percolation |l2| , or longest increasing subsequences of 
random permutations m, m- 

Reading the random input matrix column by column adds a dynamical perspective to random 
systems (with j in all three cases on Fig. [^playing the role of time). This direction has been 
substantially developed in, e.g., IS], m, 0 . 

The geometric version (also sometimes called ^HropicaF) of the RS and the RSK correspon¬ 
dence^ has also been employed in the study of stochastic systems [56] , [22] , [60] , [57] . The 
systems one obtains at this level are related to directed random polymers in random media, in 
particular, to the O’Connell-Yor, log-Gamma, and strict-weak random polymers introduced in 
m, and |2H], 122], respectively. Each such polymer model can be viewed as a positive 
temperature version of a certain last-passage percolation-like model. 

In the stochastic systems mentioned above, the RSK and related constructions provide a way 
to observe and understand their integrability. The integrability property refers to the presence of 

^The geometric RSK maps arrays of positive real numbers into other such arrays in a birational way, and is 
obtained from the classical RSK by a certain “detropicalization”, see |44) . [53] . 
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concise and exact formulas describing observables, which allows to study the asymptotic behavior 
of such systems, and also gives access to exact descriptions of limiting universal distributions, 
such as the Tracy-Widom distributions which are features of the Kardar-Parisi-Zhang (KPZ) 
universality class m, m, US] m- 

The classical RSK is deeply connected to Schur symmetric functions |52l Ch. I], while the 
geometric RSK is relevant to the Whittaker functions [l9], [28]. Both families of functions 
are degenerations of more general Macdonald symmetric functions depending on two parameters 
{q,t) |52l Ch. VI]: the Schur functions correspond to q = t, and the Whittaker functions arise in 
the limit as t = 0 and q 1, |38j . 

In the recent years, there has been a progress in understanding analogues of the RS corre¬ 
spondences at other levels of the Macdonald hierarchy: g-Whittaker {t = 0 and 0 < g < 1) [59] . 
[65], [T6] and HalKLittlewood ((7 = 0 and 0 < t < 1) [18]. At these levels, the correspondences 
become randomized, that is, the image of a deterministic word (as on Fig. center) is no longer 
a fixed pair of Young tableaux, but rather a random such pair. Because of this randomness, an 
appropriate language for describing the correspondences seems to be that of Markov dynamics 
on two-dimensional interlacing integer arrays (these arrays are in a natural bijection with semi¬ 
standard tableaux, see Remark o below for more detail). The dynamics which are analogues 
of the RS correspondences evolve in continuous time according to the j axis on Fig. [^ center. 
These dynamics act nicely on certain families of probability distributions on interlacing arrays, 
namely, the Macdonald processes [7] . 

The g-Whittaker level is relevant to integrable one-dimensional particle systems such as (con¬ 
tinuous time) gr-TASEP and the stochastic g-Boson system [69]. [7]. [12], [TO], [29], and (contin¬ 
uous time) (7-PushTASEP (= (/-deformed pushing TASEP) [23] q In particular, continuous time 
Markov dynamics on interlacing arrays constructed in [59] and [16] are two-dimensional exten¬ 
sions of, respectively, the g-TASEP and the g-PushTASEP. That is, the latter one-dimensional 
processes are Markovian marginals of the dynamics on two-dimensional interlacing arrays]^ 


In the present paper we advance further at the (/-Whittaker level, and introduce four q- 
randomized RSK correspondences, or, in other words, four discrete time Markov dynamics on 
interlacing arrays which act nicely on g-Whittaker processes (these are Macdonald processes with 
t = 0). These dynamics unify, generalize and extend all of the above RSK-type constructions: 


• When (/ = 0, our four g-randomized correspondences become usual or dual, row or column 
classical RSKs (four classical correspondences in total). The input matrix M in the usual RSKs 
has Mij G {0,1,2,...}, and in the dual RSKs one has Mij G {0,1}. When one takes Mij to be 
independent geometric (for usual) or Bernoulli (for dual) random variables and applies a suitable 
classical RSK, the shape of the resulting random Young diagram is distributed according to 
the Schur measure [63] Jj Similarly, our (/-randomized correspondences applied to (/-geometric or 
Bernoulli random inputs (note that the Bernoulli input needs not to be (/-deformed) give rise to 
(/-Whittaker distributed random Young diagrams. The latter property is an instance of “acting 
nicely” on (/-Whittaker processes (see also (2.17) in ^for more detail). 


^These systems are in fact quantum integrable in the sense of the coordinate Bethe ansatz in, m, 0. [in]. 

^The two-dimensional dynamics at the Hall-Littlewood level [18], however, do not seem to lead to any new 
one-dimensional integrable particle systems. 

^In the present paper, the word “geometric” is attached to two separate concepts — the geometric RSKs, 
and the geometric and q-geometric random variables. To avoid confusion where it can occur, we will call the 
correspondences the geometric (tropical) RSKs. See also Remark I 


8.10 
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• In a limit from discrete to continuous time, our g-randomized RSKs turn into the (simpler) 
(^-randomized RS correspondences introduced and studied in [59] . m- 

• The two discrete time (;-TASEPs (associated with (^-geometric or Bernoulli random vari¬ 
ables) studied by Borodin-Corwin [8] arise as one-dimensional marginals of our two “column” 
dynamics on interlacing arrays. In a similar way, our two “row” dynamics lead to discrete time 
(^-PushTASEPs — new integrable particle systems in the KPZ universality class. 

• In a scaling limit as (? /' 1, the dynamics on interlacing arrays associated with the g-geometric 
random input (these are two out of our four g-randomized RSK correspondences) converge to 
geometric (tropical) RSK correspondences. The latter correspondences (which are deterministic 
birational maps between arrays of positive reals) are relevant to the log-Gamma |70j . [22], |60j 
and strict-weak |25| random lattice polymers. 

In §1.2| below we describe one of our four dynamics in detail, and in §1.3| we briefly discuss 
other dynamics and results. 
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Eigure 2. Left: An interlacing array A; we require that A) ^ G '^>o- Right: A config¬ 
uration of particles corresponding to an interlacing array of depth N = 5 (right). 


1.2. g-randomized row insertion with (/-geometric input. Discrete time Markov dynamics 
(i.e., the (/-randomized RSK correspondences) which we construct in the present paper live on the 
space of integer arrays A (see Eig. [^. Neighboring levels of the array satisfy certain inequalities 
which we call the interlacing property (see (2.1) for the definition). Each level = (A^'^^ > 
... > A^^) of an array can be viewed as a partition (equivalently, a Young diagram |52l I.l]), so 
A is a sequence of interlacing Young diagrams. 
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Eigure 3. A semistandard Young tableau corresponding to the array on Eig. right. 

Remark 1.1. Each A can be also viewed as a semistandard Young tableau of shape filled 
with numbers from 1 to N. Here “semistandard” means that in this filling, numbers increase 
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weakly along rows, and strictly down columns. Then each is the portion of the semistandard 
tableau hlled with numbers from 1 to j, see Fig. 

Let us now define the (g-randomized) operation of inserting a word w = 

(i.e., the word has mi ones, m 2 twos, etc.) into an array A. The result is a new, random array u. 
At the first level we have +mi. Then, sequentially at all levels j = 2,..., N, given the 

existing change —)• at the previous level and the old state A^-^^ at the current level, 

construct the new state as follows. Each move — Ap f = 1, ..., j — 1, is randomly 

split into two pieces r) ^ Y and the piece r) ' is added to the new move of the upper 

right neighbor while the piece is added to the new move of the upper left neighbor 

(i) f?) 

A);_(_P Moreover, A) ^ receives an additional move of size mj. All these splittings and moves at 
level j happen in parallel. That is (here and below 1... stands for the indicator). 




(i) 

- \ 
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‘-i<j 


^0-1)] 
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U>1, 




(see Fig.|^. To complete the definition, it now remains to describe the distribution of the splitting 


of the move y 


0 - 1 ) 


- A, 


0-1) _ .^0-1) , tfO-i) 


= rr ' + iy ■ This is a certain g'-deformed version of the Beta- 


binomial distribution, namely, r, 
Ap with probability 
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. (g^-";g-i), (g°;g-^)c-r- (g'^g"^) 


( 1 . 1 ) 


where 


» = A® - aA‘>, 


^m—1 




{u]q)m = (1 — u){l — uq )... (1 — uq^ ^) are the g-Pochhammer symbols, and we adopt the 


convention Aq = -|-oo. The quantity 


is simply equal to — A^ 


The quantities (1.1) define a probability distribution in r for a < 6, c < 6 (these conditions 


follow from the interlacing). Moreover, this distribution is supported on {0,1,..., c} n {c — a, c — 
a-|- 1,...,6 — a — 1,6 — a}, which in fact ensures that the new array v is also interlacing (see 


lemma 6.2 for details). 
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Figure 4. Splitting of the move at level j — 1 and its propagation to the level j. Here 
we are using the particle interpretation of interlacing arrays as on Fig. right. 


Remark 1.2. The interlacing array A plays the role of the insertion tableau in our g-randomized 
RSK correspondence (cf. Remark 1.1). One can readily define an accompanying recording tableau 
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in the same way as it is done for the classical RSK correspondences. In the present paper we will 
not focus on recording tableaux. 

Now, let us take the insertion word w = ^ ^ to be random itself. More precisely, 

let rrij, j = 1,..., N, he independent g-geometric random variables: 

A; = 0,1,2,..., 0<a<l. (1.2) 


a 


Prob(mj = k) = - - :-{a;q)oc, 

[q;q)k 

Inserting this random word w into an array A dehnes one step of a discrete time Markov chain 
on interlacing arrays. We denote this Markov chain by Qrow[«]- 

Theorem 1.3. Start the Markov dynamics Qrow[Q:] from the interlaeing array with all Ap^(O) = 0. 
Then the distribution of the array A(T) after T steps of this dynamics is given by the q-Whittaker 
process: 

Prob(A(r) = A) = {a;q)l 


\TN 

)co 




Here Px/fi and Qx are the g-Whittaker polynomials, see ^for more detail. Theorem 1.3 follows 
from Theorem |6.4| which we prove in 


Remark 1.4. In fact, we can (and will) consider a more general situation when the parameters 
a in (1.2) may depend on j and on the time step as ataj. Then the g-Whittaker process above 
takes the form 

NT . 

g)oo jPA(i)(®i)^A( 2 )/A(i) (“ 2 ) • • • PxW/x(N-^)i(^N)Qxwio:i^ “ 2 , • • • ,ar)- 

We omit the dependence on f and t in Introduction. 

Let us now describe three degenerations of the dynamics Qrow[a]: 

• For q = 0, the splitting distributions (1.1) become supported at a single r G {0,1,..., c}, so 
the randomness in the insertion disappears, and the insertion itself turns into the classical RSK 
row insertion (we recall its definition in 14.3). The g-geometric random variables mj (1.2) become 
geometric, and the ( 7 -Whittaker polynomials in Theorem 1.3 turn into the Schur polynomials. This 
justifies our treatment of the dynamics Qrow[a] as the g-randomized row RSK correspondences. 

• Fix 0 < (7 < 1. When a 0 in (1.2) and one rescales time from discrete to continuous time 
(see 16.7 for details on this scaling), the random input matrix turns into N independent Poisson 
processes running in parallel (i.e., we are passing from the left to the center situation on Fig. [^. 
Then in the splitting distributions one has c = 0 or 1, and the dynamics Qrow[Q:] turns into a 
continuous time dynamics on g-Whittaker processes which was introduced in m- The latter 
continuous time dynamics should be viewed as a g-randomized row RS correspondence. 

• Let q = e~^ and a = with e 0 and 0 > 0. Define the positive random variables 
fli(i.f) via the following scaling: 


^(^) = (^ + J - 2/c + l)e ^loge ^ + e ^ log (.R;^(t, e)). 


If the quantities evolve under the dynamics Qrow[a] started from all A^^^(O) = 0, then the 
rescaled quantities W^{t,e) converge to certain ratios of partition functions in the log-Gamma 
lattice polymer model (see ^8.1| and Theorem 8.7 in particular for details). Moreover, under this 


scaling the randomness in the splitting (1.1) disappears, and the g-randomized insertion described 
above turns into the geometric RSK insertion. 
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Remark 1.5. It is worth noting that there is also a strong connection between the geometric 
(tropical) RSK and representation theory, cf. [5], |19j . At the g-randomized level this connection 
does not (yet) seem to be present. 


(i) 

When restricted to the rightmost particles A) , j = of the interlacing array, the 

dynamics Qrowfct] induces a marginally Markovian evolution which we call the (discrete time) 
geometric g'-PushTASEP. This is a new integrable particle system in the KPZ universality class. 
In the shifted coordinates Xi{t) := — i (so xn < ... < xi), the evolution of this system 

during time step t —)■ t + 1 looks as follows. Sequentially for i = 1,2,..., A", each particle Xi 
jumps to the left by mj + r) , where rrii is an independent g-geometric random variable (1.2), 
and is a random variable with distribution 


I c) 


g-(g“;g 


)c—r 


(q 


- 1 . 


(g ^)c 
g ^)r(g g ^)c—r 


a = Xi-i{t) - Xi{t) - 1, 

C = Xi-i{t + 1) - Xi-i{t) 


(this is simply the splitting distribution 0 with b = +oo). Note that if c > a, then 
chosen according to the above distribution will be at least c — a. See Fig. 


a = Xi-i{t) - Xi{t) -1 = 7 

j_i) c = Xi-i{t + 1) - Xi-i{t) = 4 



Xi{t) Xi-l{t + l) Xi-l{t) 


Figure 5. The discrete time geometric g-PushTASEP. 


In a continuous time limit as a \ 0, the geometric g-PushTASEP turns into the continuous 
time g-PushTASEP of [16], |23j. The g-moments of the form E g^(and more general 
such moments) of both g-PushTASEPs are given in terms of nested contour integrals. For the 
geometric g-PushTASEP only finitely many such moments exist (i.e., the expectation is infinite 
for sufficiently large fc), and for the continuous time g-PushTASEP the moments grow too fast 
and also do not determine the distribution of Xn{t). However, it is still possible to write down 
a Fredholm determinantal formula for the distribution of Xn{t) for both processes (started from 
the step initial configuration Xi(0) = —i) using the theory of Macdonald processes [7|, see jQ] 
Theorem 3.3]. We refer to ^for further details. 


[a] discussed in 11.2 above, we 


1.3. Other dynamics and results. Besides the dynamics Qm 
introduce three other dynamics on g-Whittaker processes: 

• QpoJct] ((6.4 and Theorem 6.11). At g = 0 this dynamics becomes the classical RSK column 
insertion applied to a geometric random input ((4.3). In a scaling limit as g 1, 

turns into a geom etric (tropical) RSK associated with the strict-weak lattice polymer introduced 
in |25| (Theorem 8.8). In a continuous time limit, turns into the g-randomized column 

RS correspondence introduced in |59|. Under the leftmost particles X^p of the interlacing 

array evolve according to the discrete time geometric g-TASEP of |8|. 
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• Q?ow[/3] (^5.1 and Theorem 5.2). At (? = 0 th is dy namics becomes the dual RSK row insertion 
applied to a Bernoulli random input Q?c 5 w[/ 3 ] (?4-3). In a continuous time limit, Qrow[P] turns 
into the g-randomized row RS correspondence m- Under Qtow[/3], the rightmost particles 
of the interlacing array evolve according to a new particle system, the discrete time Bernoulli 
( 7 -PushTASE P (D efinition 7.1). 


aJoiM (>5.4 


and Theorem 


5.7). At q = 0 this dynamics becomes the dual RSK column 
insertion applied to a Bernoulli random input (^4.3). In a continuous time limit, Qrow[/3] 

turns into the q-randomized column RS correspondence of [59]. Under Qrow[/3], the^leftmost 

f?) I—I 

particles of the array evolve according to the discrete time Bernoulli q-TASEP | 8 ] 


Remark 1.6. We do not attempt a full classification of g-randomized RSK correspondences as it 
was done for the RS correspondences by solving certain linear equations for transition probabilities 
in |16j . Similar equations for the discrete time situation seem to be much more involved, and 
in this paper we demonstrate particular solutions to these equations which lead to discrete time 
Markov dynamics (see also 12.6.1 for further discussion). 

We however believe that the four dynamics we construct are the most “natural” discrete time 
dynamics on g-Whittaker processes having all the desired properties and prescribed degenerations: 


• The update in the dynamics is sequential, from lower to upper levels of the interlacing 
array. 

• The dynamics act nicely on g-Whittaker measures and processes. 

• The continuous time limits (a or /3 —)• 0) of the dynamics coincide with continuous time 
RS dynamics of or m- 

• For q = 0, the dynamics degenerate to the ones related to the classical RSK correspon¬ 
dences. 

• In the q 1 limit, the (a) dynamics converge to the ones related to the geometric 
(tropical) RSKs. 


The dynamics Q™ 


and [/3] are related to each other via a straightforward procedure we 


call complementation (15.3) which shortens the proofs for Moreover, one can say that this 

procedure provides a direct link between the column and the row g-randomized RS correspon¬ 
dences of |59| and |T6| (which are continuous time limits of [/3] and Qmw [/3], respectively). This 
also provides a direct coupling between the Bernoulli g-TASEP and ( 7 -PushTASEP (Proposition 


7.3). 


We employ the discrete time Bernoulli processes to obtain a Fredholm determinantal formula 
for the continuous time g-PushTASEP and for its two-sided extension, the continuous time q- 
PushASEP (the latter formula was conjectured in |23|), see Theorem 7.10 See also a related 


discussion in the end of 11.2 


1.4. Outline. In ^we recall the necessary background on Macdonald and g-Whittaker symmet¬ 
ric functions and g-Whittaker processes, and also write down and discuss main linear equations 
which must be satisfied by our Markov dynamics on interlacing arrays. In ^ we discuss two 
particular types of Markov dynamics, namely, push-block and RSK-type dynamics, and explain 
the differences between them. In ^we illustrate our main definitions and concepts in the q = 0 
situation, when the g-Whittaker polynomials reduce to the simpler Schur polynomials, and the 
dynamics on interlacing arrays are relevant to the classical RSK correspondences. In ^and ^we 


®In contrast with Q?ow[«] and Q|?oi[a], there is (yet) no known polymer-like limits of Q?ow[d] or 
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explain in detail the constructions of four discrete time RSK-type dynamics on interlacing arrays, 
and prove that these dynamics act on the g-Whittaker processes in desired ways. In ^we dis¬ 
cuss moment and Fredholm determinantal formulas for our one-dimensional interacting particle 
systems. In ^we consider scaling limits as g 1 of our two (a) dynamics on interlacing arrays, 
and show that they turn into the geometric RSK correspondences associated with log-Gamma or 
strict-weak directed random lattice polymers. 
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for numerous discussions which were extremely helpful. LP would like to thank Sergey Fomin, 
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2. Macdonald processes and Markov dynamics 

In this section we collect main notation and definitions related to Macdonald processes used 
throughout the paper, and also write down and discuss linear equations satisfied by Markov 
dynamics on g-Whittaker processes which we aim to construct. 

2.1. Preliminaries. A signatur^ of length N > 1 is a nonincreasing collection of integers 
A = (Ai > ... > X]\f) G Z'^. We will work with signatures which have only nonnegative 
parts, i.e., Atv > 0 (in which case they are also called partitions). Denote the set of all 
such objects by GTj^. Let also GT"'' := IJv=i the understanding that we identify 

A U 0 = (Ai,..., Aat, 0, 0,..., 0) G (M zeros) with A G GTj!^ for any M > 1. 

We will use two ways to depict signatures (see Fig. [^: 

(1) Any signature A G GTj^ can be identified with a Young diagram (having at most N rows) 
as in [521 I-l]- 

(2) A signature A G GTj^ can also be represented as a conhguration of N particles on Z>o 
(with the understanding that there can be more than one particle at a given location). 

We denote by |A| := A* the number of boxes in the corresponding Young diagram, and by 
i{X) the number of nonzero parts in A (which is finite for all A G GT^). For /U, A G GT"*" we will 
write /i C A if (after possibly appending p and A by zeros) we have pn < A* for all i G Z>o. In 
this case, the set difference of Young diagrams A and /x is denoted by A//U and is called a skew 
Young diagram. 

Two signatures /x, A G GT"*" are said to interlace if one can append them by zeros such that 
/X G GTj^_^ and A G GT)^ for some N, and 

Ai > /xi > A 2 > /X2 > ... > Xn-1 > Pn-1 > Aat- (2.1) 


®These objects are also sometimes called highest weights, cf. [74], as they are the highest weights of irreducible 
representations of the unitary group U{N). 
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In terms of Young diagrams, this means that A is obtained from /r by adding a horizontal strip 
(or, equivalently, that the skew diagram X/p, is a horizontal strip which is, by definition, a skew 
Young diagram having at most one box in each vertical column), and we denote this by /r A. 



I-1-•-••-1-•-1 —> 

0 1 2 3 4 5 6 


Figure 6. Young diagram A = (5,3,3,2) G GT^, and the corresponding particle 
configuration. Note that there are two particles at location 3. 


Let X' denote the transposition of the Young diagram A. For the diagram on Fig. we have 
A' = (4,4,3,1,1). If X/pis a horizontal strip, then X'/p' is called a vertical strip. We will denote 
the corresponding relation by p' -<y A'. 


2.2. Macdonald polynomials. Probability measures and Markov dynamics studied in the 
present paper are based on Macdonald polynomials. Here let us briefly recall their definition 
and properties which are essential for us. An excellent exposition and much more details may 
be found in [52l Ch. VI]. We also refer to [71 §2] for a discussion of specializations of Macdonald 
polynomials. 


Definition 2.1. Let g, t G [0,1) be two parameters. Consider the first order g-difference operator 
acting on functions in N variables: 


N 


xn) := X] n • • • ’ Xi-i,qXi,Xi+i, xn)- 


i=l jjLi 


Xi — Xi 


This operator preserves the space Q{q,t)[xi,... ,X]\[]^^^'^ of symmetric polynomials with coeffi¬ 
cients which are rational functions in g and t. 

Eigenfunctions of are given by the Macdonald symmetric polynomials P\{xi ,..., xat | q,t) 
indexed by A G GT]y, with eigenvalues 

+ ... + q^^-H + q^^)Px 


(which are pairwise distinct for generic q,t). The polynomials P\ are homogeneous, and form a 
linear basis for Q{q,t)[xi,... 

The Macdonald polynomials are stable in the sense that for any A G GT)y, 

PAuo(a;i,.. • ,®Ar,0 I g,t) = Px{xi,...,XN \ q,t)- 


Therefore, one may speak about Macdonald symmetric functions Px{xi,X 2 , ... | g, t) in infinitely 
many variables, indexed by arbitrary A G GT"*". These are elements of the algebra of symmetric 
functions, which may be viewed as a free unital algebra Sym = Q{q,t)\pi,P 2 , ■ ■ ■] generated by 
the Newton power sums pk{xi,X 2 , ■ ■ ■) = In other words, symmetric functions can 

be viewed as usual polynomials in pi,p 2 , ■ ■ ■■ Note that Px{xi, ■ ■ ■ ,xn) = 0 if i{X) > N. The 
Macdonald symmetric functions admit an equivalent alternative definition: 
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Definition 2.2. Let be the scalar product on Sym defined on products of power sums 

P\ = PXiP\2 • • • as 

{P\,P^l)q,t = zx{q,t) := ( 

^i>l 

where A = means that A has mi parts equal to 1, m 2 parts equal to 2, etc. 

The Pa’s form a unique family of homogeneous symmetric functions such that: 

(1) They are pairwise orthogonal with respect to the scalar product (•, ■)g^f 

(2) For every A, we have 

Px{xi,X 2 ,... \ q,t) = Xi^X 2 ^ ... + lower monomials in lexicographic order. 

The dependence on the parameters {q, t) is in coefficients of the lexicographically lower 
monomials 0 


iix) 

n 

2 = 1 


A, 


1-9 

1 — P 


Set bx{q, t) := 1/(Pa, Px)q,t\ this is an explicit quantity determined via the shape of the Young 
diagram A. Then the symmetric functions Qx{- \ q,t) ■= bx{q,t)Px{- \ q,t) are biorthonormal 
with the Pa’s: {Px,Qij.)q,t = 1 x=il- 

Definition 2.3. The skew Macdonald symmetric functions Qx/fj,-, GT"*", are defined as 

the only symmetric functions for which {Qx/^i-,Pu)q,t = {Qxi PiiPu)q,t for any u G GT"''. The “P” 
versions are given by Px/fi = {bij,{q,t)/bx{q,t))Qx/ii- These skew functions are identically zero 
unless fi P X. 


The skew Macdonald symmetric functions enter the following recurrence relations: 

Pa(xi, ... ,XAr) = ^ Px/i,{xi,...,XK)Pt^{xK+i,---,XN), A G GT^, K>1 (2.2) 


(and similarly for the Qa’s). This may be viewed as an alternative definition of the skew Mac¬ 
donald polynomials Pa/^ in finitely many variables. If iL = 1 in (2.2), then the summation is over 

the interlacing signatures /r A. In this case Pa/^(xi) is proportional to by homogeneity 

(cf. ( |2.4[ ) below), and (2.2) is also sometimes referred to as the branching rule for the Macdonald 
polynomials. 

From now on let us set the second Macdonald parameter t to zero. Then Pa(- | q, 0) are known 
as the q-Whittaker functions, i.e., the g-deformed gl„ Whittaker functions, cf. [38] and [T] §3]. 


Remark 2.4. Other notable degenerations of the Macdonald polynomials include the HalG 
Littlewood polynomials (for q = 0, t > 0), and the Schur polynomials (for q = t). We refer to 
[52] and |l3| for details. 


We will use g-binomial coefficients and g-Pochhammer symbols 


n\ _ {q-,q)n 

kjq' {pq)k{q;q)n-k 


{a\q)m. 


(l-a)(l-ag)...(l-ag™"^), 

< 1 , 

(1 — ag“^)“^(l — aq ~'^)~^... (1 — aq^)~^, 


m > 0; 
m = 0; 
m < 0 
(2.3) 


7 9 

Lexicographic order means that, for example, X\ is higher than const • X 1 X 2 which is in turn higher than 
const • X 2 ■ 
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to record certain explicit g-dependent quantities related to g-Whittaker functions]^ We have 


P\/tiixi I g,0) = 


\M-M 


V'A//. = V’A//.(g) := l/.AhA n ’ 


2=1 


Qx/ti{xi \q,o) = (Px/f^x^i^ 4>x/fM = 4>x/ii{q) ■= i^hA 


£(X) 

-n 

-/"i i=i 


~ '^ 2 + 1 / , 


(2.4) 

(2.5) 


Definition 2.5. A specialization of the algebra of symmetric functions Sym is an algebra mor¬ 
phism Sym —)> C. This is a generalization of the operation of taking the value of a symmetric func¬ 
tion at a point. We will deal with specializations A = (ci; (3; 7 ), where a = (oi > 02 > • • • > 0) 
P = {Pi > /32 > • • • > 0 ), 7 > 0 , and + Pi) < which may be defined via the generating 

function corresponding to signatures (n) G GT^: 




n=0 


1 -|- PiU 

'n(a,n;g)oo 


( 2 . 6 ) 


Under these specializations, we have Px/f^{A) > 0 for any /i, A G GT"*" {nonnegativity). The 
Kerov’s conjecture (see [l3l §2.9.3]) states that the specializations of the form A = (a;/ 3 ; 7 ) 
exhaust all nonnegative specializations. 

Remark 2.6. The specialization with all Pi = 0 and 7 = 0 is the same as assigning values to 
the formal variables, xj = aj. We will refer to this as the pure {a) specialization, and to the 
parameters aj as the usual parameters. 

If we go back to the case of the nonzero t parameter, then the corresponding specialization 
with all aj = 0 and 7 = 0 would send Pxf/ip I qP) to Qx'//i'{Pi^ P 2 , ■ ■ ■ \ t,q), the value of the 
usual specialization into {Pi, P 2 , ■ ■ •) with q and t swapped. (Formula ( 2 . 6 ) for nonzero t contains 
an additional factor 9)oo') Hence we will refer to such specializations as pure (/3) 

specializations, and to the parameters Pi as the dual parameters (though setting t = 0 as we do 
in the rest of the paper eliminates the “full” dual nature of these parameters). 

Finally, 7 will be called the Plancherel parameter, and the corresponding specialization can be 
defined as a limit of the specializations with, e.g.. Pi = ^ / L, i = 1,... ,L (and all other parameters 
zero), as L —)• 00 . 

The present paper mostly deals with specializations with 7 = 0 . A treatment of the case 
q; = /3 = 0 , 7>0 may be found in 


Let A U B denote the union of specializations (a generalization of concatenating the sets of 
variables). Formally it is defined aspfc(AuB) = pfc(A)-|-pfc(B), k>l. An obvious generalization 


of the recurrence relation (2.2) allows to express P\{A U B) through Px/n{A) and P^(B). Thus, 
we can equivalently say that the specialization into usual parameters is completely determined by 


(2.4) (or (2.5)) and (2.2). Similarly, the specialization into dual parameters is determined by the 
same recurrence ( 2 . 2 ), but with a different one-parameter formula: 


Qx/i,{Pi I ^, 0 ) = 


|A|-H 


Px/i, = P'x/M ■= l^^vA 


n 


(1 - q 




2>1: Xi—fii, Ai+1—/ii+1+1 


(2.7) 


^In the g-Pochhammer symbol, m may be +oo since 0 < q < 1. Note also that in all cases, (a;g)m = 

(Q.,^)oo/(n^ ,q)oo 
















(j-RANDOMIZED RSK CORRESPONDENCES AND RANDOM POLYMERS 


13 


We will also need Cauchy identities for ( 7 -Whittaker symmetric functions recorded below. Sim¬ 
ilar identities (involving t) also exist for the general Macdonald symmetric functions. 


^ Px{ai,.. .,aN)Qx{A.) = n(ai; A).. .n(aAr; A); 


AeGT+ 


P,/A(A)Q,/,(B)=n(A;B) Qx/^{B)P,/^{A). 

xeGT+ AteGT+ 


( 2 . 8 ) 

(2.9) 


In ( |2.9[ ), n(A; B) is given by 

n(A;B) =exp( 


Pn{A)pniB] 


n=l 


( 2 . 10 ) 


For the proofs see |52[ VI.(2.6) and VI.7, Example 6 ]. This definition agrees with (2.6) when 
one of the specializations is into a single usual parameter. Note also that n(A U B; C) = 

n(A;C)n(B;C). 

Finally, we will need the Fieri rules: For any r > 1, 




A: A//i is a vertical r-strip 


Q{r)PfJ. = 4’X/iiP\ 

A: A//i. is a horizontal r-strip 

( 2 . 11 ) 


(an r-strip means a strip consisting of r boxes). Here is in fact equal to the r-th 


elementary symmetric function er{xi,X 2 , ■ • •) = X) 




Xi^ ... Xi^ (note that ei = pi), and the 


Q{r)^ are the quantities entering the generating function (2.6). 

2.3. g- Whittaker processes. The (depth N) q-Whittaker processes are probability measures 
on sequences of interlacing signatures A = (A^^^ -<h • • • -<h where G GTA. Such 

sequences are sometimes referred to as Gelfand-Tsetlin schemes, or patterns, they first appeared 
in connection with representation theory of unitary groups |37) n We will depict sequences A as 
interlacing integer arrays, and also associate to them configurations of particles {{X^^\k): k = 
1,... ,N, j = I,... ,k} on N horizontal copies of Z. See Fig. Let us denote the set of all 
interlacing arrays A of depth N by 

The g-Whittaker process depends on a nonnegative specializatiorp^ A = {a; (3; 7 ) (Def¬ 


inition 


2.5) and on additional parameters a = (oi,...,aAr) with aj > 0, satisfying a^aj < 1 


for all possible i and j (this ensures the finiteness of the normalizing constant n(a; A) in (2.14) 
below). The probability weights M^(A) of interlacing arrays A may be defined via the generating 
functiorF^ 


E 


( 5 ) 


|aW| 


( 5 ) 


|A(2)|-|A(1)| 


Un \ 

On J 


|aW|-|A(^-i)| 


n(u; A) 
n(a; A)’ 


( 2 . 12 ) 


and omit the word 


®This justifies the notation “GT” we are using. 

^^In the rest of the paper, we will speak only about nonnegative specializations, 

“nonnegative”. 

^^In ( 2.12| |, n(-u; A) = Ilfni; A) ... Itfuiv; A), and similarly for the denominator (cf. (2.6l, (2.10l). Here the 
Aj’s are regarded as constants, and the Uj’s as variables. 
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plus a certain q-Gibbs property requiring that the quantities 

Mi(A) 


_ (2.13) 

(“i)^A(2)/A(i) (®2) • • • Px(N)/X(N-1) (oat) 

depend only on the top row and not on ..., Note that setting u = a turns 

( 2 . 12 ) into an identity stating that the sum of all probability weights is 1 . 


Remark 2.7. It is natural to call the property involving quantities (2.13) “g-Gibbs” because 
for g = 0 and ai = ... = oat = 1 it reduces to the following Gibbs property: The conditional 
distribution of the interlacing array A under M^(A)|g=o aj=i obtained by hxing the top row 
A(^) G is the uniform distribution on the set of all interlacing arrays A G GT^^^ with 

fixed top row A^'^^ (note that the latter set is hnite). For general g and a, the conditional 
distribution will not be uniform, but instead each interlacing array will have the conditional 
weight proportional to T’a(i) («i)^a( 2 )/a(i) ( 02 ) •.. Px{n)/ x(N-i)iaN)- 

By the Cauchy identity ( |2.8[ ) and the fact that the g-Whittaker polynomials form a linear basis, 
definition (2.12)“(2.13) is equivalent to 

^a(-^) = |^^gr^-fA(i)(®l)^A( 2 )/A(i)(« 2 ) • ■ • ^A(iV)/A(Jv-i)(aAr)(5;^(iv)(A), (2.14) 


which is a more traditional definition of the measure (first given in [7j, and earlier in |64j in 
the Schur case). To see this, one also has to note that ^ ^ 7 —v is equal to the 

product of in the left-hand side of ( 2 . 12 ) (provided that the A^-^^’s satisfy the 

interlacing constraints). 

The marginal distribution of the top row under is the q-Whittaker measure 
which is defined by either of the following equivalent ways: 

Px{u) n(u;A) 


E 

AeGT+ 


Pa (a) n(a;A)’ 

Px{d)Qx{A-) 


MM^(A) = 


n(a; A) 


(2.15) 

(2.16) 


2.4. Markov dynamics. One of the main goals of the present paper is the construction of 
Markov dynamics preserving the family of g-Whittaker processes. More precisely, we will deal 
with infinite matrices Q[B] (with rows and columns indexed by interlacing arrays) such that 

m1Q[B] = mIub, u G GtW (2.17) 

A 

(the second formula is simply an expansion of the matrix notation in the hrst formula). It suffices 
to consider three elementary cases for the specialization B which is added by the dynamics: 

(1) B = (a) is a specialization into one usual parameter a. 

(2) B = (/3) is a specialization into one dual parameter /3. (2-18) 

(3) B is a specialization with ct = /3 = 0 and 7 > 0 . 

Indeed, applying a sequence of the above elementary steps one can get a general specialization 

B (if the number of parameters a* or fij is inhnite, this also requires a relatively straightforward 

limit transition). 
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Remark 2.8. Note that setting all parameters in a specialization to zero leads to an empty 
specialization 0 . The corresponding ( 7 -Whittaker process M 0 is simply a delta measure on the 
configuration = 0 for all k,j. Note also that Q[ 0 ] is the identity matrix. 


zero 


The third case in (2.18) leads to continuous time Markov dynamics, in which the parameter 7 


plays the role of time. These continuous time dynamics were studied in detail in |16j (see also 
[59]). They are simpler than the discrete time processes (corresponding to the fist two cases in 


(2.18)) considered in the present paper. 

We will thus not focus on continuous time dynamics, and will deal with construction of matrices 
Q[a] and Q[/3] whose elements Q[a](A —)• u) and Q[/3](A —)■ iz) are transition probabilities from A to 
u (where A, 1 / G in one step of the discrete time. These matrix elements are nonnegative, 

and Q[q;](A —)• iz) = 1 for all A (and similarly for the second matrix). It is also helpful to view 
Q[a] and Q[/3] as {Markov) operators acting on functions in the spatial variables A (e.g., these 
operators act in the space of bounded functions). 

Adding a specialization B = (a) or (/3) to A as in (2.17) corresponds to multiplying the 


right-hand side of ( 2 . 12 ) by 


N 

n 


{oittj] C[)oo 

9)00 


N 


or 


1=1 


1 + Puj 




(2.19) 


respectively, since n(u; A)n(M; a) = n(u, A U (a)) and n(u; A)n(u; /3) = !!(«, A U (/3)). (Factors 
containing aj correspond to normalization, and it is the dependence on Uj in these expressions 
which is crucial.) The problem of finding Markov operators Q[a] and Q[/3] can thus be infor¬ 
mally restated as the problem of turning (by virtue of ( 2 . 12 )) the multiplication operators in the 


variables u (2.19) into operators acting in the spatial variables A. 


A similar problem of turning multiplication operators (2.19) into operators acting in the spatial 
variables A G may be posed for the generating function for the g-Whittaker measures (2.15), 


(2.16). In this case, the problem of finding the corresponding matrices IP[a] and IP[/3] (with rows 
and columns indexed by signatures A G GT)!^) has a unique solution: 

Proposition 2.9. There exist unique transition matrices 1P[q;] and lP[/3] which add specializa¬ 
tions (a) or (/3), respectively, to the q-Whittaker measure for an arbitrary nonnegative 

specialization A, in the sense similar to (2.17): 


Their matrix elements are given by 

N 


T[a](A -^u) = Yl{aaj-,q)oc^^^^^K/xa^’'^ 
J=i 


I3aj Px{d) 


( 2 . 20 ) 

( 2 . 21 ) 


where (fu/\ and explicit quantities given in (2.5) and (2.7), respectively. 


Transition operators T[a\ and T[/3] were introduced in [7], see also [6] for a similar construction 
for the Schur measures (cf. ^4. 1| below). 
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Proof. Let us consider only the case of (/3), the case of (a) is analogous. 

Multiply both sides of (2.15) by dehnition of the g-Whittaker mea¬ 

sures, the right-hand side can be rewritten as 

n(u;Au(/3)) 


E 


PAS) 


"(SAUl/S)) PAaY 

tN 


In the left-hand side, use the well-known property + Puj) = Yl^=o^r{ui^ ■ ■ ■ of 

the elementary symmetric functions [52( I.(2.2)] together with the first Fieri rule (2.11) to write 

N 

Pxiu)llil + f3u,)= ^ P,(n)</,/3H-|A|. 

j=l i/: 


(In the (a) case, one needs to use the generating fnnction (2.6) and the second Fieri rule.) Then 
the left hand side of (2.15) multiplied by O^i becomes 


N 


AeGT+ 


Px{a) 


Collecting the coefficients by Pi,{u)/Pu{a), one can rewrite this as 

E ^ E MM1(A)3>[/31(A^0. 

iAeGT+ ^ 


where the operator (P[/3] is given by (2.21). 

Since P^{u)/P^{a) are linearly independent as polynomials in u, 

^^iu(/3)H= E MM1(A)T[/3 ](A^z.) = j;MMi(A)T[/3](A^i.) 

AI A 

for all u G GT)!^. To show uniqueness suppose there is another operator T[/3]' that satisfies 

MMlT[/3]' = 


Fick Ao and uq, such that T[/3](Ao —>■ t'o) / iP[/3]^(Ao —)• t'o)- For any specialization A, 

^ MMi(A) (t[/ 3](A ^ Po) - m'iX ^ ^o)) = 0. 

AeGT+ 

Take A to be a pure specialization into usual parameters (ai,... ,0^) and multiply both sides 
by n(a; A) to get 

^ PA(a)QA(«i,...,aiv) (t[/ 3](A ^ Yo) - T[/3]'(A ^ i/o)) =0 

AeGT+ 

for any ai,... ,ai\f > 0 (in fact, this sum is only over A -<v t'o and thus is finite), which contradicts 
the fact that Qx{o:i ,..., oat) are linearly independent as polynomials in ai,..., oat. □ 
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It follows from (2.9) that both operators lP[/3] and '^[a] are stochastic, i.e. for any A G 

nm ^^)= E j’[«](a^^) = i. ( 2 . 22 ) 




Remark 2.10. If N = 1 in Proposition 


2.9 


then both dynamics IP [a] and IP[/3] (living on 


Z>o = GT(|’') are rather simple. Namely, under both dynamics, at each discrete time step the 
only particle G Z>o = GT]'" jumps to the right according to 

( 1 ) the q-geometric distribution with parameter aai, i.e., Paai{n) := (aai; ( 7 ) 00 ^^^^, n = 
0, 1 , 2 ,.. .j^in the case of dynamics IP[a], or 

( 2 ) the Bernoulli distribution with parameter /3oi in the case of dynamics lP[/3]: the parti¬ 
cle jumps to the right by one with probability /3oi/(l -|- /3ai), and stays put with the 
complementary probability 1/(1 -b /3ai)j^ 

More generally, one can show that under the dynamics on GTj^, the quantities evolve 

as follows. For 1 P[q;], at each discrete time step is increased by the sum of N independent 

g-geometric random variables with parameters aai,..., aoN- For 1P[/S], at each discrete time step 
|AW| is increased by the sum of N independent Bernoulli random variables with parameters 
/3ai,... ,l3aN- To see this, use (2.22) to write 




PxiS) 




i^SGT’ 


i+ 




[aoj ; q)oo 

N 

: JJ(1 -t- POj) 


for any A G GTj^. Substituting au instead of a (or j5u instead of (5) in these equalities leads to 


E 


N 


{aaj]q)c 


ueGTi 


^)niH-iAi = rr ^ 

[acLjU-, q)oo ' 


E ^[/3](A 


N 


//GGT^ 


^ l + POj ■ 

J=l J 


The observation follows, since both left hand sides are probability generating functions of |z^| — |A| 
in the formal variable n, and the right-hand sides expand as probability generating functions of 
sums of independent g-geometric or Bernoulli random variables. 


We will call the dynamics 7[a] and lP[/3] the univariate dynamics, and the corresponding dy¬ 
namics on interlacing arrays Q[a] and Q[/S] (which we aim to construct) the multivariate dynamics. 
In a way, multivariate dynamics on arrays A = (A^^^ • • • “<h A*-^^) stitch together univariate 

dynamics on all levels , j = 1,... ,N: Namely, started from a g-Gibbs distribution, the multi¬ 
variate evolution of the array A reduces to the corresponding univariate dynamics on each of the 


1 9 

The fact that this is indeed a probability distribution follows from the g-binomial theorem. 

1 

This parametrization of Bernonlli random variables will be used throughout the paper. 
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levels j = 1,..., A^. This fact follows from the proof of Theorem 2.13 below, see also m 
§2.2] for a related discussion. 


Instead of the case of univariate dynamics (driven by identity (2.15)), the problem of con¬ 
structing multivariate dynamics (involving identity ( 2.12| )) has a whole family of solutions. This 
phenomenon was known in the Schur (q = 0) case for some time, with the presence of the RSK- 
type (e.g., see |53], [55]) and the push-block [T3| dynamics (see ^ below for more detail). A 
similar phenomenon was investigated in |16j for continuous time dynamics increasing the param¬ 
eter 7 in the q-Whittaker processes. In that simpler continuous time setting, a classihcation result 
was established in the latter paper. 


Remark 2.11. Since the g-Whittaker polynomials Px{a) entering (2.20) and (2.21) are not given 
by an especially nice formula, transition probabilities of the univariate dynamics are harder to 
analyze. On the other hand, RSK-type multivariate dynamics which we construct in the present 
paper turn out to have simpler transition probabilities. Note also that multivariate dynamics 
on g-Gibbs distributions can be used to “simulate” the univariate ones, cf. the above discussion 
about “stitching”. 


2.5. Main equations. Here we write down linear equations whose solutions correspond to mul¬ 
tivariate discrete time Markov dynamics on g-Whittaker processes. Let us hrst narrow down the 
class of dynamics Q on interlacing arrays which we consider. 


Definition 2.12. A dynamics Q on interlacing arrays will be called a sequential update dynamics 
if its one-step transition probabilities from A to iz, A, iz G have a product form 

Q(A ^ iz) = 1 [i(A(^) ^ p(^))U2(A(2) ^ 1/(2) I ^(1) ^ ^(1)^ ^ ^ ^ ^ I A(^-^) ^ z/(^-^)), 

(2.23) 

where Uj’s are conditional probabilities of transitions at levels J = 1,..., satisfyinj^ 

Uj(A(^) ^ z/(^) I A(^-i) ^ z/(^-^)) >0, ^ IIj(A(^( ^ z/(^( I A(^-^) ^ z/(^-^() = 1. (2.24) 

i/0')eGT+ 

In words, the transition A —)• iz looks as follows. First, update A(^'* —)• z^(^( at the bottom level GT]*" 
according to the distribution IIi. Then for each j = 2,..., N, given the transition A(-^“^( —)• 
at the previous level, update —)• at level GTA according to the conditional distribution Uj. 
We see that the evolution of several first levels A(^\ ..., A(*^( of the interlacing array does not 
depend on what is happening at the upper levels A(^^^\ ..., A('^(. 


This setting of sequential update dynamics is not too restrictive as it covers all previously known 
examples of dynamics on Macdonald (in particular, on g-Whittaker and Schur) processes, cf. |16j . 
For a sequential update dynamics it suffices to describe the evolution at any two consecutive 
levels j — 1 and j. 


Theorem 2.13. A sequential update dynamics Q defined via (2.23) 
q-Whittaker processes 

if 


(2.24) preserves the class of 
and adds a new usual parameter a to the specialization A if and only 


Uj{X ^ 1/ I A ^ z/)(aaj)l^l '""'Va/a'/’p/a = (o:aj;q)oo'ip^/p(f>iy/x (2-25) 

AeGT+_j 


agreement, for j = 1 we mean —>■ \ 
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for any j = 1,2,... ,N and any \,v ^ v G such that the four signatures X, v, X, v 

are related as on Fig. left (in particular, the above summation is taken only over X satisfying 
X u, X X). For j = 1 we take X = v = 0 in this equation and it becomes equivalent to 
Ui = !P[a] at level GT^,^ (as in Remark 2.10). 

Similarly, a dynamics Q preserves the class of q-Whittaker processes and adds a new dual 
parameter (3 to the specialization A if and only if 

Uj(X ^ 1/ I A ^ (2-26) 


AeGT' 


i-1 


for any j = 1,2,... ,N and any A, G GT^, v G such that the four signatures X, v, X, v 

are related as on Fzg. right (in particular, the above summation is taken only over X satisfying 
X v, X ^)- For j = 1 we take A = = 0 in this equation and it becomes equivalent to 

Ui = IP[/3] at level GT^ (as in Remark 2.10). 


(a) 


3 - 1 


A jy 

f Y 

A 


(/3) 


A jy 

f Y 


i -1 


time 


time 


Figure 7. Squares of four signatures on two consecutive levels relevant to conditional 
transition A —?■ z^ on the upper level given the transition A —)• u on the lower level, under 
dynamics Q[q;] (left) and Q[/3] (right). Note the similarity to blocks in Fomin’s growth 
diagrams (about the latter, see [30], |3T|, |33|, [33]1. 


The proof of these equations was already established in |16[ §2.2] using a more general frame¬ 
work of Gibbs-like measures. However, for the sake of completeness, we reproduce it here in our 
particular setting of the ( 7 -Whittaker processes. 


Proof. Let us consider only the case of adding (a), as the case of (/3) is analogous. 

The fact that a sequential update dynamics Q dehned via (2.23)-(2.24) preserves the class of 
( 7 -Whittaker processes and adds a new usual parameter a to the specialization A means that 


^ n(a; A) F\W {oi)Px(2) /A(i) («2) • ■ • P\(N) /A(JV-i) {aN)Q\(N) (A) 


Ui(A(i) ^ z.(i))lt2(A(2) ^ 1.(2) I A(i) ^ z.(i)).. .l[Ar(A(^) ^ z.(^) I A(^-i) ^ z.(^-i)) = 

n(a- A^U («i)^i.( 2 ) /uW (« 2 ) • • • T),(iv) {aN)Q„(N) (A U (a)) for every u and A. 

(2.27) 


Using (2.2), we can rewrite (2.27) as 


N 




yd) I \{3 1) jyw -/ 




V’AO)/A(i-i) Qaw("^) = 
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Y{{aaj;q)ooip^U)/^U-i) \ ^ ^''^4>u(n)/u for every i/ and A. 


vi=i 




Since Q\{A) are linearly independent as polynomials in ui,...,U]\f for a specialization A into 
usual variables (ui,... ,un), this is equivalent to saying that 


N 

Y, n U,(aO> ^ I A«-‘> ^ = 

A:A(^)=Ai=l 

N 

=‘/’iaW/aW ]J(aaj;g)ooV'i.0)/,,0-i) for all IZ and A. (2.28) 
i=i 

For the proof in one direction, suppose that Ui = lP[a] at level GT]*", and Uj{X^^^ —>■ | 

AO-1) —7- satisfy (|2.25[) for 2 < j < N. Then we can show by induction on k, that 


^ nit,(AW ^ I ^ = 

A:A('“)=Ai=l 

k 

~ 4’u'^k) j-yik) n {aaj-,q)oo'ip^U)/uU-->-) fo^ all 1 < A: < A, i/ = -<h A G GT+. 

i=i 

(2.29) 

Base for /c = 1 follows from the fact that Ui = lP[a] at level GT]'', while the inductive step follows 
from (2.25). So (2.28) holds. 


For the other direction, suppose that (2.28) (and hence (2.27)) holds. For 1 < A: < A by 
summing (2.27) over ..., and applying (2.9) we get 


YjT-- „ . aa -Pa(i)(Q'i)-Pa(2)/a(i)(Q2) • . . Px(k)/x(k-i)iak)Q\(k){A) 

^ mai, IVj 

lti(A(^) ^ z/(i))U 2(A(2) ^ 1/(2) I ^(1) ^ ^(1)^ ^. .ltfc(A(^( ^ p(*^( I A(^-^) ^ z/(^-^)) = 

= - , .. Puw{ai)P^( 2 )i^(i){a 2 ).. .Pj,(fc)/^(fc-i)(afc)(5^(fe)(AU (a)) for every i/and A, 

ll(ai,...,afc; AU (a)) ' ' 


which implies (2.29). For A: = 1 it means that Ui = lP[a] at level GT(|^, while for k > 2 using 
(2.29) for both k and k — 1 implies (2.25). □ 


In a continuous time setting, there also exist linear equations governing multivariate dynamics. 


cf. [161 §2.4]. In fact, the latter equations arise as small a or small /3 limits of (2.25) or (2.26), 
respectively. Markov dynamics on g-Whittaker processes corresponding to solutions to these 
continuous time equations were constructed in [59] , [l6| , [l8| . 


2.6. Discussion of main equations. Let us make a number of general remarks about the main 


equations of Theorem 2.13 
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2.6.1. The paper |16j contains a classification result in continuous time setting, which was 
achieved by further restricting the class of dynamics by imposing certain nearest neighbor in¬ 
teraction constraints. Under these constraints, putting together continuous time linear equations 


(which look similarly to (2.25) and (2.26)) with fixed A and P in a generic position, at level j one 
arrives at a system of j linear equations with 3j — 2 variables. Solutions of such a system admit 
a reasonable classihcation. 

It remains unclear how to impose (preferably, natural) constraints on solutions of discrete time 
equations (2.25) or (2.26) so that the family of all solutions would admit a reasonable description. 
Indeed, for example, in the case of a usual parameter ( |2.25 ), the number of variables is infinite 
while the number of available equations is finite. Therefore, in ^ and ^ below we devote 
our attention to constructing certain particular multivariate discrete time dynamics satisfying 


equations (2.26) and (2.25), respectively. 


2.6.2. Note that summing (2.25) or (2.26) over i' G GT^ leads to the skew Cauchy identity with 


both specializations being into one parameter (cf. (2.9)): 

I 


E 

AeGT+ 




n(a,;B) ^ 

^ ^ iyeGT+ 


Pu/uiaj)Q^/\(B), B = (a) or (/3). 


(2.30) 


Identity (2.30) may also be interpreted as a certain commutation relation between the uni¬ 


variate Markov operators IP[q:] or (P[/3] (of Proposition 2.9) and Markov projection operators (or 
linksf^ 




A E GT+, A E GT+_1, 


in the sense that 




(2.31) 


and similarly for lP[/3]. Indices j and j — I in T[a] above mean the level of the interlacing array 
at which the transition operator of the univariate dynamics acts. 


One can thus say that each solution to the main equations (2.25) or (2.26) (and, therefore, 
each discrete time Markov dynamics on g-Whittaker processes) corresponds to a refinement of 
the skew Cauchy identity (2.30) (or of the commutation relation (2.31)). 


Remark 2.14. When B = (a) is a usual specialization, one may check that all quantities entering 


both sides of (2.30) can be viewed as generating series in q, a, and a,- with nonnegative integer 


coefficients. It would be very interesting to understand whether there is a bijective mechanism 


behind identity (2.30) similar to the one existing in the classical {q = 0) case (see also the 
discussion after Theorem 4.7). We are very grateful to Sergey Fomin for this comment. 


2.6.3. The parameters ui,..., aj-i (but not afi essentially do not contribute to the main equa¬ 
tions (2.25), (2.26): they enter the equations only as a requirement that u E GT^_^ and 


A,j/ E GT+. Thus, equations (2.25), (2.26) essentially depend on two specializations: a spe¬ 


cialization into one usual parameter A = (aj) which corresponds to increasing the level number, 
and a specialization B = (a) or (fi) which corresponds to time evolution. This allows to think of 
diagrams as on Fig.j^ as well as of main equations, for any specializations A and B (see Fig. [^. It 
suffices to consider three elementary cases for A and B as in (2.18). This yields 9 possible systems 


15. 


These links in fact determine the g-Gibbs property (2.131; e.g., see [161 §2] for more detail. 
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time 


Figure 8 . A square of four signatures corresponding to arbitrary specializations A 
and B. Notation A A means that > 0, and similarly for ^b- When the 

specialization A is into a single usual or dual parameter, ^a reduces to or -<v, 
respectively. 


of equations for dynamics. If one of the specializations is pure Plancherel (case (3) in ( 2.18| )), 
then the corresponding Markov dynamics on g-Whittaker processes were essentially constructed 
in m,m- This leaves four systems of equations in which both A and B are specializations into 
a single usual or dual parameter. In this paper we address two of these four cases corresponding 
to A = (aj), which in particular give rise to two new discrete time g-PushTASEPs (as marginally 
Markovian projections of dynamics on interlacing arrays, see 35.2| and 


2.6.4. In fact, one can define the quantities <i ^\/for the general Mac¬ 

donald parameters {q,t) (see [521 Ch. VI]), and thus write down the corresponding main linear 
equations for any specializations A and B. (In particular, for t 7 ^ 0 the right-hand side of the 
identity (2.6) defining a specialization should be multiplied by '?)oo-) It is not known 

whether there exist other solutions to the main equations for general {q, t) yielding honest Markov 
dynamics (i.e., having nonnegative transition probabilities) except the push-block solution (see ^ 
below for the definition). We do not address this question in the present paper. 

There is a rather simple transformation of the main equations for general {q, t) (related to 
transposition of Young diagrams) which interchanges q ^ t and swaps usual and dual parameters 
in both specializations A and B [18]. This transformation relates the g-Whittaker {t = 0) and 
the Hall-Littlewood {q = 0) settings. 

The remaining two cases of the (( 7 -Whittaker) main equations mentioned above (corresponding 
to A = ( 6 ), a specialization into a dual parameter) should thus be thought of as discrete time 
versions of the continuous time equations of [18] (relevant to the HalELittlewood setting). As 
such, (conjectural) solutions to the former equations leading to discrete time dynamics on inter¬ 
lacing arrays are unlikely to produce new marginally Markovian TASEP-like particle systems in 
one space dimension (see also discussion in [Tb] §8.3]). In the present paper, we do not address 
these two remaining cases corresponding to the Hall-Littlewood setting. 


3. Push-block and RSK-type dynamics 

3.1. Push-block dynamics. There is a rather straightforward general construction (dating back 
to an idea of m) leading to certain particular multivariate dynamics. Namely, ass ume t hat the 
conditional probabilities IXjAA —)• u j A —?■ u) entering the main equations (Theorem 


2.13) do not 


depend on A. Then each equation (corresponding to fixed A,u G GT^, and u G GT^_^) contains 
only one unknown Vij{X —)• u j u). With this restriction the main equations admit a unique 
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solution. Let us consider the case of a usual parameter a (2.25). Observe that the left-hand side 
of (2.25) takes the following form (where signatures satisfy conditions on Fig. left): 


llj(A 




E 

AeGT+_^ 


I aa,' 


= Uj(\ 

= Uj(A 


ly ly' a 


|A|-|i.i -hl+|p| 


^A/A(«i)<3p/A(«) 


aggt; 


j-i 


ly I ly) (aaj;q)oc ^ (aoi)'""' '""'V’x/p'/’x/a, 


xeGTJ 


where we have used the skew Cauchy identity (|2.30|). Then (|2.25|) yields the solution 

1I,(A 


u \ u) = 




SxeGT+ ('^®i)^^^V'x/i/0x/A 


(3.1) 


In ( |3.1[ ) as well as in the above computation, it should beA-<hi^, F-<hZ^ and X x, u x, 
see Fig. left. 

Similarly, the solution of (2.26) not depending on A looks as 


II j (A —y ly \ — 


(3.2) 


The signatures have to be related as on Fig. right, i.e., \ v and \ x^v x. 


Definition 3.1. We will call the multivariate dynamics dehned by (3.1) or (3.2) the {discrete 
time) push-block dynamics on g- Whittaker processes adding a specialization (a) or (/3), respec¬ 
tively. About the name see [4.2 below. We denote these dynamics by Qp)((°[a] and 

The construction of push-block dynamics can be equivalently described as follows. Recall the 
commutation relation between the univariate dynamics T and the stochastic links Aj_j^ (2.31). 
Then one can say that the multivariate dynamics chooses zz at random according to the distribu¬ 
tion of the middle signature in a chain of Markov operators 

■pU) 

A- ly — - -)• u, 


conditioned on the hrst signature A and the last signature u. Denominators in formulas (3.1) and 
(3.2) reflect this conditioning. 


The push-block dynamics (in the Schur case) hrst appeared in m, see also ^4.2| below. For 
analogous dynamics in continuous time and space in which the univariate dynamics is the Dyson’s 
Brownian motion see m- As shown in m, the latter dynamics is a diffusion limit of the discrete- 
space dynamics from m- 


3.2. RSK-type dynamics. Let us now dehne an important subclass of multivariate dynamics 
which is central to the present paper. 


Definition 3.2. A multivariate sequential update dynamics Q (which corresponds to conditional 
probabilities Uj{X ^ u \ X ^ 1 ?) satisfying 
called RSK-type if 

'U.j{X ^ u \ X ^ u) = 0 unless |z^| — |A| > |p| — |A|, for all X,iy ^ GT^, X,iy ^ 


(2.24) and the main equations (2.25) or (2.26)) is 
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In the above definition, \9\ — |A| is the total distance traveled by particles at level j — 1, and 
similarly |i/| — |A| is the total distance traveled by particles at level j. Informally, under an RSK- 
type dynamics all movement at level j — 1 must propagate further to level j (and, consequently, 
to all upper le vels o f the array). 

under an RSK-type dynamics the quantity |A^-^)| — (for any j = 


By Remark 


2.10 


1,..., N) at each step of the discrete time is increased by adding a g-geometric random variable 
with parameter aaj (in the case of Q[a]), or a Bernoulli random variable with parameter 13aj (in 
the case of Q[/3]). 

Remark 3.3. This feature of RSK-type dynamics separates them from the push-block dynamics 
of §3.1[ Indeed, under a push-block dynamics movements at level j —1 generically do not propagate 
upwards because the quantities Uj (A —)• | A —)• p) do not depend on A. More precisely, the only 
steps at level j — 1 that can propagate to level j correspond to the situation n X. Then a part 
of the movement A —?• is mandatory, as it is dictated by the need to immediately (i.e., during 
the same time step of the multivariate dynamics) restore the interlacing between the levels j — 1 
and j. 

RSK-type dynamics on g-Whittaker processes that we construct in ^ and ^ give rise to 
discrete time g-TASEPs and g-PushTASEPs as their Markovian marginals. On the other hand, 
discrete time push-block dynamics do not seem to produce any TASEP-like processes Note 
also that in general the denominator in (3.1) or (3.2) does not seem to be given by an explicit 


formula, so the discrete time push-block dynamics are not easy to work with (cf. Remark 2.11). 
This provides an additional motivation for constructing and studying RSK-type dynamics. 

4. SCHUR CASE 

In this section we discuss the Schur (g = 0) case, and explain how in this case the RSK-type 
multivariate dynamics are related to the classical Robinson-Schensted-Knuth correspondences. 

4.1. Univariate dynamics in the Schnr case. When g = 0, the g-Whittaker polynomials P\ 

and Qx reduce to the simpler Schur polynomials sx- In particular, we have ipx/^i = 

and i’x/u ~ l/iAvA- Univariate discrete time dynamics on the first level = Z>o look as in 


Remark 2.10 with the understanding that the g-geometric distribution in the case of (P[a] has to 
be replaced by the usual geometric distribution Paai(n)|q=o = (1 — aai)(aoi)”, n = 0 , 1 , 2 ,.... 

Remark 4.1. The continuous time dynamics on GT(|^ increasing the parameter 7 of the spe¬ 
cialization is the usual Poisson process which can be obtained from either of the discrete time 
dynamics lP[a] or lP[/5] in a small a or small (3 limit, respectively. In fact, this observation is also 
true in the general g > 0 case. 

The univariate dynamics tP[a] and T[/S] at any higher level GT)^, N = 2,3,... (described in 
a g = 0 version of Proposition 2.9), can be obtained from the N = 1 dynamics via the Doob’s 
h-transform procedure. Informally, to get the dynamics of N distinct particles (xi > ... > xn) 
on Z>o (this state space is the same as GT~^ up to a shift Xi = Xi + N — i), one should consider 
the dynamics of N independent particles xj each of which evolves according to the corresponding 


16 . 


The continuous time push-block dynamics on g-Whittaker processes has lead to the discovery of the continuous 
time g-TASEP in [7]. A continuous time RSK-type dynamics on 5-Whittaker processes was later employed in m 
to discover the continuous time g-PushTASEP, a close relative of the g-TASEP (see also (5.6 below). In fact, q- 


Push TAS EP and g-TASEP can be unified to produce another nice particle system on Z, namely, the q-PushASEP 
(see (7.5 below), which also extends to a certain dynamics on interlacing arrays m- 
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N = 1 dynamics, and then impose the condition that the particles never collide and have relative 
asymptotic speeds ai,..., ajsf, respectively. This conditioning gives rise to the presence of the 
factors Su{a)/sx{a) in transition probabilities (cf. Proposition 2.9). We refer to, e.g., [JH], [17] . 
|55| . |62j for details on noncolliding dynamics. 


It is worth noting that the Dyson’s Brownian motion coming from NxN GUE random matrices 
m arises via a similar procedure by considering noncolliding Brownian particles. One may thus 
think that the univariate dynamics 7[a\ and IP[/3] on GT^ are certain discrete analogues of the 
Dyson’s Brownian motion. 


4.2. Push-block dynamics in the Schur case. Setting q = Q greatly simplifies formulas (3.1) 


and (3.2) thus leading to nice push-block multivariate dynamics on interlacing arrays. They were 


introduced and studied in 

Due to the sequential nature of multivariate dynamics (2.23), we will consider evolution at 
consecutive levels j — 1 and j. Assuming that the movement A —>■ i/ at level j — 1 and the old 
configuration A at level j are given, we will describe the probability distribution of G GT^ 
corresponding to ^[^(A —)• i/ | A —)• i?). 


O + T 3 


1 -tO 


2 + F2 


2-tl 


7 + yi 


block 


+1 


A -I- - A) 


A + - A) 


2-tl 

?=0ri 


4-tl 


Figure 9. An example of a step of Qpj^ [/3] at levels 4 and 5. Here A = (7, 2, 2,1,0), 
A = (4, 2,1,1), and D = (5, 3,1,1). The move A 2 = 2 —)• 1/2 = 2 -|- 1 on the upper level 
is dictated by the corresponding move A 2 = 2 —>■ U 2 = 2 -|- 1 on the lower level (due to 
the short-range pushing mechanism), so no further move of 1 x 2 is possible. The particle 
A 4 = 1 cannot move because it is blocked by U 3 = A 4 . All other particles are free to 
move (including A 3 which was blocked before the movement at the lower level), and their 
jumps Ti, I 2 ; I 3 are independent identically distributed Bernoulli random variables with 
P(yi = 0) = l/(l + /3a,). 


Let us first focus on the case of Qpb°[ 


^^(A —>■ u I u) = 


(see Fig. 9| ' In this case (3.2) simplifies to 


(/3a,-)l"llp 




SxeGT+ 

i.e. for any u', v" G GT^, such that v u', A -<v u' and v u", A -<v u", 

UjjA^n' I u) _ xW\-y'\ 

It is clear that the only dynamics with such property fits the following description. During one 
step of the dynamics, each particle Aj, 1 < i < j, can either stay, or jump to the right by one, 
according to the rules: 


simplify pictures, here and below we will display interlacing arrays of integers (cf. Fig. [^, but will still 
speak about particles jumping to the right. 
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(1) {short-range pushing) If Ui = Xi + 1, then the move Aj —)■ Pj = Aj + 1 is mandatory to 
restore the interlacing (which was broken by the move A —)• z^) during the same step of 
the discrete time. 

(2) (blocking) If A* = Ui-i, then the particle A* is blocked and must stay, i.e., r'j is forced to 
be equal to A*. 

(3) (independent jumps) All other particles A* which are neither pushed nor blocked, jump to 
the right by 0 or 1 according to an independent Bernoulli random variable with probability 
of staying 1/(1 + (daj). 


o + n 


1 + 0 


2 + I3 


2 + 2 +y2 7+1 +n 


block 



2 + 2 


4 + 4 


A + (z^ - A) 

A + (z^ - A) 


Figure 10. An example of a step of at levels 4 and 5. The move Ai = 4 —)■ = 

4 + 4 forces Ai to move to the right by 1, and similarly the move A2 = 2—>-U2 = 2 + 2 
forces A2 to move to the right by 2 (short-range pushing); note that these forced moves do 
not exhaust all possible distance traveled by Ai or A2. The particle A4 = 1 is blocked by 
U3 = A4 and thus cannot move. All other parts of the movement A —)■ u are determined 
by independent identically distributed geometric random variables 4/, 1 < i < 4 with 
parameter aaj, where each variable is conditioned to stay in the maximal interval not 
breaking the interlacing: Ti > 0 (i.e., no conditioning), 0 < ^2 < 4:, 0 < T3 < 2, 

0 < 44 < 1. 


see 


By the same explanation the dynamics Qp^^[Q:] at two consecutive levels looks as follows 
Fig. 10). Each particle Aj, 1 < i < j, independently jumps to the right by a random distance 


which has the geometric distribution with parameter aaj conditioned to stay in the interval from 
(i>i — Aj)+ := max{0, u, — A*} to Uj_i — Aj (with the agreement that Aq = +cx))|^ This conditioning 
corresponds to the denominator in (3.1). 


4.3. RSK-type dynamics in the Schur case. Let us now discuss four discrete time multi¬ 


variate RSK-type dynamics Q?, 


=0r 
row 


Q?. 


=0r 






on Schur processes. The former 


two dynamics arise from the row RSK algorithn^^ applied to geometric or Bernoulli random 
input, respectively (cf. Remark 4.5 below). Similarly, the latter two dynamics correspond to 

[71] for 


the column RSK algorithm applied to the same random inputs. We refer to 
relevant background and details on RSK correspondences (though descriptions of dynamics below 
in this subsection may serve as equivalent definitions of RSK algorithms). See also |16( §7] for a 
“dictionary” between interlacing arrays and semistandard Young tableaux viewpoints. 

Let us first recall two elementary operations of deterministic long-range pulling and pushing 
from |16j (in the language of semistandard Young tableaux they correspond to row and column 
bumping, respectively). 


18 


Due to the memorylessness of the geometric distribution, this description is equivalent to what is illustrated 


on Fig. 


19. 


10 


The row RSK is the most classical version of the Robinson-Schensted-Knuth algorithm. 
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Definition 4.2. (Deterministic long-range pulling, Fig. 11) Let j = 2,...,N, and signatures 
A, P G GT^_^, A G GT+ satisfy A -<h A, z? = A -|- e*, where e* = (0, 0,..., 0,1,0,..., 0) (for some 
i = 1,..., j — 1) is the ith basis vector of length j — 1. Define u G GT^ to be 


ly = pull(A I A —)• z/) := 


A -|- ej, 

A -|- ej+i, 


if A^ — A^, 
otherwise. 


Here e^ and e^+i are basis vectors of length j. 

In words, the particle A* at level j — 1 which moved to the right by one generically pulls its 
upper left neighbor Aj+i, or pushes it upper right neighbor A^ if the latter operation is needed 
to preserve the interlacing. Note that the long-range pulling mechanism does not encounter any 
blocking issues. 



1 + 1 3 


7 

V • 

+1 

\ 


• 

• 

• 


2 + 1 

4 



Figure 11. An example of pulling mechanism for z = 2 at levels 2 and 3 (i.e., j = 3). 
Left; A 2 = A 2 , which forces the pushing of the upper right neighbor. Right: in the 
generic situation A 2 < A 2 the upper left neighbor is pulled. 


Definition 4.3. (Deterministic long-range pushing. Fig. 12) As in the previous definition, let 
J = 2,..., A", A, u G GT^_;^, A G GT^ be such that A A and u = A -|- e*. Define u G GT^ to be 

ly = push(A I A —)■ u) := A -|- e^, where m = max{p: 1 < p < z and Xp < Ap_i}. 


In words, the particle A* at level j — 1 which moved to the right by one, pushes its first upper 
right neighbor Am which is not blocked (and therefore is free to move without violating the 
interlacing). Generically (when all particles are sufficiently far apart) Am = Aj, so the immediate 
upper right neighbor is pushed. 


Remark 4.4 (Move donation). It is useful to equivalently interpret the mechanism of Definition 


gives the particle Aj at level j a moving impulse. If A* is blocked (i.e., if A* = Aj-i), this moving 
impulse is donated to the next particle Aj_i to the right of Aj. If Aj_i is blocked, too, then the 
impulse is donated further, and so on. Note that the particle Ai cannot be blocked, so this moving 
impulse will always result in an actual move. 


4.3 in a slightly different way. Namely, let us say that when the particle A* at level j — 1 moves, it 


Let us now describe four RSK-type dynamics on Schur processes. Under each of the dynamics 
the interlacing array A is updated sequentially (cf. ( |2.23[ )) at each level j = 1,N. At each 
step of the discrete time corresponding to an update A —>■ new randomness is introduced via N 
independent random variables Vi,..., Vn, which are either geometric random variables (belonging 
to Z>o) with parameters aai ,..., aa^ in the case of Q^dw [a] and Bernoulli random 

variables G {0,1} with parameters /3ai,..., (doN in the case of and These random 

variables are resampled during each time step. 




















g-RANDOMIZED RSK CORRESPONDENCES AND RANDOM POLYMERS 


28 



A + — A) 
A + - A) 


Figure 12. An example of pushing mechanism for i = 3 at levels 4 and 5 (i.e., j = 5). 
Since the particles A 3 = A 2 and A 2 = Ai are blocked, the first particle which can be 
pushed is Ai. 


Remark 4.5. We see that all randomness in each of the four RSK-type dynamics can be organized 


into a matrix {yj^^)i<j<N,t=i, 2 ,... (with appropriate distribution of the 


W) 


s). Such matrices 

containing nonnegative integers are usually thought of as inputs for classical Robinson-Schensted- 
Knuth correspondences. 


Under each of the four dynamics, the particle at the first level of the array is updated as 
= A^^^ + Ui. Then, for each j = 2,..., N, assume that we are given signatures A, u G 


A G GT^ satisfying relations as on Fig. 
(/3) of the dynamics). Let us represent t 


7 (note that these relations depend on the type {a) or 




le movement A —?■ u at level j — 1 as 


u-X = Y^ aei 

2 = 1 


in the case of Qrow[al and QLuial 


Ci G Z>o 
Ci G {0,1} in the case of Qn 


row I 
ro=w[ 


g=0r 

•"col 

=0r 


and 


(recall that e^ is the ith basis vector of length j — 1 ). Also denote |c| := 

Depending on the dynamics, we will construct the signature v G GT^ (which also fits into 
relations on Fig. as follows: 


• (Qrow[a], Fig. |l3[ ) First, do |c| operations pull (Definition 4.2) in order from left to right, 
starting from position j — 1 all the way up to position 1. In more detail, let — 1,0) := A and 
for p = 1 ,..., Cj-i let 

hU - 1,P) := Pull(M(j - - 1) I A + (p- l)ej_i X + pej-i), 

then let fj,{j — 2 , 0 ) := /i(j — 1 , Cj_i) and for p = 1 ,..., Cj -2 let 

/i(j - 2 ,p) := pu\\{fi{j - 2,p - 1) I A + Cj-iej-i + (p - l)ej -2 A + Cj-iej-i +pej-2), 

etc., all the way up to //(l,ci) := pull(p(l,ci — 1) | u — ei —)• u). (Clearly, if some Cj = 0, then 
the steps corresponding to p,{i, •) should be omitted.) 

After these |c| operations, define u := p(l,ci) + Vj-ei. That is, let the rightmost particle at 
level j jump to the right by Vj (which is a geometric random variable with parameter auj). 

• (Qrcw[/3], Fig. 14) First, dehne /r(l,0) := A + I^ei. That is, let the rightmost particle at level 
j jump to the right by Vj (which is a Bernoulli random variable with parameter /3aj). 

After that, perform |c| operations pull (Definition |4.2[ ) in order from right to left, starting from 
position 1 all the way up to position j — I (details are analogous to the above dynamics Q?c 1 ^[q;]). 
Then set u := p,{j — 

• (QS [a]. Fig. 15^ First, the leftmost particle Xj at level j receives Vj moving impulses (here 
Vj is a geometric random variable with parameter auj). Each moving impulse means that Xj 
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1 + 1 


4+1 



3 + 2 


6 


6 6 + 3 9 + 1 + 



6 + 1 7 + 3 


A + (z^ — A) 
A + (z^ — A) 


Figure 13. An example of a step of at levels 4 and 5. Propagation steps 

(represented by numbers on arrows) are performed from left to right, according to pull 
operation. After that, the rightmost particle at level j jumps to the right by Vj. 



Figure 14. An example of a step of Q?ow [/3] at levels 4 and 5. Propagation steps are 
performed from right to left, according to pull operation. Above: Vj = 1, below: Vj = 0. 


tries to jump to the right by one, and if it is blocked (i.e. 
is donated to Aj_i, etc. (see Remark 4.4) 
moving impulses by — 1 , 0 ). 

After that, perform |c| operations push (Definition |4.3| 


if \j = Aj_i), then the moving impulse 
Denote the signature at level j arising after these Vj 


in order from left to right, starting 
from position j — 1 all the way up to position 1 (details are analogous to the above). Then we 
set u := /i(l, Cl). 


1 + 2 4 + (Uj-2) + l 6 7+1 + 1 11 + 2 + 1 

A) 


A) 

3 + 2 6 6 + 3 9+1 

Figure 15. An example of a step of levels 4 and 5. We have Vj = 3, which 

means that initially the particle A 5 jumps to the right by 2 and the particle A 4 jumps 
by 1 (because of move donation). After that, propagation steps are performed from left 
to right, according to push operation. 
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(Q^oi*^[/ 5]) Fig. 16) First, perform |c| operations push (Definition |4.3|), in order from right to 


left, starting from position 1 all the way up to position j — 1 (details are analogous to what is 
done above). Let fj,{j — be the signature at level j arising after these |c| operations. 

After that, let the leftmost particle at level j receives Vj moving impulses (here Vj is a Bernoulli 
random variable with parameter /3aj). That is, if Vj = 0, then set v := fi{j — 1, cj-i). Otherwise, 
if Vj = 1, the jth particle at level j tries to jump to the right by one. If it is blocked, the impulse 
is donated to the (j — l)th particle at level j, etc. In this case, denote by u the signature at level 
j arising after this moving impulse. 


3 6 + Vj 6 6 + 1 8 + 1 



3 6+1 6 7+1 


A + (r" — A) 


A + (i/ - A) 


Figure 16. An example of a step of at levels 4 and 5. Propagation steps are 

performed from right to left, according to push operation. We have Vj = 1, and the 
jump of the rightmost particle at level j is donated to the right. 


The above four rules of constructing the signature u G GT^ complete the description of the 
RSK-type dynamics and respectively. 

Remark 4.6. By the very construction, at each step of any of the four above RSK-type dynamics 
the quantity is increased by Vi + ... + Vn, as it should be (cf. the discussion before 

Remark 3.3). 


9=0 r 


=0r 


], and 


described above 


Theorem 4.7. The RSK-type dynamics Q?ciw [a], Q?, 

satisfy q = 0 versions of the main equations of Theorem 2.13\ and hence act on Schur processes 
by adding a new usual parameter a or a new dual parameter (3, respectively (as in ( 2.17 | ) ). 

Proof. This statement follows from bijective properties of RSK correspondences (briefly discussed 
below in this section), or, equivalently, it may be regarded as q = 0 degeneration of our main 
results about RSK-type dynamics on g-Whittaker processes (Theorems 5.2, 5.7, 6.4, and 6.11). □ 

Each of four RSK-type dynamics described above gives rise to a certain bijection between sets 
{X,X,r',Vj} and {A,u, u} (at each time step and at each level j of the interlacing array). In 
more detail, each of the dynamics Qrciw [a] and (see Fig. 13 and 

between the following sets: 


15) produces a bijection 


{A, A, u; A +h A -<h U {V,- G Z>o} ^ {A, u, u: A u +h u}. 

(Fig. 


(4.1) 


Similarly, each of the dynamics Q?, 
between the sets 


ro=w[ 


and 


14 


and 


16) establishes a bijection 


{a. A, u: A +h A -<v u} U {kj- G {0,1}} <—)■ {A, u, u: A -<v u +h u}. 


In (4.1) and (|4.2|) we have A, u G GT7 ^ and A, u G GT7, as usual. 


(4.2) 
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in E 


The understanding of RSK correspondences via bijections as in (4.1) and (4.2) was presented 
It also implies that fixing A and v and taking generating functions of both sets in (4.1), 


by weighting elements of the left set by and of the right set by (under 

the bijections, these powers are equal to each other), one recovers the skew Cauchy identity (2.30) 
for B = (a). Similarly, (4.2) leads to (2.30) with B = (/3). This observation agrees with the 
understanding of multivariate dynamics as refinements of the skew Cauchy identity (^2.6.2 ). 


Remark 4.8. In RSK-type dynamics on g-Whittaker processes considered in ^and ^ below, a 
part of new randomness at each step also comes from independent random variables Vi,, Vjv 


(having g-geometric or Bernoulli distribution, cf. Remark 2.10). Moreover, for > 0 the bijective 


mechanisms (4.1), (4.2) will be q-randomized (i.e. will no longer be deterministic bijections). 

and the 


[a 


This would lead to four q-randomized RSK correspondences: the row and column 
row and column (/3). In fact, for g > 0 the step-by-step nature of the g = 0 case (when push 
or pull operations are performed one at a time) will be broken, and certain series of push or pull 
operations will be clumped together and g-randomized as a whole. This will make the dynamics 
at the g-Whittaker level more complicated. 

Each of the four RSK-type dynamics possesses a marginally Markovian projection (onto the 
leftmost or the rightmost particles of the interlacing array) leading to a certain discrete time 
particle system on Z. Namely, Q?ciw[a] and Q?(jw[/3] give rise to the geometric and Bernoulli 


PushTASEPs, respectively, on the rightmost particles j = I,..., N. Similarly, 

geometric and Bernoulli TASEPs, respectively, on the leftmost particles 


The g-deformed dynamics of 
systems. 


and ^ below would lead to g-deformations of these four particle 


Remark 4.9. By imposing some reasonable nearest neighbor constraints on discrete time mul¬ 
tivariate dynamics, one may seek a full classification of solutions of the main equations of Theo¬ 
rem 2.13 in the Schur (g = 0) case. Such classification in continuous time setting was obtained 


in m- We do not pursue this direction here. 


and QlJ 


ADDING A DUAL PARAMETER 


5. RSK-TYPE DYNAMICS Q?o\ 

In this section we explain the construction of two RSK-type dynamics on g-Whittaker processes 


adding a dual parameter jd to the specialization (in the sense of (2.17)). For g = 0, these dynamics 


degenerate to (/3) dynamics on Schur processes arising from row and column RSK insertion. We 
also discover that for 0 < g < 1, the row and column dynamics Q?ow[/3] and Q^oJ/^] related 
by a certain transformation (we call it complementation). Moreover, in a small fd limit the 
complementation provides a direct connection between continuous time RSK-type dynamics on 
g-Whittaker processes introduced in |59] (column version) and [E] (row version). 


5.1. Row insertion dynamics Qrow[/3]. Let us now describe one time step A —?■ iz of the mul¬ 
tivariate Markov dynamics Qrow[/3] on g-Whittaker processes of depth N. A part of randomness 
during this step comes from independent Bernoulli random variables Vi,... G {0,1} with 


^^Starting multivariate dynamics from initial condition = 0 for all 1 < i < j < N and considering all 


levels of an interlacing array, bijections (4.11 and (4.21 extend to bijective correspondences between certain integer 


matrices and pairs of semistandard Young tableaux (in agreement with the well-known understanding of RSK 
correspondences). 
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parameters /3ai,... ,l3aN, respectively (these random variables are resampled during each time 
step). 

The bottom most p article of the interlacing array is updated as + Vi (as it should 

be, cf. Remark 2.10). Next, sequentially for each j = 2,... ,N, given the movement A —)• i/ at 
level j — 1, we will randomly update A —>■ at level j. To describe this update, write 

i-i 

v — \ = Cje,, Ci G { 0 ,1}, e, are basis vectors of length j — 1, 

^=1 

and say that numbers {k,m), where l<k<m<j — 1, form island{k,m) if 

Cfc_i = 0 (or A: = 1), Cfc = Cfc+i = ... = Cm = 1, and Cm+i = 0 (orm = j-l). 

That is, all particles that have moved at level j — 1 split into several disjoint islands. Also denote 
for any i = 1,j — 1: 


U = fj(p, A) := 


l-q 




gi = gi{v,\) ■.= l-q 




I _ qUi-i-ui+1 ’ 

(by agreement, let Pq := +oo). Note that all these quantities are between 0 and 1. 


(5.1) 


2 + 1 3 


3 + 1 5 + 1 

5 + 1 

T + Vj 

• 

• 



• • 


• 



• A + - A) 








+1 


+1 


+1 


+ 1 




\ 


7 - 


\ 


-< 

\ 



•. 



• 





\ + {v-X) 


2 + 1 3 + 1 

5 

5 + 1 

6 + 1 


2 + 1 3 + 1 



5 5 


6 + 1 

8 + y,- +1 

• 

• 



• • 


• 



• A + - A) 






_ 1 




+1 +1 





+1 


+1 




s 



/ 

7 - 





•. 

• 


• 




• 

A + (z^ — A) 


2 + 1 3 + 1 

5 

5 + 1 

6 + 1 



Figure 17. An example of a step of Qrow[/3] at levels 5 and 6 . There are two islands, 
( 1 , 2 ) and (4,5), moving at level j — 1. Above: Vj = 1, and the probability of the 
displayed transition is 1 • (1 — f4)g5 = 1 — q (note that here the particle A4 = 3 cannot 
be chosen not to move because f 4 = 0). Below: Vj = 0, and the probability of the 
displayed transition is (1 — fi)(l — g 2 ) • f 4 = (note that here the particle A 4 = 5 must 
be chosen not to move because f 4 = 1 ). 


The update A —)■ u at level j goes as follows (see Fig 
to the right by Vj, i.e., = Ai + Vjei. Then, independently for every island(A:,m 


17). First, the rightmost particle jumps 

of particles 


that have moved at level j — 1 , perform the following updates: 


(1) If Vj- = 1 and A; = 1 (i.e., the particle Ai has already moved, and the island contains the 
first particle at level j — 1 ), then move the particles A 2 ,..., Am+i at level j to the right 
by one with probability 1 . 

( 2 ) liVj = 1 and A: > 1, or Ij- = 0 (i.e., island(A;,m) does not interfere with the movement of 
Ai coming from Vj, or there is no independent movement of Ai), then island(A:, m) triggers 
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the movement (to the right by one) of all particles Afc, 
which does not move is chosen at random: 

• Afc is chosen not to move with probability 


, Am +1 except one. The particle 


h = 


I — qi^k-i-i'k+i ’ 


each As, k + 1 < s < m, is chosen not to move with probability 


(1 - ffc)(l - gfc+i) ... (1 - gs-l)g5 = 


gAj,— + 1 

1 _ + l 


(^l — q 


As —l^s+1 


• Am +1 is chosen not to move with probability 
(1 - ffc)(l - gfc+l) ... (1 - gm.-l)(l - gm) = 


qXk-f>k+l _ 


1 _ qi^k-i-i^k+i 


(5.2) 

); 

(5.3) 

(5.4) 


Probabilities (5.2), (5.3), and (5.4) are nonnegative, and their sum telescopes to 1. 

This completes the description of the (/3) row insertion RSK-type dynamics Q?ow[/3]. Clearly, thus 
defined conditional probabilities Uj, j = 1,..., A^, for this dynamics satisfy (2.24[). 


Remark 5.1. The g'-deformed probabilities (5.2), (5.3), and (5.4) ensure that mandatory pushing 
and blocking mechanisms (built into Definitions 4.2 and 4.3) work automatically: 

• If A<j = — 1 for any k < s < m, then the particle A^ cannot be chosen not to move. 

This agrees with the mandatory pushing of A^ by the move of A^ = A* which is necessary 
to restore the interlacing. 

• If Afc = Uk-i (i.e., Afc is blocked), then ffc = 1, so Afc must be chosen not to move. This 
means that in this dynamics no move donations ever arise (cf Remark 4.4). 


Theorem 5.2. The dynamics Qrow[/3] defined above satisfies the main equations (2.26), and hence 
preserves the class of q-Whittaker processes and adds a new dual parameter fi to the specialization 
A as in (2.17). 

Proof. We need to prove ( |2.26 ) for any fixed j = 2,... ,N and A, G GT^, u G GT+_^, where 
(cf. Fig. 0 right). For a subset I C {1, 2,..., j — 1}, set 

Uj := (1 + fiaj)Uj{X —7- I A —)• i >)-— , , 

where \ = v — ^ ^ GT^_^ is obtained from v by shifting back (by one) all particles 

with indices belonging to I. By agreement, if I is such that A does not satisfy A ^ 

(cf. Fig. right), then Uj = 0. With this notation, the desired identity (2.26) turns into 

= 1. (5.5) 


Note that the denominator (1 + fiafi coming from the Bernoulli distribution of Vj will always 
cancel the corresponding factor in all 17/’s. 

First, let us consider a particular case when u = A + XI^fc ezi be., the movement A —)■ z/ involves 
a consecutive group of particles from k to m, where 1 <k<m<j — 1. There are four subcases: 
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1. If A; > 1 and m < j, then necessarily Vj = 0, and (5.5) becomes 

m—1 


^[k — 1,771 —1] ^ ^ ^[/c—l,s—l]u[s+l,m] ^ (5.6) 

s=k 

(here and below by [fc — 1, m — 1], etc., we mean the corresponding interval of indices). See Fig. 18 

Am + 1 ... As+i + 1 As + 1 ... Afc + 1 Afc_i 


'^m+l 


• # 

Am +1 Am + 1 


^ ^ ^ ^ ^ 


As +1 + 1 


As 


#■ 

As-i + 1 


Afc-i + 1 


Figure 18. Situation corresponding to the s-th term in (5.6). 


Using (2.4), (2.7), we have (as before, here and below in the proof we agree that uq = +oo) 


U[k-l,m-l] ffc-l(u, A) - 


/Am-Am + l\ ( Afc_i^ . ,7 _,7 4_1 

V \m-Vm ) q ^Afc_i-I/fe_i+l/^ 7 — q '““2 fc-l + 


/Am,H-l /^k—1 _ q 




V’A/A/b^/i^ 

\ gAfc_i—i^fc-1 + 1 \ gAm —^m + 1 _ g^k-l—^k \ g^fc-2 — 1+ 1 

1 — q^k — 2~^k — l~^^ \ — gAm —Attt,_|_i + 1 2 gA/c —— 1 + 1 \ q^k — l~^k 

I — m — i^m + 1 

m—ATTT.-1-1 + 1 

Also for any fc < s < m — 1, 

^[fc—l,s—l]U[s+l,m] 

= ffc_i(u. A) • (1 - f5+i(u, A))(l - gs+2(u, A))... (1 - gm(u, A)) 


(A™-Pr+l)g {\-iT)q ^Xk-l-Pk-l + l\ (l-g^fe-2-^fc-l + l)(l-+-^-+l + l) 


/Arn, + 1 —Att^-I-iA /As —As_|_i\ /Afc_i—A; c~1 \ 

^ Atti + I —f'rn 'Q \As + 1 —Ps'' O ^A^_i —Pfc_i/^ s. 


]_ _ Q^k — 1 ^k 


q \-^k-l — ^k-lq 

'-s/-' 

bA/A/V’^/i^ 

^ ^Afc_l —Pfc_l+1 ^As-^l —Ps-^l + 1 gPs—Ps-l-l+1 


gEIl.+ 2 (^-G + l) 


\ — g^^fc- 2 —^fc -1 + 1 \ — i^s + 1+1 

_ q^'f^ A/Ti-i-i _g,As+i I's _^A;j._i A;, _ q^k—2 1+1^^]^_^s+i+i^ 


— gf A 771 + 1 Atti-I-I As-j-1 — gA;^ —1 —1+1 

^ ^1 _ g^m~A77i-^i > 


1 _ gAfc_i—Afc 


1 - q^ m + 1 —A 777 -I-I 




and 


U[k,m] = (1 - ffc(u, A))(l - gfc+i(u. A))... (1 - gm(u. A)) 
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f + ( k — 1 ^ \ - - I 1 


+ l —1 \ — q^k — 1 ^k 

V Atti + I —t^rn, //7 V Ai._i —r'L. — l ' 


gXk—i>k-\-^ _ q^k—i—^k-\-^ 

_ q^m Am + 1 


-0^ ,t/'0^ m 
^u}X' ^vjX 

Arn,-|-l _ ^Afc — i Afc t'fc + l 


gEIl.+i(A-P,+l) 


X — grAm + 1 A^jt,-!-! — q^k — 1 ^k — q^k — 


^h—^ —At. 


\ — grATTi + l —Am+l 




The summation in (|5.6|) thus telescopes and gives 1 as desired (similarly to the sum of expressions 
d^ , d^ , and d^y 

2. If A: > 1 and m = j, then also necessarily Vj = 0, and there is only one I, namely, [A;—1, j —1], 


contributing to (5.5). We have 


( A(j_i—Afe ■x 

W-i-i>k-i+iL I - q‘'k-2-i'k-i+i 


U[k-ij-i] = ffc-i(z^, A) 


At._i —Afc 


/A/c—1 ^k r\ 

^ A(j_i—Pfe_i / „ 


\ — gAfe_i—Afe _ qU],_2 — Uk-l + l 


Q 


1 — qk'k-2 — ^k-l + ^ \ _ Q\k-l—l'k-l + ^ \ — Q^k-l—^k 


= 1 , 


q 


so we see that (5.5) holds. 

3. If A: = 1 and m < j, then Vj can be either 0 or 1, and (5.5) now looks as 

m— 1 

Un ,m] ■E f^[l,s-l]U[s+l ,m] + {ajf3) ^U[i 

,m—1] 1* 

s=l 

This identity is established similarly to the subcase 1. Namely, one readily sees that 

I — Am+i i/.+i'i 

Ur, 1 = - - - + 

I _ q\„^-\m + l + l^ 




_ qi^m Am + l^^j^ _ + 


1-q^ m Am + 1+1 

1 — m — i^m + 1 


,E™.+i(E-P,+l). 


g" 

and the sum of these quantities telescopes and gives 1. 

4. If A; = I and m = j, this means that necessarily Vj = 1, and the only term that enters (5.5) 
is = Paj, so the desired identity also holds. 

We have now established the desired identity in the particular case y = A + 
general case there could be several consecutive groups of particles forming the move A —)• p at 
level j. Let there be gaps of at least two not moving particles between neighboring moving groups. 


Then, by the product nature of the quantities V' and ip' (2.4), (2.7), as well as by the independence 


of propagation for different islands at level j — 1, cf. Fig. 17, the sum in the left-hand side of (5.5) 


can clearly be represented as a product of sums corresponding to individual groups of moving 
particles. Each such individual sum is the same as in one of the subcases 1-4 above, and therefore 


is equal to I. This implies (5.5) in the case when moving groups at level j are sufficiently far 
apart. 
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1 ... ... ... Afc + 1 



Am-1 + 1 As + 1 As-1 + 1 Afc + 1 

Figure 19. Two islands at level j corresponding to a single island at level j — 1. 



As As_i + 1 

Afc + 1 

• • >. 

• .. 

.• • 

As As—1 + 1 

Afc + 1 


Figure 20. Two configurations giving the same contribution as the one on Fig. 19 


Finally, it remains to check (5.5) in the case when there could be moving groups at level j 


separated by one not moving particle. Consider two such neighboring groups. The only configu¬ 
ration of moves at level j — 1 (corresponding to these two groups at level j) that could prevent 


the sum in (5.5) to be of product form is given on Fig. 19 However, one readily sees that the 


contribution of this configuration is the same as the product of contributions of two conhgura- 


tions on Fig. 20 Indeed, factors involving the quantities are already in a product form, and 
the remaining factors (coming from Uj and the quantities '!/>') are 

1 — Q^k-l — ^k 


Afc _ Afc \ \ \ \ 


1 - q 


Afc_i —\k 


{1-q 


Ai;-1—Aj 


)(1 — gAa-As+l) 


(1 - 0(1 - gfe+i)... (1 - gs-i)gs on Fig.[T^ 


1 - q 


As—As 


l-q 


As—1 As 


q 


At.—Az. 


1 —At. 


l-q 


Ac_i As 


l-q 


As Ac 


+1 


1 — q^>^- 


—At 


-9 




1-0 


— Q^k — l At 


U on Fig.[^ left (1 - ft){I - gt+i)... (1 - gs_i) on Fig.|2^ right 

Note that we have expressed everything in terms of signatures A and A because the signatures u 


differ on Fig. 19 and Fig. 20 


Therefore, in the last remaining case we can still rewrite (5.5) in a product form. This completes 
the proof of the theorem. □ 

Remark 5.3 (Schur degeneration). If one sets q = 0, then in a generic situation (when particles 
at levels j — 1 and j are sufficiently far apart from each other) all quantities b and gj become equal 
to one, see (5.1). One readily sees that the dynamics Q?ow[/3] reduces to the dynamics Qrow 


on 


Schur processes. The latter dynamics is based on the classical Robinson-Schensted-Knuth row 


insertion (14.3) 


5.2. Bernoulli g-PushTASEP. One can readily check that under the dynamics Qmw [/3] we have 
just constructed, the rightmost N particles A) of the interlacing array evolve in a marginally 
Markovian manner (i.e., their evolution does not depend on the dynamics of the rest of the 
interlacing array). Namely, at each discrete time step t ^ t + 1 the bottommost particle is 
updated as -|- 1) = A^^^(t) -|- Vi, and for any j = 2,N: 
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(j—l) 

• If ' has not moved, then the rightmost particle at level j is updated as 

• If A^ ' has moved to the right by one, then the same particle is updated as 

+ 1) = + Vj + (1 - Vj) ■ ^0-1), 

where pushing by happens with probability 1 — fi = which depends 

only on the rightmost particles of the array. 

(Recall that the V^’s are independent Bernoulli random variables which are independently resam¬ 
pled each step of the discrete time.) This evolution of the rightmost particles A^'^^ 1 < j < .A^, 
leads to a new interacting particle system on Z which we call the {discrete time) Bernoulli 
q-PushTASEP. We discuss this process in detail in ^ below. 


+1 +1 


+1 +1 


+1 +1 

• 4P 4P 


+1 

• ^ 


+1 +1 


/ vr/// >-r/ / .r/// 

« V. >• . i . 4 « m' V. >• • # # V. i . 4 « 


+1 +1 +1 


+1 


+1 +1 


+1 +1 


+1 


Figure 21. Complementation of propagation rules turning the dynamics Qrow[/3] (with 
move propagation given by thin solid arrows) into (corresponding to thick dashed 

arrows). 


5.3. Complementation. Let us take another look at propagation rules employed in the d efini - 
tion of the row insertion dynamics Qrow[/7] on g-Whittaker processes (see the beginning of 15.1). 
These rules state that, generically, an island of moving particles at level j — l splits (at random) 
into two moving islands at level j separated by exactly one staying particle (either of two moving 
islands at level j is allowed to be empty). Now consider the pattern of staying particles at levels 
j — l and j. We see that an island(A:, m) (where k < m) of staying particles at level j — l always 
gives rise to an island(fc -|- l,m) of staying particles at level j, plus one more staying particle 
somewhere to the right of k (but to the left of the next staying particle at level j). The latter 
staying particle (whose index is chosen at random) is precisely the one separating the two moving 


islands at level j. See Fig. 21 


The transformation of propagation rules of Qrow[/3] that we just described informally in fact 
leads to a new RSK-type multivariate dynamics on Whittaker processes. Let us work in a more 
general setting: 


Definition 5.4 (Complementation of a dynamics). Assume that Q is a multivariate sequential 
update dynamics on g-Whittaker processes adding a specialization (/3). For j = 2,...,N and 
signatures A, u G GT^, A, u G GTj"_i satisfying conditions on Fig. right, let Uj{X ^ v \ A — )• u) 
be the corresponding conditional probabilities. Assume that the dynamics is translation invariant, 
i.e., that the values lIj(A — )• u | A —)• u) do not change if one adds the same number to all 
coordinates of all four signatures. 

For S a sufficiently large positive integer, define the complement conditional probabilities as 


U' (A ^ u I A ^ u) := (^[5 - A] ^ [S + 1 - u] 


[S' — A] —>• [S -I- 1 — , 
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where 


[S - A] := {S-Xj>S- Xj-i > ... > S - Ai) 

is the complement of the Young diagram A in the j x S rectangle, and similarly for [S' + 1 — i^], 
[S —A], and [S + 1 —p] (hence the name “complementation”). Note that these four new signatures 
also satisfy conditions on Fig. right. 

Let us denote by Q' the dynamics on interlacing arrays corresponding to IX'-, j = 2,..., iV. Note 
that due to translation invariance, the complement dynamics does not depend on S provided 
that S is large enough. 


Lemma 5.5. Let S be sufficiently large. For /r G /r G GT^ such that g, have 

'^[S-u]/[S-p.] ~ 

For X G GTt such that /r -<v x, we have 


'^'[s+i-x\/[s-F\ ~ 




Proof. A straightforward verification using dehnitions (2.4), (2.7). 


□ 


Proposition 5.6. If Q is a multivariate sequential update dynamics adding a specialization {jl), 
then so is the complement dynamics Q'. 


Proof. One can show that the complement dynamics Q' satishes the same main equations (2.26) as 
the original dynamics Q. Indeed, Lemma 


5.5 


ensures 


that the coefficients '^v/v'^'y/x 


do not change under complementation, and powers of (aj/3) also transform as they should: 
= ([5 - A] ^ [5 + 1 - I/] [S - X] ^ [S + I - i?' 


= [a 




— A] — [5+1—— [5—A] + [5+1— 


u. 


— A] —>■ [S' + 1 — [S — A] —[S + 1 — . 


This establishes the main equations for the complement dynamics. 


□ 


5.4. Column insertion dynamics 


coiiYj. Clearly, the row insertion dynamics Q?ow[/3] on q- 
Whittaker processes is translation invariant (in the sense of Definition 5.4), so one can define the 
complement dynamics. Denote it by describe (in an explicit way) the evolution of 


the interlacing array under this new dynamics during one step of the discrete time. See Fig. 22 
for an example. 

As before, a part of randomness comes from independent Bernoulli random variables Vj G {0,1} 
with P{Vj = 0) = 1/(1 + /doj), j = 1,..., N. The bottommost particle of the interlacing array is 


updated as = A^^^ + Vi. Sequentially for each j = 2,..., N, given the movement A —)• at 

level j — 1, we will randomly update A —)• at level j. Let us denote (as usual, T'o = +oo) 


1 _ 

= f/(p,A) :=-4- 

; I — qi'k-i-i'k+i' 


g's = g'siLffi)--=l-q^^-^-^P (5.7) 


The update X ^ v looks as follows: 

(1) Consider a pair of moved particles (Ar,Afc) at level j — 1, where 0<r</c<j — 1, 
such that the particles ++i,..., A^-i in between did not move (by agreement, r = 0 
corresponds to Xk being the rightmost moved particle at level j — 1). Regardless of the 
value of Vj , each such pair of moved particles at level j — I triggers the move (to the right 
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Figure 22 . An example of a step of levels 5 and 6 . Above: Vj = 0, and the 

probability of the displayed transition is 1 — fg = (g + + q + q^)- Below: Vj = 1, 

and the probability of the displayed transition is (1 — gg)(l — fg) = cp'. Note that in 
the latter case the particle A 3 = 6 cannot be chosen to move because it is blocked by 
A 2 = A 3 which is not moving; this agrees with fg = 0 . 


by one) of exactly one particle A^, r + 1 < s < k, between them at level j. If r + 1 = k, 
then there is only one choice s = k, so Xk must move. Otherwise, the moving particle As 
is chosen at random (independently of everything else) with the following probabilities: 

• If s = k, then As is chosen to move with probability 


f/ ^ ^-q’' 

Ifr + l<s<A:, then As is chosen to move with probability 

Vk-l—^k 


(1 - fD(l - gfc-l) ■ • • (1 - gs+l)gs = 


q 


1 _ qVk-i-Vk+^ 


qUk-i-Uk+i 

ql^i^ 


—Ai+l) _ ^Ua-\—\s 


If s = r + 1, then As is chosen to move with probability 


(1 - fD(l - gLi) • • • (1 - gr+3)(l - gr+ 2 ) = 




—Ai. 




\ — qi^k-i-Vk+1 


(5.8) 


(5.9) 


qT.’f=r+iVi-V+i) _ ( 5 . 10 ) 


Clearly, these probabilities are nonnegative, and their sum telescopes to 1. 

(2) ff Vj = 1, then in addition to the moves described above, exactly one more particle at 
level j is chosen to move (to the right by one). Namely, let Am be the leftmost moved 
particle at level j — 1. If m = j — 1, then the additional moving particle at level j is Xj, the 
leftmost particle. If m < j — 1, then one of the particles As with m+1 < s < j is randomly 
chosen to move (independently of everything else) with the following probabilities: 

• If s = j, then As is chosen to move with probability 


g' 

• Ifm + l<s<j, then As is chosen to move with probability 

(1 - gjxi - si-i)... {1 - e;+i)g; = (1 - /--"■), 

• If s = m + 1, then As is chosen to move with probability 

(1 - gjxi - si-i). ■ ■ (1 - e;„+3)(i - g:„+2) = 


-Ai+i). 

1 


■Ai+i) 


(5.11) 

(5.12) 


(5.13) 
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The sum of these probabilities also telescopes to 1. 

This completes the description of the (/3) column insertion RSK-type dynamics 


Theorem 5.7. The dynamics Q^oJ/^] defined above satisfies the main equations (2.26), and hence 
preserves the class of q-Whittaker processes and adds a new dual parameter fi to the specialization 
A as in (2.17). 


Then the desired state- 
□ 


Proof. One can readily check that Q^oJ/^] is the complement of Qmw 
ment follows from Theorem |5.2| and Proposition |5.6[ 

Remark 5.8. Similarly to Qrow[fi] (cf. Remark 5.1), probabilities (5.8)-(5.13) employed in the 
definition of Q^oJ/^] ensure the mandatory pushing, blocking, and move donation mechanisms 
(described in Definitions 4.2 and 4.3 and Remark 4.4). Namely, observe that 

• If Afc = Afc for some k and A^ has moved at level j — 1, then = 1, which means that A^ 
is chosen to move with probability 1. 

• If As = As-i, and As-i has not moved, then = f( = 0, so according to (5.8), (5.9) the 
particle As at level j cannot be chosen to move. If, moreover, As has moved at level j — 1, 
then this move will trigger some other particle to the right of As at level j to move. In 
other words, the moving impulse coming from As —Ps = As + I will be donated further 
to the right of As- 


Remark 5.9 (Schur degeneration). When ^ = 0, one readily sees from (5.7) that generically (i.e., 
when particles at levels j — 1 and j are sufficiently far apart) we have R = gs = 1- This implies 


that the dynamics degenerates to the multivariate dynamics on Schur processes. 

The latter is based on the classical Robinson-Schensted-Knuth column insertion (14.3). 


5.5. Bernoulli g-TASEP. Under the dynamics 


coU 


the leftmost N particles A^^ of the 
interlacing array evolve in a marginally Markovian manner. Indeed, one can readily check that at 
each discrete time step t —)■ t + 1 the bottommost particle is updated as A^^^(^ + 1) = A^^^ {t) + Vi, 
and for any j = 2,..., N: 

(j—l) 

• If A)_^ ' has moved, then the leftmost particle at level j is updated as 


(j—l) 

If A:_i ' has not moved, then the same particle is updated as 


X^/\t + l) = X^f{t) + Vf, 


i-1 


X^p{t + l) = X^f>{t) + Vj-l, 


- xOT 


^ ^ is chosen to move’ 


where Xj is chosen to move with probability gj = 1 — q ^ ^ ^ ' which depends only 

on the leftmost particles of the array. 

This evolution of the leftmost particles A^^ 1 < j < N, is the {discrete time) Bernoulli q-TASEP 
which was introduced and studied in [8] . 

5.6. Small fi continuous time limit. If one sends the parameter fi to zero and simultaneously 
rescales time from discrete to continuous, then both dynamics Qtow[/ 3] and turn into certain 

continuous time Markov dynamics on q-Whittaker processes. At the level j = 1 (cf. Remark 2.10), 
this limit transition coincides with the one bringing the (one-sided) discrete time random walk to 
the continuous time Poisson process. In continuous time setting, at most one particle can move 
at each level j = 1 ,..., A during an instance of continuous time. 
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The continuous time limit of Qtow[/ 3] looks as follows. Each rightmost particle of the 
interlacing array has an independent exponential clock with rate Uj. When the clock rings, the 
particle jumps to the right by one. 

(j—l) 

There is also a jump propagation mechanism present: If at level j — l some particle Am ' has 
moved (to the right by one), then this move instantaneously triggers the move of the upper left 


(?) 

neighbor with probability fm = 


X — q'^m 

1-q 


21 


or the move of the upper right neighbor 


Am^ with the complementary probability 1 — f^. This dynamics was introduced in [16] (Dynamics 
8). Under it, the rightmost particles of the array also evolve in a marginally Markovian manner. 
This leads to the continuous time q-PushTASEP on Z [T6l §8.3], [23] . 

The continuous time limit of looks as follows. Each particle Afc, 1 < /c < j, at level j 

has an independent exponential clock with rate 


ajgj, 

ajil - g')(l - ... (1 - g'k+i)g'k, 

[a,(l-g')(l-g'_i)...(l-g§)(l-g'2), 


k = j; 

1 <k <j] 
k = l. 


These quantities correspond to (5.11 )-(5.13) with u = X (because if an independent jump occurs 


at level j then at level j — l there could be no movement). When the clock of Xk rings, this 
particle jumps to the right by one. Note that the move donation mechanism described in Remark 


4.4 follows from the above probabilities. 

There is also a jump propagation mechanism: If a particle Xk has moved at level j — 1, then it 
triggers the move (to the right by one) of exactly one particle A*, 1 < s < A:, at level j, where s 
is chosen at random with probabilities 


fL 


ki 


(l-f')(l-g'_,). 


• •(i-gUi)gs> 

..(i-g^)(i-g^). 


s = k] 

1 < s < k; 
s = 1. 


The above probabilities are given by (5.8)-(5.10) where n differs from A as = A + e^. This 

Under it, the leftmost particles of 


dynamics on g-Whittaker processes was introduced in 
the interlacing array evolve in a marginally Markovian manner as a g-TASEP. This continuous 
time particle system was introduced in [7]. See also, e.g., m, [in], m for further results on the 
g-TASEP. 

Thus, the two continuous time dynamics on g-Whittaker processes (or, in other words, q- 
randomized Robinson-Schensted correspondences) introduced in [59] and [16] are the /3 -A 0 
degenerations of and Qtow[/ 3], respectively. On the other hand, complementation (15.3) 

provides a straightforward link between the two latter discrete time dynamics. 

6 . RSK-TYPE dynamics Qrow[a] and Q^oJo] adding a usual PARAMETER 
In this section we explain the construction of two RSK-type dynamics Qrow]®] and oa 


g-Whittaker processes adding a usual parameter a to the specialization (as in (2.17)). Por q = 0, 


^^Note that here we are rewriting fm through the particle coordinates before the move at level j — l (cf. (5.11) 
so that the probabilistic meaning is clearer. One could also similarly rewrite all other formulas for fm,gm,fm,gm 
above in this section, but for the purposes of checking the main equations of Theorem |2.13| it is more convenient 
to use the notation involving the coordinates P after the jump. 
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these dynamics degenerate to (a) dynamics on Schur processes arising from row and column RSK 
insertion. As in the case of Qtow[/ 3] and dynamics, in a small a limit the dynamics Q?ow[a] 

and QpoJ®] degenerate to continuous time RSK-type dynamics from j59] (column version) and 
m (row version). 

6.1. The g-deformed Beta-binomial distribution. We will use the following quantities: 

Definition 6.1. Let y G {0,1, 2, ...} U {-|-oo}, and s G {0,1, ..., y}. Recall the g-notation from 
ra. Let 


I y) ■= C 


iv/(\Q)si^-,q)y-s {q;q) 


iv,q)y 

If y = -|-oo, the limits of the above quantities are 

,{v/(-,q)s i^-^q) 


{q‘,q)s{q;q)y-s 


I +Oo) = e" 


( 6 . 1 ) 


( 6 . 2 ) 


{q-,q)s {v,q)oo' 

An important property of the quantities (6.1) and (|6.2[) is that for all y G {0,1, 2,...} U {+oo}, 


we have 


I 2/) = 1- 


(6.3) 


s=0 


This statement may be rewritten as the g'-Chu-Vandermonde identity for the basic hypergeometric 
series 201 • For the proof and more details see [36], [21]. Recall that in general the unilateral basic 
hypergeometric series j4>k is defined via 


7 0fc 


ai 

bi 


-,q,z 


oo 


(® 1 , ■ ■ ■ i CLj] q)n 

^ (6i,.. .,bk,q]q)r 


((-!)”«<»>) 


1+k-j 


(6.4) 


where (ci,..., Cm] q)n = YYiLiici] q)n- Later on in this section to prove some identities we will 
need to apply transformation formulas for certain hypergeometric series. 

Therefore, for all values of the parameters (g,^,y) for which \ y) is well-dehned and 

nonnegative for every 0 < y < s, (6.1) defines a probability distribution on {0,1,...,y}. One 

I . Another choice of parameters 


such family of parameters is0 <g<l,0<r/<^<l, cf. 
leading to a probability distribution which we will use is ^q-i^qo.^qb{- \ c), where a < b, c < b are 
nonnegative integers. 

The distribution appears (under a simple change of parameters) as the orthogonality 

weight of the classical g-Hahn orthogonal polynomials |46| §3.6], and is also related to a very 
natural g-deformation of the Polya urn scheme unj. As such, be regarded as a g- 

deformed Beta-binomial distribution, since the latter is the orthogonality weight for the Hahn 
orthogonal polynomials, and also arises from the ordinary Polya urn scheme. We can also directly 
see by taking g = e~^,^ = letting e —)> 0-|-, that (pq^^ q{s \ y) converges to 

r(a -h /3) r(Q; + y- s)T{/3 + s) /y'' 


r(a)r(/3) r(a + ^ + y) 

which is the probability of s under the beta-binomial distribution with parameters y, a, (3. 
Let us now record two straightforward observations which we will be using below. First, 


n 

kjq-. 




(6.5) 
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Second, if o < 6, c < 6 are nonnegative integers (6 might also be +oo), then for any s G {0,1,..., c} 
one has 


Indeed, in this case 




('5 I ls=max{c—a,0} ■ 


^ (g^;g-^)c W, 


( 6 . 6 ) 


If a > c, as g —)• 0 this converges to 1 for s = 0 and to 0 for s > 0. If a < c, as <7 —)• 0 this 
converges to 0 for 0 < s < c — a, since (g'“; q~^)c-s vanishes, to 0 for s > c — a, since a positive 
power of q tends to 0 , and to 1 for s = c — a. 

6.2. Row insertion dynamics Qrow[a]- Let us now describe one time step A — 1 / of the multi¬ 
variate Markov dynamics Qrow[a] on g- Whittaker processes of depth N. A part of randomness a 
time step comes from independent g-geometric random variables Vi,..., G Z>o with parame¬ 
ters aai,..., auN, respectively (these random variables are resampled during each time step). 

The bottommost particle of the interlacing array is updated as -|- Vi. Next, sequen¬ 
tially for each j = 2,..., N, given the movement A = AL“^1 —> p = at level j — 1, we will 

randomly update A = A^^ u = at level j. To describe this update, write 


P - A = y^^Cjej, 

i=l 


Ci G Z>o, Qi are basis vectors of length j — 1. 


Note that by interlacing, it must be that Cj < Aj_i — Aj. 

Sample independent random variables ITi,..., ITj_i, such that each Wi G {0,1,..., Cj} is 
distributed according to 


I Ci), where ^i := and r^i := 


(6.7) 


(this is a probability distribution because Xi — \i < Ai_i — Ai and c* < Ai_i — Ai, cf. 16.1). We 
will use the conventions Aq = -|-oo and 7/1 = 0. Define a sequence of signatures 

A = 77(0),/x(l),... ,p(j - 1) 


Via 


H{i) := n{i - 1) -F Wj-iej-i + {cj-i - Wj-i)ej-i+i for 1 < i < j - 1 

(where ej are basis vectors of length j). Finally, define := p(j — 1) -|- VjCi, this is our new 
signature at level j. 

In words, each zth particle on the (j — l)-st level which has moved by Cj, must trigger a total 


of Ci moves (to the right by one) at level j (this the RSK-type property, see Definition 3.2). Each 


such particle at level j — 1 independently from the others, in parallel, splits contribution from 


its jump between its nearest neighbors on the level j, according to the distribution (6.7). After 


this pushing, the rightmost particle on the j-th level additionally performs an independent jump 
according to the g-geometric distribution with parameter aaj. Clearly, thus defined conditional 


probabilities Uj, j = I,..., N, for this dynamics are nonnegative and satisfy (2.24). See Fig. 


23 


One must verify that the interlacing properties (as on Fig. left) are preserved by this dy¬ 


namics: 


Lemma 6.2. If X A -<h A and Uj(X —)• z/ | A —>■ z^) > 0, then z/ z^ and X z^. 
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Proof. Observe that for a <b and c < b 

Pq-i^qa^gb{s I c) = 0, if s > 5 — a or c — s > a. (6-8) 

Apply this for a = Xi — Xi, b = Aj_i — Xi, c = Ci to get Ci — Xi + Xi <Wi < Aj_i — Aj. 

Since Vi = Xi + Wi + Cj-i — Wj-i, we have 

^ Xi A^—1 A^ A Ci—\ Wi—i — A^—1 A Ci—i 1 — ^i—1 hhi—1 ^ ^i—1 

and = Aj A Ci < Ai A Wj < r',, so v Ah i'- Moreover, we can also write 

t'i < Aj A Ai_i — Ai A Ai_i — Ai_i = Ai_i, A* < Vi 
which implies that Ah A. □ 

This verification completes the description of the (a) row insertion RSK-type dynamics Q?ow[a]. 

Remark 6.3. (Schur degeneration). If one sets q = 0, then the dynamics Qrow[a] reduces to 
the dynamics on Schur processes based on the classical Robinson-Schensted-Knuth row 

insertion (14.3). To see this, observe that (6.6) implies 

I ls=max{ci-Ai+Ai,o} as g —)• 0, 

that is, each Wi becomes equal to max{ci — A* A Ai,0} in the q \ 0 limit. Therefore, the 
update A —is reduced to applying Ci operations pull at positions i from j — 1 to 1, plus an 
additional independent jump of the rightmost particle according to the geometric distribution 
with parameter aoj. 


0 + 2 


6 + 1 


9 + 2 


15 + 5 


25 + 5 35 + 1 + V,' 




2 + 3 


12 + 4 


19 + 5 


29 + 4 


A + (i^ - A) 

A + (z5 - A) 


Figure 23. An 
6, with Vi = 


example of a step of Qrowla] at levels 5 and 

3. The probability of this update is equal to 

{oiaep 


I ^7 (2 | 4)(,Pg-l ^6 glO (2 | 5)<,Pq-l ^6 q(1 | 4)(q;£16, ( 7)00 (^q-^q)^ ' hlote 

that, e.g., ,j7(0 | 4) = 0, which ensures the mandatory pushing (by at least 1) of 

As by the move of A3. 


Theorem 6.4. The dynamics Qrowla] defined above satisfies the main equations (2.25), and hence 
preserves the class of q-Whittaker processes and adds a new usual parameter a to the specialization 


A as in (2.17). 


(cf. Remark 2.10). 


Proof. We will prove (2.25) by induction on j. Case j = 1 is straightforward because A is empty 


Assume now that (2.25) holds for signatures A, u having length j — 1, and let us prove this 


identity for A, z/ or length j. The idea is to expand each term in the sum in the left-hand side 


of (2.25) with respect to what happens to the leftmost particle on the {j — l)-st level (and its 


neighborhood), and then use the inductive assumption and the fact that the ip's sum to 1. 

For a signature = (/ii > ... > //m) we denote by fj,~ the signature (/ii > ... > fim-i) obtained 
by deleting the smallest part of /r, and by A [s]im the signature {fii > ... > Am-i > hm A s) 
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obtained by adding s to the smallest part of /j, (for s < /Jm-i — Mm)- To simplify certain notations 
below, also denote 


Vj(A —>• p I A —)• z^) := 'lXo(A —)• I A —)■ v) - — -r-. (6-9) 

{aaj-,q)oc 


Temporarily let t stand for Cj-i = z?j_i — Aj_i which is the move of the leftmost particle on 
the (j — l)-st level. In order to have at least one nonzero summand in the left-hand side of (2.25), 
we need to have (see Fig. 24): 


• t > Vj — Xj, since the jump of the leftmost particle on the j-th level happens due to 
contribution of a part of the jump of the leftmost particle on the {j — l)-st level. 

• t > Vj — Xj + Vj-i — Aj_i, since + Xj \ t) > 0 implies by 

— Aj < Aj_i — Uj-i +1. 

• t < Uj-i — Xj, since we must have Aj_i > Xj. 

• t < Vj-i — Aj_i -|- Vj — Xj, since the contribution from the jump of the leftmost particle 
on the {j — l)-st level is split between particles Xj and Xj-i at the level j. 


(6.8) that 


Denote the interval of t satisfying the above inequalities by I. We also must have 


t < Xj -2 — Xj-i + Vj — Xj, since + Aj | t) > 0 implies by (6.8) that 

t — Vj + Xj < Xj -2 — Aj_i. For j = 2 this last inequality should be omitted. 


We will use the notation A := A“. Denote by J{t) the set of signatures A of length j — 2, such 
that A -<h z^, A A“, and Xj -2 > t -|- Xj-i — Vj + Xj. For j = 2 this set consists of just the empty 
signature 0. If A is such that A Pj A A and 'Vj(A —)• z/ | A —)• z^) 7^ 0, then A“ G Uze/ 


3 


A. 


3 - 1 



t 


Aj-i 

o 



Aj-2 ^ 3-2 


Ai 


Ai 

o 


Figure 24. We expand sum with respect to the jump t = Cj-i = Mj-i — Aj_i of the 
leftmost particle on the (j — l)-st level. Note that the signatures v, X, v are fixed, while 
the positions of particles A* vary in the sum. 


The left-hand side of (2.25) divided by the right-hand side of the same equation is equal to 


EE 

is/ AeJ(t)^ 




^3 ^33 g-l 


A,-2-A,-i . „-l 


{q 


\j-2-yj-i+t. „-l 


g )t Uj+Xj 


Q )t 


‘P 




X Vj_i(A + [t ~ zAjf -|- Aj]ini —y 12 I A —>12 ) • 


/p-'/’i/-/A--|-Aj]ini 
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q \ uj-\j )q W ■' " jy 



t&I 


E(a-., 

^ ^ Vj_i(A + ~ ~\~ Aj]im —^ ^ I A —y ly ) • 


xeJ(t) 


^A -\-[t—Uj-\-Xj]i^/X^i' /X 
/p- 4^1/-/X- + [t—Uj+Xj]ini 




tel 


Above Vj_i and Vj have the same value of the parameter a = aj. We have also used the fact 
that 


hence 'Vj{\ 
X- 


|i^| — |A| — \iy\ |A| — \iy | — |A \~X t^j — Xj — \iy | A |A | — ^j—i Aj_i 
= \iy I ~ |A [t ~ Ajjiml ~ I “I" |A I, 

A- I A —)• z?) involves the same power of aaj as 'Vj_i(A“ + [t — fj + Aj]im 


V ). Also, ( 6 . 8 ) implies that ip^-i 


-i-'i 




only for t ^ I, hence one gets 1 after summing these quantities over t ^ I. 

This concludes the proof, and also establishes Theorem 1.3 from Introduction. 


IS nonzero 


□ 


6.3. Geometric g-PushTASEP. Under the dynamics Qrow[a] we have just constructed, the 
rightmost N particles A^'^^ of the interlacing array evolve in a marginally Markovian manner (i.e., 
their evolution does not depend on the dynamics of the rest of the interlacing array). Namely, at 
each discrete time step t —)■ t + 1 the bottommost particle is updated as A^^^ (^ + 1) = A^^^ {t) + Vi, 
and for any j = 2 ,..., A if we let gapj(t) = X^^\t) — A^'^ be the gap between the rightmost 
particles on the (j — l)-st and the j-th levels at time t, then 

X<f\t + l)=X[^\t) + Vq+Wj,t 

for an independent random variable Wj^t distributed according to 

The random variable Vj (recall that it has the g-geometric distribution with parameter aaj which 
is resampled during each time step) represents an independent jump of A^'^^ The variable Wjj 

(j) (i—l) 

represents the pushing of A) ' by the move of A) A 

This evolution of the rightmost particles A) , 1 < j < A, leads to a new interacting particle 
system on Z which we call the {discrete time) geometric q-PushTASEP. 


6.4. Column insertion dynamics Description and discussion. Let us now describe 

one time step A —)> zz of the multivariate Markov dynamics Q^oiC®] Whittaker processes of 
depth A. As in the previous case, the bottommost particle of the interlacing array is updated as 
z/j^^ = A^^^ + X for a g-geometric random variable X with parameter aai. Next, sequentially for 
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each j = 2,..., N, given the movement A —)• at level j — 1, we will randomly update A —)• at 

level j. To describe this update we write, as usual, 

j-i 

z/ — A = ^ ^ CjCj, Cj G Z>o- 

i=l 

All randomness during this update comes from a collection of 3j dependent random variables 
Xi,..., Xj,Yi,... ,Yj, Zi,..., Zj (they are resampled during each time step), and 

voluntary jump push from push from the “stabilization fund” 

(It will be convenient to let i represent the position of the particle counted from the left.) Observe 
that Yi must be identically zero. The “stabilization fund” means the leftover push from the first 
i — 2 particles from the left at level j — 1 (i.e., from Aj_i,..., Aj_i+ 2 ) (in particular, Z\ and Z^ 
are identically zero). 

Let us hrst formally define the distribution of all the parts of the jumps: 

(1) Set 9i := 1. For i from 1 to j sample Xi according to 


and set 


Xi ~ ‘Pq^aajdi,oi' I ^j-i -^j-i+l) 


9i+i := 9iq^^-‘ ^ 3 -^+^ 


( 6 . 10 ) 


The jump Xi comes from the input Vj, see Remark 6.7 below. Here the convention 
Ao = +00 applies when i = j. 

(2) Set Yi := 0. For i from 2 to j — 1 take Yi = y with probability 

~ “ ^j-l + l ~ ~ U \ ^j-i ~ ^j-i+l ~ ^i)- (6-11) 

Finally, set Yj := ci. 

(3) Set ri = r 2 = 1 and Zi = Z^ := 0. Set := For i from 3 to j — 1 take Zi = z 

with probability 


Yi 

and set 


i+1 h) Z \ Xj—i Aj_i-|_i Xi Yi^ 


( 6 . 12 ) 


n+i := riq''^-^+^ 


ld=s. 


Finally, let Zj := log^r^. 

Remark 6.5. For fixed s,u,d> 0 (possibly u = oo) and D ^ oo observe that 

_ J I (s-d)(D-d)/ s -U (g^+^ ^)D-d 

‘Pq-^,q»,q^+o[D d \ D) q [Q , Q 

Therefore, the definitions of Zj = log^ rj and Yj = ci are consistent with the definitions of Yj 
and Ij (i < j), respectively. In words, the consistency for Zj means that the stabilization fund is 
depleted for the push of the rightmost particle on the j-th level. The consistency for Yj means 
that the whole value of the jump of the rightmost particle on the (j — l)-st level is transferred to 
the rightmost particle on the j-th level via immediate pushing. 


Lemma 6.6. If X zz, A -<h A and Uj(A —)• z/ | A —>■ z/) > 0, then z/ zz and A z^- 
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0+1 3 + 2 + 1 8 + 1 + 3 + 2 17 + 3 + 4 + 1 30 + 0 + 4 + 3 


2 + 4'. 7 + 5'. 



15 + 5'. 


E 


E 


28 + 4 


stabilization fund 


Figure 25. An example of a step of Q^oiC®] levels 4 and 5. 


Proof. It is straightfor-ward from the definition of the dynamics Q^oJ®] that Vj-i+i < Xj-i < 


min(7'j_j, Xj-i). Also for 2 < i < j (6.11) together ■with ( 6 . 8 ) implies that Xj-i — Xj-i^i — Xi — Yi < 
Xj-i — Xj-i+i — Cj-i+i, hence Uj-i+i > Aj_i+i + Xi + Yi > It follo'ws that the interlacing 

properties are preserved. □ 

In the rest of this subsection we will describe the column insertion dynamics in words, and also 
discuss its various properties. The (rather involved) proof that this dynamics acts on ( 7 -Whittaker 
processes in a desired way is postponed to the next subsection. 

There are two stages of the update of particle positions Xj, Xj-i, ..., Ai, performed in order 
from left to right, which we will describe below. 

During the first stage of the update, the particles at level j level make voluntary jumps in 
order from left to right. The value A* of the voluntary jump of Xj+i-i depends on the previous 
jump Xi-i, where 2 < i < j. Indeed, this dependence comes from the parameters 9i (note that 
they are nonincreasing in i), see (6.10). Note that unlike the Qrow[ct] case, in which all random 
movements not coming from pushing are restricted to the right edge, in the case of any 

particle might make a voluntary jump. 

Remark 6.7. The random variable Xi + .. .+Xj has the g-geometric distribution with parameter 


auj, as it should be by Remark 2.10 and the discussion of { 3.2 This is seen by applying inductively 
the next Lemma. 

Lemma 6 . 8 . Let A and B he random variables such that A is distributed according to oi' I 
and B given A is distributed according to I (where b might he +ooj. Then A + B 

is distributed according to Pq^a,oi' I 0 , + b). 

Proof. Indeed, we have 

y 

Prob{A + B = y) = Prob{A = s)Prob{B = y — s| A = s) 

s=0 

= a^(a; q)a-s q),-y+s ( 

g—Q V'®/ a A 

y 

= '^ay{a-,q)a+b-y(f 


y - s 


,a(y-s) (g°;g ).(g ;9 )y-s ^-s[y-s) f y 


s=0 


(g;g) 


y 
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— (ot , q) a+b— 


a + b 

y 




'^+b{y - s\y) 


<? s=0 


= ^q,afl{y \a + b), 

which establishes the desired statement. 


□ 


The second stage of the update consists of pushing, in order from left to right. We start an 
initially empty stabilization fund, which will collect impulses not immediately used for pushing, 
and will be a source of the pushes Zi. The value of the stabilization fund just before the movement 
of Xj+i-i is logg ri (by agreement, ri = r 2 = 1 always). For each i ranging from 2 to j, the 
following three steps happen: 


(1) The particle Xj+i-i gets a push Tj fro m its lower left neighbor The size 

of this push (distributed according to ( 6 . 11 )) is at most Cj-i+i. 

(2) Then Aj+i-j gets a push from the stabilization fund (if it is not empty) of size 
not exceeding th e curr ent value of the stabilization fund. This push is distributed 
according to Zi ( 6 . 12 ). 

(3) Finally, the amount of pushing not used in (1) above, i.e., — Yi, is added 

to the stabilization fund. 


(6.13) 


One can also think that the above two update stages are performed together for each particle 

Xj, Ajr' — l, ..., Ai. 

Proposition 6.9. One can switch the order of the lower left neighbor pushing and stabilization 
fund pushing (i.e., steps (1) and (2) in (6.13)J without changing the dynamicsf^ 

Proof. Fix k = 2,... ,j. Suppose that after the voluntary displacement stage the distance from 
the /c-th particle from the left at level j (denote this particle by P) to Aj_|_i_fc is h := Xj-k ~ 
Aj_fc_|_i — Also set I := — Aj_fc+i, b := Aj_fc — Xj^k+h and let the current size of 

the stabilization fund be R. 

If the steps (1) and (2) in (6.13) are not interchanged, then the probability that P jumps by 
s > 0 is 

S 

‘Pq-kqt,qi>{h “ V \ h)(pg -1 ^gR^Q{h - s \ h - y). 

y=0 


If the steps (1) and (2) in (6.13) are interchanged, then the same probability is given by 

Y Pq-hq^,oi^ -S + y\ h)<.Pg-i gi gb{h - S | ^ - S T I/). 
y=0 

q(R+e)(h-s)^gb-i.q-l~^^_^ A' 


After dividing each of these two expressions by 
following identity we need to verify: 

. 1 —n Vy/ Q-I 


{q^',q~Oh-s 


C). 


we arrive to the 


^^Here for the version with interchanged steps (1) and (2) the distributions of the jump components are changed 
as Zi ~ “ ^j—i+i ~ ~ ' I Aj—i — Aj—— Xi) and Yi ~ “ Aj—i+i — 

Xi Zi • I Xj—i Aj—i+i Xi Zi)> 
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E 

y=0 


yj q-^ 


q^^{q^;q %{q^]q ^)s-y{q'' ’"^^;q)s-y (6.14) 


We are very grateful to Christian Krattenthaler for providing us with a proof of the g-binomial 
identity (6.14), which we reproduce below. 

First, use a transformation formula for 3(^)2 series |36l (III.12)]: 


3</'2 

Sending 6 —)■ 0 we obtain 
3</>2 


q-^Ac 

d, e 


]q,q 


{elc]q)r 

{e]q)n 




q "',c,dlb . bq 


1—n !^ 1 


d^cq^ ^/e 


r g-",o,c 1 

(e c;q)n „ , 

q '^jC dq 

i d,e 

- . . C 292 

{e-,q)n 

d,cq^-^le 


(6.15) 


Multiply both sides of ( 6.15[ ) by c '^{d; q)n{e', q)n to obtain 

{q~'^-^q)y 


E 


E 


,=0 

( 9 "”; 9 )t/(e/c; q) 


{c-q)y{dq^ ;q )n-yieq'^ ]q )n-yq^ = 


y^oAq)yicq^ ”/e;Q)j/ 
This equality can be rewritten as 

n 

^ / n\ 

c 

\ 1! f 

2/=0 


(c; q)y{dq^-^-q-X-y{-^fq^^^-^'^'\dq/e)y. 


^ r cy-^{c-^-q-\{dq^-^;q-X.y{eq^-^-q-^\ 


n-y 


E 

y=0 


n 


VJ q-^ 


{e-/c-,q)n-y{c ^)y{dq'^ ^]q ^)n-y{dq'^ 


Now make the substitution n := s, d := q^~^^ c := q e := q^~^^ ^ ^ to arrive to (6.14). □ 

Remark 6.10. (Schur degeneration) If one sets q = 0, then the dynamics Q^oJ®] reduces to the 
dynamics Q^^[a] on Schur processes based on the classical Robinson-Schensted-Knuth column 
insertion ( ^4.3[ ). Indeed, observe that 

lim q^q,uqt,o{s I g) = 1^=0 for t > 0 , and liin \ g) = (I-u + uls=g)A. 

q—>-0 ^ ^ q—>-0 ^ 

Thus, the first update stage (voluntary movements) reduces to the propagation of the impulse 
the leftmost particle receives (which has geometric distribution with parameter aaj). The lower 
left neighbor pushing due to ( 6 . 6 ) and the stabilization fund pushing together degenerate to 
performing Cj_i + • • • + ci operations push (Definition |4.3[ ) in order from left to right. 

6.5. Column insertion dynamics Proof. 


Theorem 6.11. The dynamics defined above satisfies the main equations (2.25), and 

hence preserves the elass of q-Whittaker processes and adds a new usual parameter a to the 
specialization A as in (2.17). 


Proof. We aim to prove the desired statement by induction on j. To apply this induction, we 
will need a more general statement. To describe it, introduce the following notation. For a 
nonnegative integer h, use ^(^(A —)■ | A —)> iz) to denote the probability that transition A —)• iz 
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on the (j — l)-st level spurs a transition A —)■ on the j-th level according to the rules of 
specified above, but modified so that Z 2 = z with probability 

“ ^jr-i — X 2 — Y 2 — z \ \j -2 — Aj-i — X 2 — Y 2 ), r 2 := q^. 

Note that the original dynamics r 2 = 1. In other words, the modification 11 j means 

that we introduce an additional impulse of size h which is distributed among particles at level j 
(except for \j), as if coming from (nonexistent) particles preceding the leftmost particle on the 
(j — l)-st level. 

Let a := |z/“| —|A“| —li^l + IAI (recall that the notation ijl~ means /x without the last coordinate). 
Under the modified probabilities as above, u — /i is a sum of voluntary movements of particles 
on the j-th level except for the leftmost one. Note also that —?■ zx | A —)• P) = 0 for h > a. 


For the purposes of the proof, let (see Fig. 26) 


a := Xj, 
c .— Aj_i, 


k := izj — Xj, 


b := Vi- 


1) 


t . 12 ■j — 1 At — 1 


d := 


;— iyj—\ At— 


j-i, 


J-2 
X ■.= X 2 


s iyj -2 — Xj- 2 , 

y ■= Y 2 . 


For a nonnegative integer H define 

H 

'Lf|^(A —)• IX I A —)• z^) := 


h=0 


^(H-h)a+h(b-t-a-k)y^h^^ ^ | A ^ P). 


(6.16) 


In particular, l(^(A —)• zx | A —)• P) = ^.^(A —)• zx | A —)• P). In general, the quantities are not 
probability distributions in zx. Their only meaning is that they come up in the inductive proof 
below. 

With all the above notation we are now able to describe and prove the generalized statement 
which we will prove by induction: 


E 

AeGT+_^ 


zx I A 


'i^vlv4>vlX 


= 1 


for any H > 0. 


(6.17) 


Here and below Vf is related to Uf as in (6.9). For H = 0 this statement gives us (2.25). 


For j = 1 we have cr = 0, so only the term h = 0 contributes to (6.16). Therefore, checking 
this induction base is the same as in the proof for Qrow[a] dynamics. 


a a + k 


J 

impulse h ^ 


c + X c + X + y c + 


Ai 


T — t b 


d — s d Ai 

^ # ••• a 


j -1 

Figure 26. We expand sum with respect to jump t = Cj-i = Pj_i — Aj_i of the 
leftmost particle on the {j — l)-st level, voluntary movement x of the second leftmost 
particle on the j-th level and push y from the leftmost particle on the (j — l)-st level. 


Let us now perform the inductive step. Denote by I the range of {t, x, y, h), such that 
t,x,y,h>0, x + y<i, h + t + x — ^>t), t<b — a — k. 
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Then we may write (see Fig. 




AeGT+_^ 


E 


'^vIv^vIX 

{Q;q)H {q-,q)c-a 


{ q ^ q ) b — a — k{qj q ' jc + i—b { q ^ q ' jd — s — b+t 


(q-, q)h{q] q)H-h (q; q)b-t-a(q; q)c-b+t (q; q)c+£-a-k (q; q)t(q; q)d-s-b 

{q!q)d—s—c {q-i q)d—s—c—x 


{t,x,y,h)&I 

„ {q',q)k{q-,q)c-a-k {q-,q)b-t-a 


{ q'l q ) c—a { q'l q ) k { q : q ) b — t — a—k { q ' lq ') x {, q ] q ) d — s — c—x { q ] q ) y { q ] q ) d — s — c — x—y 
{q iq ^y{q ^q }d—s—c—x—y {q^q)d—s—c—x—y 


{q i q )£—x—y 


X q 
X q ‘ 


H+B 

E 

r=0 


^qd s {^q] q^£_x—y{Q] Q)d—s—c—i 

—h{H—h)—y{d—s—c—x—y) — {d—s—c—i){i—x—y)+t(d—s—c—x—y) 


h{d—s—c—i)-\-(H—h)a-\-h(b—t—a—k)-\-(b—t—a—k)x-\-(b—t—a—k-\-i—x){(7—x—h) 


E 


AeGT 


_ I T _X ’^X-l'X^u-IX 

(A \X^iy )- - , 

Wv-lv-9v- jx- 


<h+t-\-x—£ / \ — 


3-2 


H + B 


{H-\-B—r)(a—£-\-t)+r{d—s—c—£) 




X E 

AeGT+_2 


_ I T _s '^X-/X^u-/X 


Vf^^{X- \X^v-) 


AeGT' 


7-2 


^v-lv-'i>v-jX- 

'^X- jX^v- jX 
1 v-(\>v-IX- 


= 1 . 


Here we have applied Proposition 6.12 (see below) with A = H, B = b — a — k, C = c — b+l, where 


r := h+ t — ^ + X is the value of the stabilization fund just before the push of the third leftmost 
particle on the j-th. level plus the value of the additional impulse in the inductive assumption. 
This completes the inductive step in proving (6.17), and thus implies the theorem. □ 


Proposition 6.12. For A, B, C,i,r > 0, such that A + B > r and B + C > i, one has 

B i i—x 


EEE 

t=0 x=0 y=0 _ 




q )£—x— 


{q-,q)r 


y I 


g-x \ “ / g-x (qi q)r+£—i 

{q-, q)A (g^+'®“'’; q~^)B-t+£-x{q'^^^\q~^)£{q^; q-^)t-x-y 


{q]q)A+B 


{qB+C.q-l)^(^qC+t.q-l)^_^ 

t(£—x—y)+{r+£—x){B—t) + {A+B—r)x 


Here and thereafter we use q-multinomial notation 


X q 


n 


= 1. 


{q]q)r 


m, kj ' {q-,q)kiq;q)m{q;q) n—m—k 
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We are extremely grateful to Christian Krattenthaler for providing us with a proof of this 
proposition. We reproduce the proof below. 

Proof. The left hand side of the equality can be expressed as a power series in , hence 

we can set a = q^, fi = q"^, j = q^ and prove a more general equality: 


B £ £—x 

EEE 

t=0 x=0 y=0 _ 


B\ . * ^)t+£-x-y 

(aq^lP; q-^)B-t+i-x{'yq^; q~^)e{r, q~^)£-x- 


{aq^;q ^)B{iq^-,q ^)t(7g*;g ^)i-x 
By first summing over y, the left hand side can be written as 




= 1 . 


B I 

EE 

t=0 x=0 . 


34>2 


q-\0,q-^+^ 


B\ {l3q^-;q-^)t+,_, 
<?■ 




(«g //^;g )B-w-x{iq]q )eir,q )e-x xob-i-xbi 
[aq^] g“^)B(7g^; g“^)r(7g*; g“^)r-a; 


We now apply transformation formula (6.15) to rewrite this as 
B i 


EE 

t=0 x=0 


202 


^-£+x 

H 5 y 


;g,/3gi+'-*-77 


/3gi-0g-V7 

{aq^ 115] q-'^)B-t+£-x{iq^] g“^)r(7; q~^)£- 


B\ %+i-x {iq]q)t 

®A-iVvg-i ( 0 g^"^;g"^)r-a; ( 7 g^"^+"';g)t 


X ^x oB—t—x^B£—£t-\-tx 


{aq^-,q ^)s(7g^;g ^)e{'yq^;q ^)e- 


X a^/3 


B t rmn{t,£—x} 

EE E 

t=0 x=0 y=0 


(g *;g)y(g ^''■,Q)y 


B\ {iq-,q)t 

(g“V7;g)y g)y(g;g)y q-A^)q-i (7g^“^+"=;g)t 


(0;g ^)tiaq /P;q ^)B-t+e-xhQ ;q )K7;g )£- 

(ag^; g-^)R(7g^; q~^)eW; q~^)e-x 


X 


B I 

EE 

f=0 j/=0 . 


101 


-£+y 


aq 


q 

1 - 


;g,ag^^* 


(-10 


(g ^;g)t/(g *;g)y 


(g ^h\(i)y{^(£ ^■,Q)yAQ)y 

B\ (0; qAtiaq^/P; g~^)R-t+H7g*; g~^)^ Rr-tt+£j/+y2/2+?;/2-tt/ 

7,-. (a?®; 9 -')b{t(®; <!-*)( P 1 <! 

The last equality is obtained by summing over x. We now use the summation formula j36l (II.5)]: 


101 


a 

c 


■,q,c/a 


and by summing over y rewrite our expression as 


{c/a;q)oo 

(c;g)oo 
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B 

E 

t=o L 


3</>2 


q I3q ^/a,q * i+f , 




(/3;g %{aq^ll3-,q ^)e ^ oB-tm-et 


t) a-^ {mVl3\q ^)e{aq^;q ^)Bi'yq^]q ^)e 


X 


We now aim to use the transformation formula |36i (III.13)]: 

(e/c; q) 


3(p2 


^ ■,q,deq^/bc 


{e-,q)r 


i4>2 


<1 ■ 
d.cq^-'^je 


(6.18) 


Applying it, we can rewrite our expression as 
B 


E 

t=o L 


302 


q \q \otq 


B\ (i/7;g)K0;g )tiaq //3;q ^)B-t+£{'yq*-,q )i _ 

* / 0-1 («9V0; q~^)B{iq^; q~^)£{q~^/i-, q)eiq; q)i 




X 


B e 

EE 

i=0 j/=0 . 
t B 

EE 

j/=0 t=y L 


/ f-' •) '1 I'l j JD \ ! 'i 7'i 

(g~^; g)y(g~^ g)y(«g; g)i;(7; g~^)^-;/(0g^; g~^)t(«g^/0; g~^)R-f+£ .s-t 
i y g-i (agV0; 9"^)^(a9^; q~^)B{iq^] q-^)f;{/3qy-, q-^)y{q; q)y 

B-y\ 0^-*(/3; q-^)t.y{aq^jfd- q-^)B-t ^ fi\ q^^^-^Kr, q-^)i-yiq^-, q'^ 

X 


t-y d q-l 


{aq^;q ^)B-y 


y/q-^ 


{iq^\q ^)t 

The fact that this expression is equal to 1 now follows by applying (|6.3l) twice. □ 


6.6. Geometric g-TASEP. Under the dynamics leftmost N particles of the 

interlacing array evolve in a marginally Markovian manner. 

Namely, let gapj(t) := Aj'y^^^(t) — {t) be the gap between the consecutive leftmost particles 

at time t. We assume gapi(t) = +oo. Then at each discrete time step t —)• t + 1 the leftmost 
particle on the j-th level is updated as 

X^\t + i) = \B){t) + Wj,t 

for an independent random variable Wj^t distributed according to <^q^aajfi{' I g3Pj(0)- 

This evolution of , 1 < J < A^, is the {discrete time) geometric q-TASEP which was 
introduced and studied in [H]. 


6.7. Small a continuous time limit. Let us send the parameter a to zero and simultaneously 
rescale time from discrete to continuous. Namely, set a := (1 — q')A, and let each discrete time 
step correspond to continuous time A. In the limit A —)• 0, both dynamics Qrowfa] and Q^oJo] 
turn into the same continuous time Markov dynamics on g-Whittaker processes as in §5.6| above. 
That is, the limit of Qrow[a] is the dynamics introduced in m. same as for Qrow[/S]- The limit 
of Q^oJ®] is the dynamics introduced in [59], same as for Q^oJ/^j- To see this, note that repeated 
g-geometric trials can be approximated by (continuous time) Poisson processes in this scaling. 
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7. Moments and Fredholm determinants 


In this section we briefly discuss moment and Fredholm determinantal formulas for the Bernoulli 
g-PushTASEP started from the step initial configuration (corresponding to a['^^(0) = 0, j = 
1,... ,N). The Fredholm determinantal formula which we extract from moment formulas in a 
way similar to |12j allow us to prove (in a small (3 continuous time limit, cf. ^5.6[ ) the conjectural 
Fredholm determinantal formula |23( Conjecture 1.4] for the continuous time g-PushASEP (which 
is a two-sided dynamics unifying continuous time g-TASEP and g-PushTASEP). 


7.1. Bernoulli g-PushTASEP on the line. In this section it will be convenient to work in the 
shifted coordinates 


:= - i, 




so that xi > ... > xn- We will think that the Xj’s encode positions of particles on the line Z 
which jump to the left. Let us reformulate the dehnition of the Bernoulli g-PushTASEP (15.2) in 
these terms. 


Definition 7.1. Each discrete time step t ^ t + 1 of the Bernoulli (^-PushTASEP consists of the 
following sequential updates (see Pig. 27): 

(1) The first particle xi jumps to the left by one with probability and stays put with 

the complementary probability 

(2) Sequentially for j = 2,..., N-. 

(a) If the particle Xj-i has not jumped, then Xj jumps to the left by one with probability 

, and stays put with the complementary probability 

(b) If the particle Xj-i has jumped (to the left by one), then Xj jumps to the left by 


with probability ^ a —, and stays put with the complementary probability 


one 


l+aj/S 


, where gap j{t) := Xj-i{t) — Xj{t) — 1 is the number of holes between the 
particles before the jump of 


(1) Update Xi{t -\- 1) first: 


Prob = 


aiP 

1 + 01/9 




(2a), (2b) Then update X 2 {t + 1) based on whether xi has jumped: 


Prob = 

r\ 


0.2 P 

1 + 02/3 


Prob = 

r\ 


oiP+q^ 

1 + 02/3 


Q # QOOOOQ # QO > 0 # 000000 # 00 > 0 # 00000 # 000 > 

X2{t) xi{t) X2{t) a;i(t-|-l) X2{t) a;i(t-|-l) 


Eigure 27. Bernoulli g-PushTASEP (on this picture, gap 2 (t) = 6). 

We will assume that the Bernoulli g-PushTASEP starts from the step initial configuration 
a^i(O) = -i, i = 

Remark 7.2. Similarly to [H], one could readily make the parameter f3 of the process depend on 
time, so that at each discrete time step t ^ t + I, a new parameter (3t+i is used. This will not 
affect the presence and the general structure of moment and Fredholm determinantal formulas. 
Below we will use a fixed parameter f3. 

^^Note that if Xj-i has jumped and Xj{t) = Xj-i{t) — 1, then the probability that Xj jumps is equal to one, as 
it should be. 
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7.2. Connection to the Bernoulli ( 7 -TASEP. The Bernoulli ( 7 -PushTASEP looks quite similar 


to the Bernoulli q-TASEP introduced in [ 8 ] (see also 15.5 above for an explanation of how the 
latter process arises from the dynamics g-Whittaker processes). 

Moreover, there exists a direct coupling between the two processes which we now explain. Recall 
that under the Bernoulli g-TASEP (we will denote its particles with tildes: xi(t) > ... > xjyit)) 
particles jump to the right by one according to the rules on Fig. 28 Let this process also start 
from the step initial configuration Xi(0) = —i, i = 1,... ,N. 


(i) Update ii(t + 1) first: 

p,.„K _ 0-10 

riou — 

r\ 

o m oooooo m oo > 

X2{t) Xi{t) 


(2a), (2b) Then update X 2 it + 1) based on whether xi has jumped: 
rroD - 


r\ 

o m oooooo m oo > 

X2{t) Xiit + l) 


P'* = its 

o » oooooo6 t 

X2{i) Xiit- 


- 1 ) 


Figure 28. Bernoulli g-TASEP (on this picture, gap 2 (t) = 6 ). 


Proposition 7.3. Let {xi(t)}i=o,i,... he the Bernoulli q-PushTASEP started from the step initial 
configuration and depending on parameters {oj} and fi. 

Then the evolution of the process {t + a;j(t)}t=o,i,,,. coineides with the Bernoulli q-TASEP 
{5i(i)}i=o,i,... started from the step initial eonfiguration and depending on the parameters {a~^} 
and 

Proof. The process {t + Xi{t)} jumps to the right, and, moreover, each of its particles makes a 
jump precisely when the corresponding gr-PushTASEP particle xfit) stays put. In particular, the 
first particle xi{t) stays put with probability 1/(1 + aifi) = (aj"^/3“^)/(l + Next, if 

xi{t) stayed put, then X 2 {t) stays put with probability 1/(1 + 02/3) = I3~^)/{1 + a^^/3“^). 

Otherwise, if xfit) jumped to the left, then X 2 {t) stays put with probability 

^ 02/3 + _ 1 - _ 02 

1 + 02/3 1 + 02/3 1 + O 2 ^ 

We see that the particles {t + Xi{t)} indeed perform the Bernoulli (;-TASEP evolution with the 
desired parameters. □ 


One can think that this coupling between the two particle systems on Z comes from the 
complementation procedure (15.3) relating the corresponding two-dimensional dynamics. 


7.3. Nested contour integral formulas for g-moments. The above coupling between the 
Bernoulli g-PushTASEP and the Bernoulli q-TASEP allows to readily write down moment for¬ 
mulas for the former process: 

Theorem 7.4. Let {xi{t)}t=o,i,... be the Bernoulli q-PushTASEP jumping to the left, started from 
the step initial configuration. Fix k > 1. For all t = 0, 1 , 2,. .. and all integers N > ni > n 2 > 

... > rifc > 0 , 
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fc(fc-l) 


(27ri) 


n 


l<A<B<k 


- ZB 

- qZB 


n n 

i=i V=i 
-1 


1 — aiZ 


-1 


l + q-^PzJ^' 


dz-i 


(7.1) 


where the contour of integration for each contains , • • ■, , and the contours {qzB}B>A, 

but not poles 0 or (—/3). 


□ 


Proof. Immediately follows from Proposition 7.3 and [U Theorem 2.1.(3)]. 

Remark 7.5. Since Xi{t) + i > —t for any i = 1, ..., and any t > 0, the g-moments in (7.1) 
admit an a priori bound. Therefore, for a fixed t > 0 they determine the distribution of the 
random variables (xi(t),... ,X]sf{t)). 


Remark 7.7. One can also establish the nested contour integral formula (7.1) directly, similarly 
to [8] (see also 


Remark 7.6. The moment formula (7.1) readily leads (similarly to [l2]) to a Fredholm determi- 
nantal formula for the g-Laplace transform of Xn{t) for any n. A similar Fredholm determinantal 
formula already appeared in [9] (based on the technique of Macdonald difference operators, cf. 
0 )- We discuss this and more general two-sided formulas in ^7.5| below (in particular, see Propo¬ 
sition 7.15 See also 17.4 for a related discussion of the case of the geometric g-PushTASEP. 

can also establisf 
). Indeed, denote 

/j(y) := (^d?/l+J/2-l-...+3/jv) Jgstep I 

'^i=0 

and, by agreement, the product is zero if yo > oQ One can directly 
show that these quantities satisfy certain linear equations in the y^-’s. For each i = 1,... ,N, 
consider the following difference operators acting on functions in y\ 

Vi 

[H'^'%f{y) ■= I yi)f{yo,yi, ■ ■ ■,yi- 2 ,yi-i + Si,yi - Si,yi+i,.. .,yN)- (7.2) 

Si=0 


N 


where (yo, ■■■,yN) e 


(7.3) 


Here the quantities are defined in (6.1). 

Also, denote by the operator which acts as [77'^’^]* in each variable yp. 

Applying operators in this order corresponds to first changing yi (by decreasing it by si), 

then y 2 (by sending S 2 to yi — si), etc., up to y^- In other words, these changes (encoded by 
si,..., Sat) happen in parallel, simultaneously with each of yi, y 2 , ■.., yN- 

One can then show that for any t = 0,1,2 ,... and any y = (yo, yi,..., yAr) G the 

quantities T(y) satisfy 

These linear equations can then be solved by the coordinate Bethe ansatz technique, because 
the action of each of the operators and reduces to the action of a free operator 

(i.e., which acts on each of the variables n* separately; note the identification of the yj’s and 
nfs in the previous footnote) plus two-body boundary conditions. This immediately leads to the 
desired nested contour integral formula. 


24, 


One should think that the variables yj encode the rii’s in (7.11: each yj denotes the number of Ui’s which are 


equal to j. 
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7.4. Remark. Geometric (^-PushTASEP formulas. There are also nested contour integral 
formulas for q-moments of the geometric (^-PushTASEP (^6.3). They can be obtained directly 
using the definition of the dynamics, similarly to the approach outlined in Remark |7.7[ The 


moment formulas (for the geometric g-PushTASEP jumping to the left) will have the same form 
as in (7.1), with the following replacement of factors: 





1 + q-^Pz 


1 


1 + /32 


-1 


n {l-aq-^z- 


However, because particles in the geometric g-PushTASEP can jump arbitrarily far to the left (at 
least as far as by independent g-geometric jumps), only a finite number of g-moments of the form 

jgstep ^ exists. Therefore, these (/-moments do not determine the distribution of 

the geometric g-PushTASEP. 

One can overcome this issue and write down a Fredholm determinantal formula for the dis¬ 
tribution of the geometric (/-PushTASEP via a certain analytic continuation from the Bernoulli 
case. This analytic continuation is performed in [9l Section 3] and heavily relies on properties 
of (/-Whittaker symmetric functions. Namely, [9l Theorem 3.3] contains a Fredholm determi¬ 
nantal expression for the distribution of the particle under a (/-Whittaker process. Our 
results on RSK-type dynamics (^ show that the same distribution arises under the geomet¬ 
ric g-PushTASEP started from the step initial configurationj^ thus [HI Theorem 3.3] provides a 
Fredholm determinantal formula for the geometric g-PushTASEP. 


7.5. Fredholm determinantal formula for the continuous time g'-PushASEP. The mo¬ 
ment formulas of Theorem 7.4 combined with certain spectral ideas allow us to establish Con¬ 
jecture 1.4 from |23j concerning a Fredholm determinantal formula for the continuous time q- 
PushASEP. 

For simplicity, from now on we assume that all parameters Uj = 1. Let us recall the dehnition 
of the g-PushASEP. 


Definition 7.8. The continuous time q-PushASEP x(r) depending on parameters L,R> 0 lives 
on particle configurations (x^r < < ... < xi) on Z and evolves in continuous time r as 

follows: 

• Each particle Xj jumps to the right by one at rate i?(l — agreement, 

xo = -1-00). 

• Each particle Xj jumps to the left by one at rate L. If a particle Xj has jumped to the left, 
then it has the possibility to push its left neighbor Xj+i with probability g^dT)-Xi+i(r)-i^ 
(Pushing means that Xj+i instantaneously moves to the left by one.) If the push¬ 
ing of Xj+i has occured, then Xj+i may push Xj +2 with the corresponding probability 

^Xi+i(r)-Xi+2('r)-l^ ancj so on. 

We are assuming that this process starts from the step initial configuration Xj(0) = —i. 


Remark 7.9. When P '\0 and one rescales the discrete time to the continuous one as ~ t//I, the 
Bernoulli g-TASEP (Fig. 28) clearly converges to L = 0, R = 1 version of the above process (this 
is the pure, one-sided g-TASEP). Under the same limit, the Bernoulli g-PushTASEP (Fig. 27) 
becomes the L = 1, R = 0 version (pure g-PushTASEP). 


other words, the geom etric g-PushTASEP provides a coupling of the quantities arising from g- 

Whittaker processes Ma (j 2.31 differing by adding usual parameters to the specialization A. 
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Theorem 7.10 f [23l Conjecture 1.4]). Let x„(r) he the n-th particle of the q-PushASEP started 
from the step initial configuration. For all C & C \ M>o, 


E 


L 


-) =det(/ + iv:^). 
) / 


Here det(/ + is the Fredholm determinant of K(^\ L?‘{Ci) —)• L?‘{Ci), where Ci is a small 
positively oriented circle containing 1, and K(^ is an integral operator with kernel 


A'c(w,u/) = ^ 


TT 


^ioo+1/2 

/-ioo+l/2 Sm^-TTSj 


(-O' 


{w]q)'^ q^w-w' 


ds. 


The L = 0 case of the above theorem corresponds to the (;-TASEP and was established in 
0 , m The i? = 0 case (one-sided g-PushTASEP) with different contours is contained in [9l 
Theorem 3.3], and it is obtained by an analytic continuation from the Bernoulli casej^ This 
analytic continuation relies on algebraic properties of g-Whittaker symmetric functions. The 
corresponding algebraic picture for the two-sided dynamics is not developed (cf. the discussion 
in [231 Appendix A]). To establish the above theorem we will utilize a certain two-sided process 
at the Bernoulli level whose /? \, 0 limit leads to the distribution of particles under the two-sided 
g-PushASEP at any given time r > 0. 


The rest of this subsection is devoted to the proof of the above theorem. 

Definition 7.11. Eix t G Z>o, and let = (X 7 v^(^) < •■■ < be a random 

vector which encodes positions of particles started from xl^^(O) = —i (i = which 

have evolved according to the Bernoulli g-TASEP (with right jumps) for [tR\ stepsj^ and then 
according to the Bernoulli g-PushTASEP (with left jumps) for [tLj steps. Both discrete time 
processes are now assumed to depend on the same parameter fd (also recall that Oj = 1 ). 


It is possible to write down g-moments of this random vector: 


Proposition 7.12. Fix k > 1. For all t = 0,1, 2,... and all N > ni > n 2 > ■ ■ ■ > Uk > 0, 



(-l)^g 


fc(fc-l) 

2 


(27ri)^ 


n 


l<A<B<fc 


- ZB 
^A - qZB 


where the contour of integration for each contains 1, and the contours {qzB}B>A, but not 
poles 0 or (—/3)^^. 


Proof. Let Tright and Tieft denote the one-step transition operators of the Bernoulli g-TASEP 
and the Bernoulli g-PushTASEP acting in x variables, respectively. Also denote 

N 

H(f,y) := 

j=0 

(by agreement, the product is zero if yo > 0 )- 

^®The contours used in [5] seem to be more suitable to the particular type of asymptotic analysis performed in 
that paper. Here we do not pursue further discussion of integration contours. 

^'^Here and below, [• • • J means the floor function. 
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Then it follows from the results of |8] and the discussion of Remark 7.7 that 


/ o transpose 


_i \ transpose 


ltL\ 

□ 


where the operators are defined in (7.2)-(7.3), and q~(y'^+"-+yN) jg the multiplication opera¬ 
tor. Observe that the operators 77^’^ applied to H do not change yi -|-... -|-yAr and hence commute 
with this multiplication operator. The inverse operators above exist on a certain space Wmax of 
functions in y, see pi] . 

Applying the Plancherel theory m (see also [231), we now see that the application of the 
product (Tright)^^^^ which is needed to compute the expectation in the left-hand side of 

(7.5), reduces to multiplication by the corresponding eigenvalue 0^=1 ) 

under the nested contour integral. This completes the proof. 

Remark 7.13. Similarly to Remark 7.5 one sees that the g-moments determine the distri¬ 
bution of the vector This implies that any powers of operators Tright and Tieft (defined 

above) commute when applied to the step initial configuration, i.e., they yield the same proba¬ 
bility distribution regardless of the order of application. 

Proposition 7.14. As /3 \ 0, the random variables X^\\j/13\) defined above converge in distri¬ 
bution to Xn(T) (the latter is the particle position coming from the continuous time q-PushASEP 
started from the step initial configuration). 

Proof. Let Gright denote the infinitesimal generator of the pure g-TASEP process (with only 
right jumps) corresponding to Definition 7.8 with L = 0, 77 = 1. Similarly, let Gieft denote the 
infinitesimal generator of the pure g'-PushTASEP (with only left jumps) corresponding to L = 1, 
77 = 0. 

R follows from Remark 7.13 that the semigroups and commute when applied 

to the step initial configuration (in the same sense as in Remark 7.13). This means that the 
semigroup of the two-sided g-PushASEP has the form 

^ right ^ left 

(when applied to the step initial configuration). Therefore, Xn^r) has the same distribution as if 
it arises by first running the pure Gright process for time t 77 , and then the pure Gieft proce ss fo r 
time tL. The latter two-part evolution clearly is the /3 \ 0 limit of Xn\\j / fi\) (cf. Remark 7.9), 
as desired. □ 


Proposition 7.15. Recall the random vector x^^\t) of Definition 


7.11 


E 


XCqxi^^P)+^-q), 


= det(7 +a: 


ihh 
c ^ 


For all C £ C \ M>o, 

(7.6) 


Here det(7 is the Fredholm determinant of : L^(C'i) —)■ L^(C'i), where Ci is a small 

positively oriented circle containing 1, and is an integral operator with kernel 


28 


We are writing “transpose” simply to indicate that the operators in the right-hand side act in variables y. 
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TT 


-i-cr 


^ ^ioo+l/2 

27ri y_ioo+i/2 sin(-7rs) 

{wq^; q)^ 1 + I3q^w\ I3{q^w)~^ ^ 

1+I3w J V ^+/3w 


(w;q) 


n 

oo 


,-l 


q^w — w' 


ds. (7.7) 


Proof. The passage from the moment formulas (7.5) to the desired Fredholm determinantal for¬ 
mula is done similarly to [12] (based on [7]). Namely, for |((| sufficiently small, the Fredholm 
determinant can be obtained by purely algebraic manipulations. Then one must show that the 
resulting right-hand side of ( |7.6[ ) is analytic in (( G C \ M>o, which follows from bounds like in 
[121 Proposition 3.6] and can be readily checked in our situation. The left-hand side of (7.6) is 
also analytic in ( because the function ( (-)■ (C;9)oo is uniformly bounded away from zero and 
analytic once C is bounded away from M>o. Therefore, the desired claim holds. □ 


Theorem 


7.10 


now follows by observing that the kernels converge, as /3 \ 0 and t = [r//3J, 


to the kernel from Theorem 7.10 because q^ and w are uniformly bounded on our contours. 
This also implies the convergence of the corresponding Fredholm determinants. On the other 
hand, by Proposition 7.14, we know that the left-hand sides of (7.6) converge to the left-hand 


side of (7.4). This establishes Theorem 7.10 


8. Polymer limits of (a) dynamics on q-Whittaker processes 

In this section we explain how the two (a) dynamics on g-Whittaker processes behave in 
the limit as q 1. This leads to discrete time stochastic processes related to geometric RSK 
correspondences and directed random polymers. 


8.1. Polymer partition functions. Let us first describe the polymer models we will be dealing 
with. They are based on inverse-Gamma random variables: 

Definition 8.1. A positive random variable X has Gamma distribution with shape parameter 
0 > 0 if it has probability density 

P{X G dx) = —^x^~^e~^dx. 

r(0) 

We abbreviate this by A ~ Gamma(0). Then X~^ has probability density 

P{X~^ G dx) = ^ x~^~^e~^^^dx, 

r(0) 

which is called inverse-Gamma distribution and denoted by Gamma”^(0). 


We recall partition functions of two models of log-Gamma polymers in 1 -|-1 dimensions studied 
previously in [ZQ], m. m, [60] . [58] . [25] (see also [56] for a continuous time version). Both 
models are defined on the lattice strip {{t,j) \ t G {0,1, 2,...}, j G {1, 2,..., n}}. One should 
think of t as time. Suppose we have two collections of real numbers Oj for j G {1, 2,..., n} and 
9t for t G {0,1,2,...}, such that Oj -|- 0* > 0 for all j and t. 


Definition 8.2 (Log-Gamma polymer [70]; Fig. 29 
equipped with a random weight dtj. 


left). Each vertex {t,j) in the strip is 
These weights are independent, and dtj is distributed 
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according to Gamma ^{9j + 6t). The log-Gamma polymer partition function with parameters 
9j,0t is given by 

(s,i)S7r 


where the sum is over directed up/right lattice paths vr from ( 1 , 1 ) to (t, j), which are made of hori¬ 
zontal edges (s, i) —)• i) and vertical edges (s, i) —)• (s, i-\-l). Extend this definition to denote 

by Rj^it) for t > k the weighted sum over all fe-tuples of nonintersecting up/right lattice paths 
starting from ( 1 , 1 ), ( 1 , 2 ),..., ( 1 , fc) and going respectively to {t, j — k1), {t, j — k + 2),..., 
{t,j). The weight of a tuple of paths is defined by taking a product of weights of vertices of 
these paths. The inequality t > k ensures that R^f^{t) is positive. 


Definition 8.3 (Strict-weak polymer 
strip is equipped with a random weight 


[58jl_jFig. 29 right). Each horizontal edge e in the 
. These weights are independent, and 


is distributed according to Gamma(0j -|- 9t). The strict-weak polymer partition function with 
parameters 9j , 9t is given by 


Riit) ■= Yl^e, 

7r:(0,l)—eSTT 


( 8 . 2 ) 


where the sum is over directed lattice paths from (0,1) to {t,j) which are made of horizontal 
edges (s,f) —)• (s -|- l,i) and diagonal moves {s,i) —)> (s -|- l,z -|- 1). The product is taken only 
over horizontal edges of the path. Extend this definition to denote by T],(f) for t > j — k the 
weighted sum over all fc-tuples of the corresponding nonintersecting lattice paths starting from 
(0,1), (0, 2),..., (0, k) and going respectively to {t,j — k -\- 1), {t,j — k -\- 2),..., {t,j). The weight 
of a tuple of paths is defined by taking a product of weights of horizontal edges of these paths. 
The inequality t > j — k ensures that L^j^{t) is positive. 


Distributions of ratios of the polymer partition functions defined above are sometimes called 
Whittaker processes (or, to be more precise, a-Whittaker processes), cf. [7]. They arise as limits 
(as q, the aj's and the afs simultaneously go to 1 ) of suitably rescaled particle pos itions in an 
interlacing integer array distributed according to the Whittaker process ((2.3), where 

A = {ai,... ,at), a = (oi,..., On). 


The convergence of ( 7 -Whittaker processes to Whittaker processes is known in the literature, see [3 
Thm. 4.2.4]. Since both dynamics Qrow[a] and Q^oJ'a] constructed in ^sample the ( 7 -Whittaker 
processes, we can employ them to give another proof of this limit transition. Moreover, we also 
establish the convergence of the corresponding stochastic dynamics. 

Let us first define the appropriately scaled pre-limit dynamics. In what follows, for e > 0 and 
9j,9t as above, we set q := e“^, aj = and at := 

Definition 8.4 (Scaled Qrow[a] dynamics). Start the dynamics Qrow[a] from the zero initial 
condition (that is, Ap^(O) = 0). Denote by rj^k{t,e) the position of the k-th particle from the 
right on the j-th level of the array after t steps of the dynamics (at each time step t —?■ f -|- 1 , 
apply the dynamics Qtow[o] with parameter a = at+i). For t > k, define the random variables 
Rl{t, e) via 

rj,k{t, e) = {t + j -2k + l)e~^ log e~^ -h log(Rl(t, e)). 


29, 


These two papers independently introduce essentially the same model. We will be using the notation of |25| 
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Figure 29. Paths and tuples of paths that contribute to the polymer partition func¬ 
tions: i?i(9) (top left), RI{9) (bottom left), F®( 8 ) (top right), ^ 3 ( 8 ) (bottom right). 


The reason for the restriction t > k comes from the fact that e) = 0 for t < k. 

We will view the collection of random variables e)} as a stochastic process R{t, e) which 

at a fixed time t becomes an array Rl.{t, e), 1 < k < j < n fov t > n, or a truncated array Rl.{t, e), 
1 < A: < min{t,j} < n for 0 < t < n. 

Definition 8.5 (Scaled dynamics). Start the dynamics from the zero initial con¬ 
dition, and denote by e) the position of the /c-th particle from the left on the j-th level of 

the array after t steps of the dynamics (again, at each time step t —)■ t -|- 1, apply the dynamics 
''^ith parameter a = ot+i). For t > j — k + 1, define the random variable T'(.(t, e) via 

^j,k{t,e) = {t-j + 2k- l)e"^log€"^ - e"Mog(L^(t, e)). 

The reason for the restriction t > j — k + 1 comes from the fact that e) = 0 for t < j — k + 1. 

We will view the collection of random variables {Lj^{t, e)} as a stochastic process L{t, e), which 
at a fixed time t becomes an array Ll,{t,e), 1 < k < j < n for t > n, or a truncated array Lj^{t, e), 
^ ^ k < j < min{n, k + t — 1} for 0 < t < n. 


Remark 8.6. Observe that for a fixed time t, the array rj^k{t, c) has the same distribution as the 
array £jj-k+i{t, e) (by Theorems 6.4 and 6.11, they are distributed as g-Whittaker processes). 


n 


Hence the (possibly truncated) arrays Rl.(t,e) and 1/Tj_^_,_]^(t, e) for 1 < A: < min{A,j} < 
have the same distribution. However, these arrays will not be identically distributed as stochastic 
processes in t since they come from different multivariate dynamics. 
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In the setting of polymer partition functions, define random processes R{t) and L{t) on (pos¬ 
sibly truncated) arrays via 

Ri{t) := Ri{t)/Ri_-^{t), for 1 < A: < min{t, j} < n 

and 

Li{t) := Li{t)/Li_-^{t), foT 1 < k < j < min{n, k + t-1}. 

They are well defined, because Ri.{t), R?j^_i{t) > 0 for t > fc and > 0 for t > j — k+1. 

We are now in a position to formulate results on the limiting behavior of dynamics Qtow [ct] and 
QL[ a\. In this section we prove the following: 

Theorem 8.7. Ts e —)> 0, the process R{t,e) of Definition 
process R{t). 

Theorem 8.8. As e —)• 0, the process L{t, e) of Definition 
process L{t). 

Corollary 8.9. The (possibly truncated) arrays R\{t) and ^-/L^j_]^j_i{t) for 1 <k < min{t, j} < n 
have the same distribution. 


8.5 eonverges in distribution to the 


8.4 converges in distribution to the 


In particular, and L\{t) have the same distribution. The latter fact was proven in 

and was used to analyze the strict-weak polymer partition function via the geometric RSK row 
insertion (see (8.2.1 below), and to establish the Tracy-Widom asymptotics for the strict-weak 
polymer. 

See also |58) for the close relation between the log-gamma and strict-weak polymers, where it 


is explained that one is the complement of the other at the level of lattice paths. To the best of 
our knowledge, the full statement of Corollary |8.9| has not previously appeared in the literature. 

Let us provide a brief outline of our proofs of Theorems |8.7| and |8.8| which are presented in the 


rest of this section. First, in (8.2 we describe the constructions of the geometric RSK correspon¬ 


dences, which will serve as e —)• 0 limits of elementary steps used in dynamics Qrowfa] and Q^Qj[a]. 
Then in §8.3| we prove a number of lemmas concerning e —)• 0 behavior of the (/-distributions from 
Finally, in §8.4| we use these ingredients to establish the desired statements. 


16.1 


8.2. Geometric RSKs. As we already know, the dynamics on (/-Whittaker processes con¬ 
structed in ^ degenerate for q = 0 into the dynamics Q'iTw [a] and based on the classical 

RSK row or column insertion, respectively. In this subsection we describe the corresponding 
geometric Robinson-Schensted-Knuth insertions, which will serve as building blocks for under¬ 
standing q 1 limits of the dynamics on (/-Whittaker processes. 

The q = 0 and q 1 pictures (i.e., the classical and the geometric RSK correspondences) are 
related via a certain procedure called detropicalization. Namely, the geometric RSK row insertion 
introduced in |44] is obtained by detropicalizing the classical RSK row insertion by replacing the 
(max, -|-) operations in its definition by (-I-, x). About the geometric RSK row insertion see also, 
e.g., p], [22], m, and [l9|. 

By analogy with the geometric RSK row insertion, one can define the geometric RSK column 
insertion, by detropicalizing the classical RSK column insertion, this time replacing the (min, -|-) 
operations by (-I-, x). 
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Remark 8.10 (Names and notation). The geometric RSK correspondences are also sometimes 
called tropical RSK correspondences HU, [53], [22], despite the fact that they come from the 
process of detropicalization. We adopt a convention of calling them the geometric RSK corre¬ 
spondences (following, e.g., |60|, |19) . [57ji. The latter name arises in connection with geometric 
crystals (see [l9| for more background). 

Note that the word “geometric” in the name of the geometric RSK correspondences should be 
distinguished from the same word in the names of the geometric g-PushTASEP and the geometric 


g-TASEP (described in )6.3 and 16.6 respectively). The former refers to detropicalization of the 
classical RSK correspondences, while the latter is attached to the g-geometric jump distribution. 

Below in this section, by A, i^,... we will denote vectors (words) with continuous components, 
and not signatures as before. To indicate the difference, we will use superscripts to denote their 
components. 

8.2.1. Geometric RSK row insertion. Consider a triangular array z^. (1 < A; < j < n) of nonneg¬ 
ative real numbers, such that a word Zk = {z ^,..., z^) either has all positive entries or is equal 
to (1,0,..., 0) (in which case we call it an empty word). 

First, define the geometric row insertion of a nonempty word a = (a^,..., a"") into a nonempty 
word A = (A^,..., A*^) as an operation that takes the pair {A, a} as input, and produces a pair of 
words {u = ,..., v"'), b = ..., 6”)} as output via the following rule: 


KK .. 

i=k 

tr = a-' 


At-Vi 


If A is an empty word, then by definition h is not produced, while 

:= (a^a^a^+^...,a^o^+^•••a^) 

is produced according to the same rule. The word b is also not produced for k = n. Observe that 
always i/t = (At -|- ioi k < j <n and = A^a^. 

a = a(0) 


z\- 


Z2 


t 

a(l) 


4 


a(2) : 


a{n — 1) 


Figure 30. Geometric RSK row insertion. 
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Definition 8.11. The geometric RSK row insertion of a word a = into an array 

^ 1 Cl r nr- 


zl is defined by consecutively modifying the words zi,..., Zn via the insertion according to the 


diagram on Fig. 30 The bottom output word a(l),a(2),... of each insertion is then used as a 
top input word for the next insertion. If after some insertion no bottom output word is produced, 
then no further insertions are performed. 




“1 


U.2 


d. 




Figure 31. Array and strip for the geometric RSK row insertion. 


The geometric RSK row insertion is related to the polymer partition functions of 18.1 in the 
following way: 

Proposition 8.12 ([53]). If we start with an array z of empty words, and consecutively insert 
into it nonempty fixed words ui,..., at, a* = (al,..., a"), via the geometric RSK row insertion, 
then in the obtained array we have 


zUt) = 


Ri{t){ai, 


for all t > k. 


Here with a slight abuse of notation we denote b y Rf {t){ai, ...,at) the same weighted sum over 
k-tuples of nonintersecting paths as in Definition\8.^ but in a strip in which each node {s,i) has 


a deterministic weight a\ (see Fig. 31). 


8.2.2. Geometric RSK column insertion. Consider a triangular array yj. {1 < k < j < n) of 
nonnegative real numbers, such that in each word = (y^,..., y^) either all entries are positive, 
or there is k < j < n, such that y^ = 1, y^ = 0 for j < i < n and y^ > 0 for A; < i < j. We again 
call (1,0,..., 0) an empty word. 

To define the geometric RSK column insertion first define the insertion of a word a = {a^, ..., a"') 
with positive entries into a word A = (A^,..., A"") as an operation that takes the pair {A, a} as 
input, and produces a pair of words { 
following rule: 


V = 


(u^...,u 

II 

+ 

a^X^ 


KK + AA-1 

for fc < j < n 



a-^— 

V 

o 

Xi-^v3 


a3v3~3^ 

if X3 = 0 and X^ 

aR 

if Al-i = 0. 


/-I 


> 0 , 
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Definition 8.13. The geometric RSK column insertion of a word into an array is defined similarly 


defined above, in order as on Fig. 30 


to the row insertion (Definition 8.11), by consecutively performing the column insertion operations 


. Vi 


. 2/2 






V2 


-.Vi 


. vKi 


• y; 


•Vn-l 


•— 



-• 

«1 

a| 

at 

w 



- w 

a\ 

• -•- 

al 

•-• 


Figure 32. Array and strip for the geometric RSK column insertion. 


Note that in this definition corresponds to in the classical RSK column insertion. 


We will need the following fact which is analogous to Proposition |8.12' 


Proposition 8.14. If we start with an array y of empty words, and consecutively insert into it 
words ai,... ,at with positive entries via the geometric RSK column insertion, then in the obtained 
array we have 


yiit) = 


Li_^{t){ai,... ,at) 


for all t > j — k + 1. 


Here again we denote by L;[(t)(oi ,... ,at) the same weighted sum over k-tuples of nonintersecting 


paths as in Definition 


weight a\ (see Fig. 32). 


8.t^ but in a strip in which each edge (s — l,i) —>■ {s,i) has a deterministic 


Proof. Our proof is similar to that of Proposition 8.12 (the latter is given in [53]). 

For a = (a^,... denote by H{a) the n x n matrix such that := a*, = 1, 

and other entries are 0. For a = ,... ,«"■), denote by Hk{a) the n x n matrix of the form 

~ ''' ’ A* > 0 for F < i < j and A* = 0 for j < i < n, 

denote by G'(A) the n x n matrix of the form 

0 
0 

0 Idn—j 

where G is the upper-triangular (j — k + 1) x {j — k + 1) matrix with 

Gp,r = ^p+k-2 ioT 1 < p < r < j - k + 1. 

Assume A^“^ = 1. 

The key to the proof is the commutation relation 

G{X)Hk{a) = Hk+i{b)G{i 2 ), 



(8.3) 
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whenever a pair of words v = ,... , 5"') is obtained by inserting a = 

{a^,..., aP) into A = (A^,..., A”). 


To check (8.3), denote its left-hand side by L and right-hand side by R. Clearly, = 1 

for 1 < i < A: — 1 and Lj^j = Rj^j = a* for j + 2 < i < n, and Lj+ij _|_2 = Rj-(-ij -|-2 = 1- Otherwise 
U,m = Ri,m = 0 unless k < i < j + 1 and k < m < j + \. Let us thus assume that the two latter 
inequalities hold. On the diagonal, for k < i < j we have 


. A* 


Lii = = Rii, \_kk = a^\^ = v'" = ^k 

’ iA~^ ’ ’ ’ 


and 


L 




. , w+i 
= = — = R 




■i+ij+i- 


Above the diagonal, for k < i < m < j + we have 


A 


m—1 


+ 




km - 


_ _1_ _ U_ _I __ — P. 

\i-l ^ lA-l lA 


since ^ + ^ = + 1) = wk, and hnally 


I — \m—1 , \mm _ _ p 

^k,m A -|-AQ. U 


This completes the proof of the commutation relation (8.3). 


By applying the commutation relation multiple times according to the geometric column RSK 


insertion (Definition 8.13), we get 

G{yn{t)) 


Giyiit)) = Hiai) 


Hiat). 


(8.4) 


Observe that the (i,j)-6atry of the right-hand side above is equal to the sum of weights of all 


directed strict-weak (as on Fig. 29 right) paths from (0,i) to where the weight of a path is 


given by the product of weights of horizontal edges, as before. Indeed, this entry is equal to 


E 

ii=i,it+i=j i=l 


^ ^ J_ _[—1 (^Aie=ie+i)- 

l<ii..,it+l<n:ii=i,it+i=j ^=1 


By the Lindstrom-Gessel-Viennot principle isa, m, the determinant of the minor of the right- 
hand side at the intersection of the first k rows, and columns from (j — k + l)-st to j-th, is 

Next, observe that for 1 < s < A:, j — A: -|- 1 < p < j, the (s,p)-entry of the left-hand side 
of ( |8.4| ) is equal to the sum of weights of directed up/right (as on Fig. 
from (A; -|- 1 — s, s) to (minj/c -|- A -|- 1 — p, A;},p) in the array as on Fig. 


29 


33 


left) lattice paths 
(the left picture if 

t > j, and the right one if j — A: -|- 1 < t < j). The weight of each path is dehned to be 
the product of weights of all nodes along the path. By Lindstrom-Gessel-Viennot principle, 
the determinant of the minor of the left-hand side of ( |8.4[ ) at the intersection of the hrst k 
rows, and columns from [j — k + l)-st to j-th, is equal to the sum of weights of all fe-tuples 
of nonintersecting paths from (A:, 1),..., (2, A — 1), (1, k) to (min{2A + t — j, k},j — A -|- 1),..., 

— 1), (min{A -|- A -|- 1 — j. A}, j). There is only one such tuple, which cov- 
k 


(min{A -|- A -|- 2 — A 
ers all points on Fig. 


33 


and has weight y™ 


i=l 
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yj 

vj-l 

y{ 

yi 

1 

y]-t 


vi~^ 

vi-\ 

yp^ 

yp^ 

yp\ 

y]^ 


1 • 

• 

• 

^ • 

• ' 

■■ • 

yP^ 







■ y\ 







• 


vt 


yi 

yp^ 

yt 


Vl 

yp^ 

k • 


• 

k • 


• 

- 1 

Sfe-l 

• 

vl 

v\ 

• 

k-i 

Sfc-1 

• 

yl 

y\ 

• 



y\ 



y\ 

1 


• 

1 


• 

1 

2 

k 

1 

2 

k 


t > j j — k+l<t<j 


Figure 33. 


Arrays used in the proof of Proposition 


l.U 


Therefore, 


i=l 


which establishes the desired statement. □ 

8.3. Asymptotics of (^-deformed Beta-binomial distributions. We will need several lem¬ 
mas about the limiting properties of the distributions \ y) (6.1). 

Lemma 8.15. Let he a hyo-valued random variable with 


Prob{X^ = j) = {a;q)c 


w 


for a = e and q = e 




Then as e —)• 0, eexpleX*^} converges in distribution to Gamma ^{0). 
Lemma 8.16. Let n(e) be a function M>o —> ^> 0 ? such that 


lime ^exp{—en(e)} = (p. 


e-S>0 


Let be a Z>Q-valued random variable with 

Prob{X^ = j) = I n{e)) for a = e~^^ and q = e“^ 

Then as e —)• 0, exp{—eX^} converges in distribution to ip + Gamma(0). 

These two lemmas were both proven in |25] (Lemma 2.1 and a part of proof of Theorem 1.4, 
respectively). In the next three lemmas, parameters of distributions which are not explicitly fixed 
are assumed to depend on e, and sometimes might also be random themselves. 
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Lemma 8.17. Fix C and 0 < cr < 1. Let be a "Lyo-valued random variable distributed 
according to | ri) with q = e~'^, Cq~"' cr, n > \oge~^—e~^C, and\ogri+2ne < logci. 

Then as e —^ 0, eV'^ —> log(l + a). 

Proof. The fact that with such parameters is indeed a probability distribution for e small 

enongh follows from ineqnalities in the statement of the lemma. Let A{e) := e“'^)oo- By [TJ 

Corollary 4.1.10], 

(e loge-i-Ce-iJ < A{e)C' 

for all e small enough and some constant C that depends only on C. As e —?■ 0, 

pr 

e\og{e~^-,e~^)^r/e\^ / log(l - e""^)dx, 

Jo 

since the left-hand side is a Riemann snm for the right-hand side integral, hence we have 

(e“";e“")oo 


elog- 


log(l — e ^)dx. 


(e e ^)oo Jq 

(Note that althongh this integral blows up at 0, it is still finite and convergence holds.) 
Fix 5 > 0. For e small enough, 

_ q~^)n-k (n\ ^ 


Prob{Y^ = k) = 


{r]',Q ^)r 


(e Le ")a 




< 


< C' exp log 2cr — ~ ® ^)dx^ 


(e Le ^)k 

< e~'^ for T large enough and k>T/e. 


Hence 

OO OO 

Prob{Y^ > T/e) < ^ ^ e"^^ 

A:=rr/e] i=0 

which can be made less than 5/2 for all e small enough by choosing sufficiently large T. 

Observe that for k < T/e and e small enough there is some constant Cq that depends only on 
C and T, such that 


^-1 ^ )k {q]q)n ^ 

° “ (^9; 9)00(7?; g“^)n (g;g)n-fc “ 


Let 


P' 1 }) ^'0 —logCJ 

/(V') := —+ (logcj)V’— / log(l — e“*)fix — / log(l — e“^)dx 

Jo Jo 

for fj >0. Then 

/'(V’) = -l-logfT — log(l — e“'^) — log(l — e“'^+*°®°'), 

which is strictly decreasing, and /'(log(l -1- cr)) = 0. Hence / attains a unique maximum at 
log(l -|- a), so one can choose Mi > M 2 > M 3 > M 4 such that 

/('i/’) > Ml for i/ E (log(l -I- cr) - 5/2, log(l -f cr) -|- 5/2) 

and 

/(if) < M 4 for ^ (log(l + a)-5, log(l -f cr) -f 5), 
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and for e small enough 

Prob{eY^ G (log(l + cr) — 6/2, log(l + a) + 6/2)) > 

and 

( T — 26 \ 

-h 2 j yl(e)e^^/^ 

Therefore, for e small enough 

Prob{eY'^ G (log(l + cr) — (5, log(l + cr) + (5)) > 1 — (5, 
and this completes the proof. □ 


Lemma 8.18. Fix C and 0 < cr < 1. Let Y^ be a X>o-valued random variable with 

Prob{Y^ = j) = I n) 

for a^a as e^0,q = e~^ and n > e~^ log e~^ — €~^C. Then as e —)• 0, eY’^ —)■ — log(l — a). 
Proof. Suppose ct > 0 and fix <5 > 0. We can write 


ProbiY- = fe) = ^ 


a 


(e ^;e ^)k (e ^;e ^)n-k (e ^;e ^)k 

< exp ^-((^logcr)T — J log(l — e~^)dx)'^ < 
for T large enough and e small enough, such that k > T/e. Hence 

OO 

Prob{Y^ > T/e) < ^ 


i=0 


which is less than 6/2 for T large enough. 

For e small enough, k < T/e, and some constant Cq that depends only on C, T, and cr, one 
can write 


^ (a;e ")n-fc(e ";e "), 


Cq {a; e ^)oo < 


Let 


(e e ^)n—k 
fp 


< Co{a]e ^)oo- 


/(V^) := (logcr)V’- / log(l-e ^)dx 
Jo 

for fj > 0. Then 

f'fif) = logcr - log(l - e~'^), 

which is strictly decreasing, and /'(—log(l — a)) = 0. Hence / attains a unique maximum at 
— log(l — fj), so one can choose Mi > M 2 > M 3 > M 4 such that 

f{'ip) > Ml for V' G (— log(l — cr) — 6/2, — log(l — cr) + 6/2) 

and 

f{'if) < M 4 for f) i{- log(l - a)-6,- log(l - a) + 6), 
and for e small enough 

Prob{eY'' G (- log(l - a)- 6/2, - log(l - a) + 6/2)) > - 1^ (a; e~^)^e^'^l^ 



g-RANDOMIZED RSK CORRESPONDENCES AND RANDOM POLYMERS 


72 


and 

( T — 2S \ 

— - -h 2 j (a;e“^)ooe^®/^ 

Therefore, for e small enough we can write 

Prob{eY^ G (— log(l — a) — S, — log(l — cr) + (5)) > 1 — (5, 


and this completes the proof for the case a > 0. 

If (7 = 0, fix arbitrary u > 0. For all large enough U, such that Uu > ^Uu — log(l — e~^)dx, 
€ small enough, and k > j, 

ProbCV^ = k) < -— < exp ( -(—like — / 


log(l — e 


'^)dx) ) < e 2 e . 


So, 

OO 

Prob{eY'^ > u) < e~^ ^-Ui/ 2 ^ 

i=0 

which can be made less than any given d > 0 for U large enough. Thus, we have convergence 
eY^ — 7 - 0 in distribution. □ 


Lemma 8.19. Fix a,a £ (0,1) such that a(l + cr) < 1. Let Y^ be a Z>Q-valued random variable 
with Prob{Y^ = j) = Pq-^,a,o{j I log(l + cr) as e ^ 0, q = e~^. Then as e —)• 0 , 

eT*^ —)■ log(l + aa). 

Proof. a{l + a) < 1 ensures that this is indeed a probability distribution for e small enough. This 
distribution looks as 

Prob{Y^ = k)= . - {aq^-^)\a-q-X-k- 

[T Q)k{Q', Q)n-k 

Let 

ri) ^log(l+CT)—■(/' 

f{'ip)'.= — / log(l — e“^)dx — / log(l — e“^)(ix + V'logo —' 0 ^ + V'log(l + cr) 

Jo Jo 

/ — log a 

log(l — e~^)dx for log(l + cr) > 0 > 0 . 

- log a—log(l+(j)+'j/; 

Then 

gb 

/'(0) = - log(l - e~^) + log(l - y-—) + log a - 20 + log(l + a) - log(l - q ;(1 + cr)e"’^), 

i + cr 

which is strictly decreasing, and /01og(l + aa)) = 0. Hence / attains a unique maximum at 
log(l + aa), so one can choose Mi > M 2 such that 

/( 0 ) > Ml for 0 G (log(l + aa) — 5j2, log(l + acr) + 5/2) 

and 


/(0) < M 2 for 0 (log(l + acr) — 5, log(l + acr) + 5). 
Hence for e small enough, 

Prob{eY^ G (log(l + aa) - 5/2, log(l + acr) + 5/2)) > - 1^ (e“0 
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and 


Prob{eY^ ^ (log(l + aa) — 5, log(l + aa) + <5)) < 


/log(l + a) -26 


+ 2 




□ 


Thus, for e small enough, we have 

Prob{eY^ G (log(l + aa) — 6, log(l + aa) + (5)) > 1 — 5, 
and this completes the proof. 

8.4. Proofs of Theorems 8.7| and |8.8| , In our proofs, we denote by 

Aie),B{e),C{e),D{e),Eie),Fie) 

possibly random positive valued functions tending to deterministic constants A,B,C,D,E,F, 
respectively, as e —)• 0. This notation will be repeatedly used in the arguments below for defining 
conditional probabilities. 


8.4.1. Proof of Theorem 8.1. We must show that for fixed {t,k,j), such that k < min{t,j} and 
j < n, the random variable e) conditioned on Rj^{T, e) —)• X'^{T) for all (T, K, J) < (t, k,j) 
in the lexicographic orderconverges as e —)■ 0 to Rl.{t) conditioned on Rj^{T) = X'^(T) for all 
(r, K, J) < {t, k,j) in the lexicographic order. Here X'^{T) are some fixed constants. In the rest 
of the proof we will always assume this conditioning. 

Right edge {k = 1). The Markovian projection of the Qrowfct] dynamics on the right edge is the 
geometric g-PushTASEP ( ^6.3[ ), hence the proof of the theorem for the right edge is the same as 
showing that suitably rescaled positions of particles in the geometric g-PushTASEP converge to 
the partition functions of the log-Gamma polymer (Definition |8.2[ ). 

a) If t = 1, then rj^i(l,e) = rj_i_i(l,e) + an independent random moveme nt d distributed 

according to I (assume ro,i(l,e) = 0 and Aj’(l) = 1 ). By Lemma 8.15 

log(R{(l, e)) = rj,i(l, e)e - j loge“^ 

= rj-i,i{l,e)e - {j - l)loge“^ + de - loge“^ 

^log(Ari(l)) + log(r) 

for an independent random variable P = distributed according to Gamma”^(0j + 9i), which 
is consistent with R{( 1 ) = x(~^{l)a{. 

b) If j = 1 and t > 2, then ri^i(t,e) = ri^i{t — l,e) + an independent random movement d 

distributed according to Tq,aiat,oi^ I Lemma 8.15 

log(.R}(t,e)) = ri,i(t,e)e - flog 

= ri,i(t - 1 , e)e - {t - 1) loge”^ + de - loge”^ 

^log(Ai(t-l)) + log(r) 

for an independent random variable T = distributed according to Gamma”^(0i + 9t), which is 
consistent with R\{t) = X^(t — l)aj:. 

c) Assume t > 2 and j > 2 . Condition on 

rjj(t -l,e) = (t + j - 2)e~^loge~ 


+ e ^logA(e), so A;( (t - 1) = A; 


30. 


That is, T < t, or T = t and K < k, or T = t, K = k and J < j. 
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- l,e) = (t + j - 3)e ^loge ^ + e Mog5(e), so Xf \t - 1) = B; 

e) = (t + j - 2)e~^ loge~^ + e~^logC(e), so x(~^(t) = C. 

The movement of the rightmost particle on the j-th level during the time step t — 1 —)■ t which 
happens due to the pushing by the rightmost particle at level j — 1 behaves as log(l + (by 


Lemma 8.17). The independent m ovem ent of the rightmost particle on the j-th level behaves as 

e“^logr + e“^loge“^ (by Lemma 8.15), for an independent random variable T = distributed 
,- 1 / 


according to Gamma + 0i). Therefore, 


C. 


log(Rl(t,e)) log^ + log(l + — ) TlogT = log((A + C')r), 

which is consistent with R{(t) = (X((t — 1) + 

/c-th edge from the right for k > 2. 

a) Assume t = k. We have rj^k{k,e) = rj_i fc(/c,e) + a movement due to pulling from the 
{k — l)-st particle from the right on the (j — l)-st level (assume rk-i^k{k, e) = 0 and X^~^{k) = 1 ). 
Condition on 

rj_i,fc_i(A: -l,e) = {j -k + l)e"Moge"^ + e"^logL»(e), so Xlz\{k - 1 ) = L»; 

rj_i,fc_i(/c,e) = (j - fc + 2 )e“Moge“^ + logL;(e), so Xlz\{k) = E] 

rj^k-i{k - l,e) = (i - A; + 2)e"Moge"^ + e"MogF(e), so Xl_^{k -1) = F. 

By Lemma 8.17, this movement times e and minus log e~^ converges to log{E/D) — log(l + E/F). 
Therefore, 


log(i7^(fc, e)) = rj^k{k, e)e - {j - k + 1) log 


-1 


log{Xl-\k)) + log{E/D) - log(l + E/F) 
7-1 


= log 


Xi-\k)EF \ 

{F + E)D J ’ 

Xl-'{k)XC]{k)Xl_,(k-l) 

XLi(k-l) + xclik)) xi-_l(k - 1) 


which is consistent with 

Riik) = 


Indeed, if we insert (via the geometric row insertion) a nonempty word b = ( 6 ^,..., 6 ”) into 
the empty word Xk = ( 1 , 0 ,... , 0 ), where b is itself the bottom output of the insertion of a = 
into Afe_i = (A^l},..., = (Xj/zlik - 1),.. ■, XJ/_-^{k - 1)), then we get 


A = vi = ui 


= 


Atl-Li 

which is the same as the expression for Rl,{k) above, 
b) Assume t > A: + 1 and k < j < n. Condition on 


j—i 

4-1 


K--\ -n + K-4 


rj^kit-l,e) = {t + j-2k)e loge +€ logA(e), so A^(t - 1 ) = A; 
rj_i,fc(A - l,e) = (t + j - 2A: - l)e"^log€"^ + e"MogB(e) 




so 
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= {t + j - 2k)e~^loge~^ + log C'(e), so = C\ 

-l,e) = {t + j - 2k + l)e"^loge"^ + e"MogZ?(e), so Xlz\{t - 1) = D] 
= {t + j - 2k + 2)e"Moge"^ + e"Mog^(e), so Xlz\{t) = E- 
-1+) = {t + j - 2k + 2)e"^loge"^ + e"MogF(e), so Xl_^{t -l) = F. 


Movement of Vj^k during time step t — 1 —)■ t due to pushing by the A:-th particle from the right 


on the (j — l)-st level behaves as e ^log(l + (by Lemma 8.17). The movement due to the 


pulling (by the [k — l)-st particle from the right on the (j — l)-st level) times e ^ minus loge 
by the same lemma tends to log(^) — log(l + %). Hence we may write 




\og{Rl{t, e)) logH + log ( 1 + ^ ) + log 


- log ( 1 + ^ ) = log 


{A + C)EF\ 
{F + E)D ) 


which is consistent with 


R+) = 


[xiit - 1 ) + xi 


xir+f]xi_,(t - 1 ) 


+ H~-\(i))+r+t-i) 


•7-1 / 


Indeed, if we insert (via the geometric row insertion) a nonempty word b = (6 ^,..., 6 ”) into a 
nonempty word A*, = (A^,..., A^) = (X^(t-l),.. ■,X^{t - 1 )), where b is itself the bottom output 
of the insertion of a = {a^~^,. ..,a+) into A^-i = (A^lJ,..., A^_;^) = {X^zlit - 1 ),... - 

1 )), then we get 




, j-i u/i-i 
—pn —-+k +^fcJTpT 


j—1 

'4-1 


Xi 


k-l^k-l 


K-i -ri+At; 


which is the same as the expression for Rl,{t) above. 


c) Finally, if j = A; and t > A; + 1, then the previous argument carries out with the exception 
that in this case the leftmost particle on the j-th level experiences only pulling of the leftmost 
particle on the {j — l)-st level, hence we should take C = 0 in the formulas from part b). 


This completes the proof of Theorem 8.7 


8.4.2. Proof of Theorem \8.^ We must show that for fixed {t, k, j), such that A:<j<t + A; — 1 and 
j < n, the random variable Lj^{t, e), conditioned on Lj,{T, e) ^ YJ{T) for all (T, K, J) < {t, k, j) 
in the lexicographic order, converges as e —)• 0 to L'(.(t), conditioned on LUT) = Y^{T) for 
(r, K, J) < {t, k,j) in the lexicographic order (here Y^{T) are some fixed constants). In the rest 
of the proof we will always assume this conditioning. 


Left edge {k = 1). The Markovian projection of the Qcoil*^] dynamics on the left edge is the 
geometric (;-TASEP ( ^6.6[ ), hence the proof of the theorem for the left edge is the same as showing 
that suitably rescaled positions of particles in the geometric g-TASEP converge to the partition 
functions of the strict-weak polymer (Definition |8.3[ ). This was already done in [25], but we still 
include this part in the proof for the reader’s convenience. 


a) lit = j = 1, then by Lemma 8.15, log(Lj(l, e)) = loge ^ —^ 1 , 1 ( 1 , e)e converges in distribution 


to logT for a random variable T = aj distributed according to Gamma(0i -|- 0i). 
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b) Assume t = j > 1. Then m = e) is distributed according to I f))- 

If we condition on 

- l,e) = e“Moge“^ - e“MogF(e), so - 1) = F, 

then by Lemma 8.16| 

log(L{(j, e)) = loge"^ - 4i(j, e)e ^ log(F + T) 

for an independent random variable T = a^- distributed according to Gamma(0j + 9j), which is 


consistent with L\{j) = a] + Yl (j-1). 

c) Let t > j = 1. By Lemma 8.15[ the quantities 

log(Ll(Le)) - log(L}(t - l,e)) = loge"^ - - h^i{t - l,e))e 

converge in distribution to log L for a random variable L = distributed according to Gamma(0i+ 
6t), which is consistent with L\{t) = alY^t — 1). 

d) Assume t > j > 1. Condition on 

4i(t - l,e) = {t - j)e~^ loge~^ - e~^logA{e), so Y^{t - 1) = A; 

-l,e) = {t - j + l)e"^loge"^ - e~^logF{e), so Y^~^{t -1) = F. 

By Lemma 8.16[ the movement of the leftmost particle on the j-th level during the time step 
t — 1 —)• t times e and minus log e~^ converges to — log(^ + r) for an independent random variable 
r = Uj distributed according to Gamma(0j + 6t). Therefore, 

log(Li(L e)) ^ log A + log(F/A + T) = log(F + AT), 
which is consistent with L{{t) = Y^~^{t — 1) + alY^^t — 1). 

Second edge from the left {k = 2). 

a) If j = 2, t = 1, then we have to look at f 2 , 2 (l, e) = -^1,1(1) e) + a jump m distributed according 

to ¥’g,a 2 ai,o("^ I Lemma 8A5 

log(L 2 (l)) - log(L}(l)) = loge“^ - me logT, 

where T = af is an independent random variable distributed according to Gamma(02 + ^i); which 
is consistent with -£ 2 ( 1 ) ~ alY^il). 

b) Assume j > 2, t = j — 1. We have fj, 2 (j ~ 1t) = + 1 ^ 2 , where mi is an independent move 

distributed according to I ^j-i, 2 {j — 2,e)), and m 2 is the push from the leftmost 

particle on the {j — l)-st level distributed according to 


<3 

Condition on 




„b-i.20-2.^) - 2 , e) - mi - m2 I - 2 , e) - mi). 


- l,e) = e ^loge ^-e MogF(e), so T/ (j - 1) = F; 
ij-ifiU - 2,e) = 2e“Moge“^ - e“MogF(e), so Y^~^{j -2) = F. 


By Lemma 


1.16 


loge ^ — mie —?■ logT, where T = a ]_i i s an independent random variable 
distributed according to Gamma(0j + Oj-i)- By Lemma 8.17 


FT, 


log(L^(j - l,e)) = 21oge - ij, 2 {j - l,e)e logF + log(l + —) = log(F + FT) 
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which is consistent with Y^{j — 1) = ^{j — 2) + a^j_iY^ ^{j — 1). 

c) Let j = 2, t > 1. Then 

e) = i 2 , 2 {t - 1, e) + e) - - 1, e) + m, 

where the jump m is distributed according to 

Condition on 


By Lemma 


^ 2 , 2 {t - 1, e) = loge“^ - e~^ log A(e), so Y^it - 1) = A] 

- l,e) = (t - l)e“Moge“^ - e“^logB(e), so - 1) = B; 
= te~^\oge~^ - e“^logC'(e), so Yl{t) = C; 

^2,i(Ce) = {t- l)e“^loge“^ - e“MogL;(e), so if (t) = E. 

8.18, me —)■ — log (l — , hence 


log(L^(t, e)) = i 2 , 2 {t, e)e - (t + 1) log e ^ ^ log ( ^ ( 1 - ^ 


which is consistent with L|(t) = if (t — 1) ^1 — I'^deed, we have 

_ f2/. 1 _ f 2/. 1 \“ 1)-^i(^ “ 


^2{t) — -^2(^ “ 




= m-iy- 


Lf(t-l) L\{t-l)Ll{t) 


= Ll{t-l) 1- 


Li(t-1)\ L\{t) 


Llit) 


d) Assume j > 2, t > j — 1. Condition on 

^j, 2 {t - l,e) = {t - j + 2)e“Moge“^ - e“MogA(e), so Y^{t - 1) = A; 

- l,e) = (t - j + l)e“Moge“^ - e“MogB(e), so if “ft - 1) = B; 

£j_i,i(t,e) = {t- j + 2)e~^ loge“^ - e~^ logC'(e), so if “ft) = C; 

f,i(t,e) = (t - j + l)e“Moge“^ - e“^logL;(e), so if (t) = E; 

f-i, 2 (t - l,e) = (t - j + 3)e“Moge“^ - e“MogF(e), so if “ft - 1) = F. 

Denote by m the independent move of the particle which is the second from the left on the j-th 
level. This move is distributed according to 


a,!(*■■=) ,o("^ I f-L2f “ C f “ f.sf “ 1> f)' 


As in the previous case, by Lemma 8.18, me —>■ — log (l — 

Thus, we see that t'j, 2 (t, e) = £j-i^ 2 {t — 1, e) — M, where M is distributed according to 

I f-i, 2 (t - l,e) - ij^ 2 {t - l,e) - m). 
(l — f gA + 1. Therefore, 


Hence by Lemma 


8.17 


Me 


log(Lft,e)) =ij^ 2 {t,e) -{t-j + 3)e Moge 


-1 


log 
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which is consistent with ^ 2 ^^) 
checks that 


1 - 


Yl(t) 


- 1) + Yi Indeed, one 












+ 1/2 {t — 1) 


L{ ^(^-1)^ 'L{ \t) 


i-l, 




^i(^) J L{ ^(t- 1) 


+ LP(t-l). 


k-th edge from the left for k > 2. 


a) Start with assuming that j - 
distributed according to ^q^a,^ai,o 


= k, t = 1. We have_^fc(l,e) 
(m I oo). By Lemma 


8.15 


4-i,fc-i(l,e) + a jump m 


log(Lfc(l)) - log(L^_J(l)) = loge ^ - me ^ logL, 


where F = is an independent random variable distributed according to Gamma(4 + 4); which 
is consistent with ^^(1) = a\Yl^Zi{^)- 

b) Let j>k,t = j — k + 1. We have — k + l,e) = mi + m 2 , where mi is an independent 
move distributed according to +’^, 0 ^ 0 : -fe+i,o(™'i I — k,e)), and m 2 is the push from the 

(k — l)-st particle from the left on the {j — l)-st level distributed according to 

- fc,e) - mi - m 2 I 4-1,fcO' -k,e)- mi). 

Condition on 


4_i,fc_i0'-/i; +l,e) = (A; - l)e Moge ^-e 4ogL;(e), so I^'-i (j - A: + 1) = 
4 - 1 ,fcO’ -k,e) = ke~^loge~^ - e"4ogF(e), so -k) = F. 


By Lemma 


8.15 


loge ^ — mie —>■ logF, where F = is an independent random variable 

distributed according to Gamma(4 + &j-k+i)- By Lemma 8.17 


log(L^(j - A: + l,e)) = A:loge -4+(j - A: + l,e)e logF +log 


EY \ 

1 + —j =log(F + EF), 


which is consistent with Y^{j — k + 1) = Yj? ^{j — k) + 
c) Assume j = k, t > k. Condition on (for 1 < i < A: — 1) 

4-1,j(A - l,e) = (A - A: + 2i - l)e“^ loge"^ - e“4ogBi(e), so Y^~^{t - 1) = Bf, 
4-1,4A+) = (a - A: + 2A)e"4oge"^ - e'^ogC^e), so Y^~^{t) = 4; 

4,j(A - l,e) = (A - A: + 2A - 2)e“^ loge“^ - e~^ log A(e), so Y^{t - 1) = A; 
4,j(A+) = (A - A: + 2i - l)e“4oge“^ - e~^ logEi{e), so Y^{t) = 4; 
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By Lemma 8.18, the independent move of the rightmost particle on the fc-th level times e converges 


to 0, while the push from the previous particles times e and minus loge ^ converges to 


-log 


nt~2 An-=ia 



which is consistent with 

Ykit) 

For t = k we take Di = 1. 


LUt-i) 


n-Zi D^ n ■=! g 

nti' E. nti' 


El — Bi 
Di 


d) Let j = k, k > t > 1. Make the same conditioning as in the previous part, but with different 
ranges of indices: k — t + l<i<k — 1 for Di, Ei, and k — t<i<k — 1 for Bi, Ci. Take 
Bk-t = Dk-t+i = 1- The independent move of the rightmost particle on the k-th. level times e 
converges to 0, while the push from the previous particles times e and minus loge“^ converges to 


which is consistent with 


-log 


[ntlt+lE^ntltB,''^) 


YkHt) 


LUt-l) 


n ■j.Ui A n g. 


{Ek-t+i 


e) Let j > k, t > j. Condition on 


Bk-t+i)- 


ij-i^i{t - l,e) = {t - j+ 2i - l)e Moge ^-e ^logi?j(e), so T/ - 1) = for 1 < i < A;; 
ij-i^i{t,e) = {t- j + 2i)e“Moge“^ - logCi(e), so = Q for 1 < i < k - 1; 

ij,i{t — l,e) = {t — j + 2i — 2)e~^ loge“^ — logHi(e), so [t — 1) = Di for 1 < i < k] 
ij^i{t, e) = {t — j + 2i — l)e“^ log log Ei{e), so Y^{t) = Ei for 1 < i < k — 1. 


The independent move of the A:-th particle from the left on the j-th level times e converges to 0. 
Denote by mi the distance from this particle to — 1, e) after the push from the {k — l)-st 

particle from the left on the {j — l)-st level. Let m 2 be this distance after pushes from other 


particles. By Lemma 


8.17 


mie —>■ log (1 + 


Ck-lDk 


By Lemma 


8.19 


m 2 e —>• log 



o-iPtOK Antic. 

A-iAnSantiA 



This is consistent with 

Liit) 


Bk + Yi{t-1) 


nti D'^nti C, El - Bi 

n-E^E^U'E^B, Di 


f) Finally, assume that j > k, j >t>j — k + 1. Make the same conditioning as in the previous 
part, but with different ranges of indices: j — t + 1 < i < k for Di, Dj-t+i = 1, j — t + 1 <i< k — 1 
for Ei, j — t < i < k for Bi, Bj-t = 1, and j — t<i<k — 1 for Ci. The independent move of the 
fc-th particle from the left on the j-th level times e converges to 0. Denote by mi the distance 
from this particle to ij-i^k{t — 1t) after the push from the {k — l)-st particle from the left on 
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the (j — l)-st level. Let m 2 be this distan ce after pushes from other particles. By Lemma 8.17 

. By lemma 8.19 


Ck-iDk A Wi=j-t Ci 


k-2 


77126 


log 1 + 




rfc-l 


rfc-2 


This is consistent with 


n A nfc-t+i a 


m = Bk+nit - 

1 L=j_t+i 1 L=j_t+i 

This completes the proof of Theorem |8.8[ 


B, 


■{Ej-t+i - Bj-t+i). 
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